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Foreword 

The  importance  of  watershed  ecosystem  research  to  the  management  of 
natural  areas  such  as  national  parks  cannot  be  overstated.  Watershed  conditions 
are  affected  by  uncertain  and  complex  interactive  environmental  trends,  many 
of  global  occurrence.  Questions  about  watershed  science,  although  often  site 
specific,  thus  require  answers  in  a  global  context.  Frequently,  input  to  the  highest 
levels  of  government  decision  making  is  needed  if  modem  society  is  to  mitigate, 
stop,  or  reverse  the  large  scale  alterations  of  watershed  ecosystems  world- 
wide. However,  our  ability  to  sensibly  and  effectively  manage  natural  resources 
in  today's  global  environment  is  often  constrained  by  our  lack  of  knowledge 
about  the  hydrologic  cycle  and  its  relation  with  the  geosphere  and  biosphere. 
We  know  that  the  current  condition  of  natural  resources  in  parks  and  other 
preserved  lands  is  subject  to  many  widespread  anthropogenic  changes  from 
acidification,  eutrophication,  toxic  substances,  boundary  encroachments, 
overuse,  species  shifts  or  extirpation,  desertification,  loss  of  biodiversity,  land 
use  change,  and  sea  level  rise.  Long-term  watershed  research  and  monitoring 
of  natural  and  remote  areas  in  the  National  Park  System  and  in  similar  reserves 
provide  important  data  on  ecosystem  processes  and  interactions  for  detecting 
spatial  and  temporal  changes  in  environmental  conditions.  These  data 
collections  allow  the  partitioning  of  cause-and-effect  relations  of  ecological 
change  in  watersheds.  They  also  serve  to  meet  reference  and  early  warning 
objectives  correlative  with  natural  ecosystem  change.  Accordingly,  the  scientific 
and  public  policy  communities  can  employ  watershed  ecosystem  information 
as  one  means  of  obtaining  early  indications  of  the  potential  effects  of 
anthropogenic  stress  and  an  improved  assessment  of  its  magnitude.  Use  of  the 
concept  to  address  several  goals  of  research  into  acid  precipitation  or 
biogeochemical  cycling  demonstrated  the  utility  of  these  integrated  watershed 
data  for  inter-ecosystem  comparison  between  preserved  and  other  watersheds. 
The  long-term  management  strategies  of  the  National  Park  Service  combined 
with  the  protected  nature  of  park  lands  placed  the  service  in  a  unique  position 
among  federal  agencies  and  among  international  conservation  organizations 
to  document  relations  between  ecosystem  effects  and  anthropogenic  influences. 

Studies  of  watershed  ecosystems  collect  long-term  baseline  data  on  the 
ecosystem  health  of  park  and  equivalent  reserves.  Research  and  monitoring  at 
the  watershed  level  by  the  U.S.  Geological  Survey  and  the  National  Park  Service 
since  1980  have  contributed  to  the  accumulation  of  important  baseline 
information  on  deposition,  meteorology,  hydrology,  ecosystem  functioning, 
and  biology  in  selected  parks  and  biosphere  reserves.  Important  information 
on  biological  diversity  and  biogeochemical  processes  have  also  been  obtained. 
Activities  have  ranged  from  needs  identification  to  reconnaissance  or  synoptic 
analyses  to  long-term  monitoring  and  long-term  ecosystem  research. 
Quantification  of  the  hydrologic  cycle  and  chemical  flux  is  a  major  objective 
of  the  watershed  program.  Such  measurements,  when  combined  with  other 


geographic  resource  data  (e.g.,  geology,  landuse,  topography,  and  historic  and 
pre-historic  records),  provide  a  better  understanding  of  ecosystem-level 
processes  and  of  how  watershed  ecosystems  respond  to  various  natural  and 
anthropogenic  stimuli.  Currently,  these  data  are  available  for  developing,  testing, 
and  implementing  state-of-the-art  methods  and  procedures  for  improved 
management  of  water  and  land  resources  at  the  national  and  international  levels. 

By  1 982,  the  National  Park  Service  had  implemented  studies  in  three  national 
parks:  Sequoia-Kings  Canyon  (California),  Rocky  Mountain  (Colorado),  and 
Isle  Royale  (Michigan)  national  parks.  Sites  were  selected  to  be  biogeographically 
representative  of  mature  ecosystems  in  relatively  remote  areas  that  receive 
varying  levels  of  atmospheric  contaminant  inputs.  Two  years  later,  the  largely 
pollution-free  Olympic  National  Park  (Washington)  was  added.  The  sites  have 
several  important  common  attributes  including  a  long-term  protected  status,  a 
core  area  (undisturbed  watershed),  research  hypotheses  and  experimental  design 
for  detection  of  ecosystem  change  and  understanding  the  mechanisms  of  change, 
basic  inventory  information,  research  or  monitoring  capabilities,  monitoring 
techniques  that  will  not  change  unless  calibrated  to  new  techniques,  sampling 
protocols  that  fit  within  the  time  frame  of  known  physical  and  biologic  events, 
and  standardized  protocols  for  data  collection,  sample  storage,  sample  archival, 
and  data  management. 

Today,  comprehensive  watershed  ecosystem  inventory,  monitoring,  and 
research  data,  a  series  of  publications  on  ecosystem  function  and  structure, 
and  to  some  extent  management  of  park  resources  that  span  more  than  15 
years  are  available  from  the  four  original  sites  and  from  Shenandoah  National 
Park  (Virginia)  where  watershed  studies  began  in  1979.  Shorter  or  discontinuous 
records  of  watershed  ecosystems  are  available  from  Big  Bend  (Texas),  Denali 
(Alaska),  Glacier  (Montana),  and  Great  Smoky  Mountains  (North  Carolina 
and  Tennessee)  national  parks  and  from  Noatak  National  Preserve  (Alaska). 
Six  parks  are  now  part  of  the  National  Park  Service-U.  S.  Geological  Survey 
Global  Climate  Change  Program:  Rocky  Mountain,  Isle  Royale,  Olympic, 
Sequoia,  Glacier,  and  Crater  Lake  national  parks.  Three  watershed  sites  are 
included  in  the  cooperative  Inventory  and  Monitoring  Program  of  the  service 
and  the  survey:  Shenandoah,  Great  Smoky  Mountains,  and  Denali  national 
parks.  Sites  range  from  the  hot  Chihuahuan  Desert  in  the  Southwest  to  the 
moist  boreal  forests  of  Michigan,  the  eastern  deciduous  forests  of  Virginia, 
Tennessee,  and  North  Carolina,  the  alpine  environment  of  California  and 
Colorado,  and  the  Alaskan  taiga-tundra. 

The  condition  of  natural  resources  within  these  natural  watersheds  has  often 
been  difficult  to  ascertain,  but  gaining  this  knowledge  has  provided  the  key  to 
understanding  the  nature  of  ecosystem  change.  The  existence  of  sites  with  a 
commitment  to  gathering  "long-term"  ecosystem  level  data  permits  research 
activities  aimed  at  testing  hypotheses  relevant  to  ecosystem  processes  and 
structure.  These  data  further  make  it  possible  to  question  existing  paradigms 
and  to  obtain  new  understandings  about  fundamental  relations  within  and 


between  naturally  functioning  ecosystems.  This  applies  equally  to  other  areas 
and  is  suggested  as  appropriate  for  a  potential  network  of  long-term  global 
baseline  research  sites.  Potentially  significant  conclusions,  common  to  most 
long-term  watershed  ecosystem  programs  and  studies,  are  thus  emerging. 

This  monograph  presents  the  results  of  the  first  13  years  of  watershed  research 
in  Olympic  National  Park.  It  is  an  example  of  the  value  of  staying  the  course  to 
increased  understanding  of  park  natural  resources  and  of  mature  boreal 
ecosystems. 

R.  Herrmann 
Watershed  Research 
U.S.  Geological  Survey 
Biological  Resources  Division 
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Abstract.  Since  the  1970s,  there  has  been  widespread  concern  about  the  effects  of 
anthropogenic  inputs  of  nitrogen  and  sulfur  on  forest  ecosystems  and  their  associated 
streams  and  lakes  in  North  America  and  Europe.  Few  studies,  however,  have  been 
conducted  in  areas  that  are  still  relatively  pristine.  Long-term  monitoring  in  these  eco- 
systems can  be  used  to  detect  global  change.  This  monograph  summarizes  our  current 
understanding  of  the  structure  and  functioning  of  two  small,  temperate  old-growth 
watersheds  in  the  Hoh  River  valley  in  Olympic  National  Park,  Washington,  in  an  envi- 
ronment with  little  apparent  anthropogenic  disturbance.  Using  the  small  watershed 
approach,  we  examined  vegetation  patterns,  biomass  distributions,  hydrology,  solu- 
tion chemistry  in  precipitation,  throughfall,  stemflow,  soil  and  stream,  and  long-term 
trends  in  precipitation  and  stream  chemistry  from  1984  to  1993.  Some  data  from  the 
1994-97  period  are  also  discussed,  particularly  NO3-N  in  precipitation  and  the  stream 
in  the  West  Twin  Creek  watershed.  The  primary  site  for  our  studies  was  the  valley-wall 
58-ha  West  Twin  Creek  watershed  (180-850  m  elevation),  located  32  km  from  the 
Pacific  Ocean.  The  second  site  was  the  125-ha  Hoh  Lake  watershed  (1250-1525  m 
elevation)  located  50  km  from  the  Pacific  Ocean.  The  Hoh  River  valley  contains  a 
wide  variety  of  species  and  community  types.  Five  plant  community  types  were  iden- 
tified in  the  West  Twin  Creek  watershed,  and  three  forest  community  types  and  a  wet 
meadow  community  were  identified  in  the  Hoh  Lake  watershed.  Elevation  (mostly 
reflecting  differences  in  temperature)  and  slope  differences  are  the  most  important 
factors  that  control  vegetation  patterns.  Western  hemlock  (Tsuga  heterophylla)  and 
Pacific  silver  fir  {Abies  amabilis)  were  the  dominant  tree  species  in  the  West  Twin 
Creek  watershed,  although  Douglas-fir  {Pseudotsuga  menziesii),  western  redcedar 
{Thuja  plicata)  and  Sitka  spruce  {Picea  sitchensis)  were  present.  Mountain  hemlock 
{Tsuga  mertensiana),  Pacific  silver  fir,  and  Alaska  yellow  cedar  {Chamaecyparis 
nootkatensis)  were  the  dominant  tree  species  in  the  Hoh  Lake  watershed.  Large  bio- 
mass accumulations  were  recorded  in  both  the  West  Twin  Creek  (average  live  tree 
biomass  1044.7  Mg/ha)  and  Hoh  Lake  (average  live  tree  biomass  1006.8  Mg/ha)  wa- 
tersheds. Coarse  woody  debris,  both  snags  and  logs,  is  a  dominant  structural  component 
in  these  watersheds.  Annual  precipitation  from  1984  to  1993  varied  from  267  to  475 
cm  at  the  Hoh  Ranger  Station.  Annual  stream  discharge  in  the  West  Twin  Creek  water- 
shed ranged  from  54%  to  89%  of  precipitation.  Precipitation  chemistry  was  dominated 
by  the  sea  salt  ions,  Na,  and  CI,  whereas  stream  chemistry  was  dominated  by  Ca  and 
SO4,  products  of  weathering.  There  was  little  change  in  precipitation  chemistry  be- 
yond 13  km  from  the  coast.  Nitrate-N  concentrations  were  higher  than  NH^-N 
concentrations  in  precipitation,  and  concentrations  were  low.  Sulfate-S  concentrations 
were  also  low.  The  chemistry  of  throughfall,  stemflow,  soil  solution,  and  the  stream  in 
the  V/est  Twin  Creek  watershed  was  markedly  different  from  that  in  precipitation.  The 
pH  averaged  5.3  in  precipitation,  5.0  in  throughfall,  4.3  in  stemflow,  5.7  in  the  forest 
floor,  6.2  in  soil  solution  (at  15  cm  depth),  and  7.3  in  the  stream.  The  concentration  of 
most  analytes  in  precipitation  and  in  the  West  Twin  Creek  stream  showed  few  trends 
from  1984  to  1993;  NO3-N  concentrations  in  precipitation  declined  slightly,  and  SO^- 
S  in  stream  water  increased.  From  fall  1993  through  1996,  however,  NO3-N 
concentrations  in  precipitation  and  in  the  stream  increased  and  pH  decreased  to  as  low 
as  4.6.  There  was  also  a  slight  increase  in  SO4-S.  This  seems  to  have  been  due  to 
across-Pacific  transport  of  pollutants  from  Asia.  The  old-growth  ecosystem  in  the  West 
Twin  Creek  watershed  is  on  the  verge  of  N  saturation,  which  has  many  implications 
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for  the  forest  health.  Long-term  monitoring  has  provided  the  opportunity  to  examine 
the  influence  of  increased  NO3-N  inputs  on  forest  streams  in  an  area  that  is  considered 
to  be  relatively  pristine. 

Key  words:  Vegetation,  hydrology,  precipitation  chemistry,  throughfall,  stemflow, 
stream  chemistry,  old-growth  forest. 
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Chapter  1.  Introduction 

Since  the  1970s,  there  has  been  widespread  concern  about  the  effects  of 
anthropogenic  atmospheric  inputs  of  sulfur  and  nitrogen  into  forest  ecosys- 
tems and  their  associated  streams  and  lakes  in  North  America  and  Europe.  To 
understand  these  effects,  several  major  research  programs  were  initiated  in  the 
1980s:  the  National  Acid  Precipitation  Assessment  Program  (NAPAP),  funded 
mostly  by  the  Environmental  Protection  Agency  (NAPAP  1993);  the  Integrated 
Forest  Study,  funded  by  the  Electric  Power  Research  Institute  (Johnson  and 
Lindberg  1992);  and  the  programs  of  Germany  and  other  European  countries. 
The  National  Acid  Precipitation  Assessment  Program  was  established  to  (1) 
develop  an  understanding  of  the  mechanisms  and  effects  of  atmospheric  depo- 
sition (acid  rain,  air  pollution)  on  forest  ecosystems,  (2)  examine  the  influence 
of  acidic  precipitation  and  other  air  pollutants  on  the  nation's  vegetation  and 
aquatic  ecosystems,  and  (3)  develop  an  objective  and  comprehensive  informa- 
tion base  for  decision  makers  (Schroeder  and  Keister  1989). 

Acidic  deposition  has  been  recognized  as  a  potentially  serious  threat  to  bio- 
logical systems  in  the  northeastern  United  States  and  in  parts  of  Europe 
(Overrein  1972;  Braekke  1976;  Wright  and  Gjessing  1976;  Maimer  1976; 
Schofield  1976;  Cronan  et  al.  1978).  The  southwestern  United  States  have  also 
experienced  acid  precipitation  events  (Lawson  and  Wendt  1982;  Olson  and 
Lefohn  1989;  Olson  et  al.  1992).  Acid  precipitation  can  have  several  detri- 
mental effects  on  forest  ecosystems.  They  include  ( 1 )  acidification  of  soils  and 
the  consequent  release  of  aluminum  ions,  which  are  toxic  to  tree  roots,  (2) 
leaching  of  important  nutrients  from  tree  canopies  and  soils  (particularly  cat- 
ions), (3)  reduction  of  tree  growth  and  increased  mortality,  (4)  the  adverse 
effects  of  heavy  metals  on  soils  and  ecosystem  processes,  and  (5)  acidification 
of  streams  and  effects  on  stream  biology.  In  addition,  excess  N  deposition  in 
some  areas  is  causing  problems  (Agren  and  Bosatta  1988;  Aber  et  al.  1989) 
including  NO3  losses  from  ecosystems  that  suffer  from  N  saturation. 

To  understand  the  effects  of  anthropogenic  atmospheric  deposition  on  natu- 
ral forest  ecosystems  in  the  United  States,  several  study  sites  were  established 
in  national  parks  in  the  early  1980s  as  part  of  NAPAP,  and  the  small  water- 
sheds approach  was  used.  The  parks  were  Sequoia-Kings  Canyon  in  California, 
Isle  Royale  in  Michigan,  Rocky  Mountain  in  Colorado,  Shenandoah  in  Vir- 
ginia, and  Olympic  in  Washington.  Sequoia-Kings  Canyon,  Rocky  Mountain, 
Isle  Royale,  and  Shenandoah  national  parks  were  in  areas  affected  by  anthro- 
pogenic inputs.  However,  anthropogenic  deposition  in  Olympic  National  Park 
was  small  and  provided  an  opportunity  to  examine  processes  in  a  relatively 
pristine  environment  or  sort  of  control  site. 

The  small-watershed  approach  has  been  used  successfully  at  Coweeta  in 
North  Carolina  (Swank  and  Crossley  1988),  Walker  Branch  in  Tennessee 
(Johnson  and  Van  Hook  1989),  and  Hubbard  Brook  in  New  Hampshire  (Lik- 
ens and  Bormann  1995)  to  study  biogeochemical  cycling  under  the  premise 
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that  chemical  fluxes  and  cycling  are  intimately  linked  to  hydrological  cycles. 
Understanding  the  biogeochemical  cycles  of  nitrogen,  sulfur,  carbon,  and  other 
elements  contributes  to  a  basic  understanding  of  ecological  processes  and  re- 
sponses of  ecosystems  to  stress.  Long-term  data  sets  can  point  out  changes  in 
chemical  cycling  patterns  and  species  composition  and  provide  baselines  against 
which  to  assess  future  conditions,  including  global  change.  We  were  particu- 
larly interested  in  determining  whether  nitrogen  and  sulfur  air  pollutants  from 
Asia  were  reaching  Olympic  National  Park,  and  if  so,  what  effects  they  were 
having  on  stream  chemistry. 

Although  there  is  considerable  knowledge  of  coastal  Pacific  Northwest  for- 
est ecosystems  (Edmonds  1982;  Franklin  and  Dyrness  1988;  Henderson  et  al. 
1989),  data  on  the  function  of  small  watersheds  are  limited,  especially  in  the 
temperate  rainforests  on  the  western  Olympic  Peninsula. 

Using  the  ecosystems  approach,  we  studied  terrestrial  and  aquatic  ecosys- 
tems in  two  small  watersheds  in  Olympic  National  Park:  low  elevation  West 
Twin  Creek  and  high  elevation  Hoh  Lake.  Both  watersheds  are  in  pristine, 
temperate,  old-growth  rainforests  in  the  Hoh  River  valley. 

From  1984  through  1993,  we  examined  vegetation  patterns,  biomass  distri- 
butions, hydrology,  solution  chemistry,  and  long-term  trends  in  precipitation 
and  stream  chemistry.  In  the  West  Twin  Creek  watershed,  we  established 
baseline  data  for  precipitation  chemistry,  canopy  interception  processes 
(throughfall  and  stemflow),  soil  solution,  hydrology,  and  stream  chemistry. 
Some  data  on  precipitation  and  stream  chemistry  and  hydrology  have  been 
assessed  through  1997. 

Chapter  2  of  this  monograph  presents  a  description  of  the  study  sites.  Sub- 
sequent chapters  present  methods,  results,  and  discussions  of  our  studies  of 
the  vegetation  communities  and  distribution  of  biomass  in  the  West  Twin  Creek 
and  Hoh  Lake  watersheds,  hydrologic  cycle  in  the  West  Twin  Creek  water- 
shed, solution  chemistry  in  the  West  Twin  Creek  watershed,  and  long-term 
trends  in  precipitation  and  stream  chemistry  in  the  West  Twin  Creek  water- 
shed. Program  publications  are  listed  in  Appendix  E. 
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Chapter  2:  Study  Sites 

Olympic  National  Park  is  located  along  a  narrow  coastal  strip  in  the  central 
portion  of  the  Olympic  Peninsula  in  the  northwestern  continental  United  States 
(Fig.  1).  The  west  side  of  the  Olympic  Peninsula  is  an  excellent  location  for 
examining  forest  ecosystems  that  are  unaffected  by  air  pollution  because  the 
inputs  of  atmospheric  particles  in  rain  and  dryfall  come  directly  off  the  Pacific 
Ocean.  It  has  been  hypothesized  that  coastal  forest  ecosystems,  such  as  those 
in  the  Hoh  River  valley,  should  be  as  clean  as  or  cleaner  than  any  ecosystem  in 
North  America. 

We  conducted  our  research  in  two  small,  pristine  watersheds  (West  Twin 
Creek— 124°  00'  36"  W,  47°  50'  00"  N  and  Hoh  Lake— 123°  AT  12"  W,  47°  53'- 
54"  N)  in  the  Hoh  River  valley  in  the  western  Olympic  National  Park  (Table 
1).  The  primary  site  for  the  study  of  hydrologic  and  forest  ecosystem  pro- 
cesses was  the  58-ha  West  Twin  Creek  watershed,  located  32  km  from  the 


20  Km 


West  Twin  Creek 
Watershed 


Figure  1.  Washington  state,  western  Olympic  Peninsula,  Hoh  River,  Hoh  Ranger  Sta- 
tion, National  Atmospheric  Deposition  Program  (NADP)  site,  and  the  precipitation 
collection  sites  and  their  proximity  to  the  Pacific  Ocean. 
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Pacific  Ocean  (Fig.  1).  Baseline  data  were  collected  on  precipitation, 
throughfall,  stemflow,  and  soil-solution  and  stream  chemistry.  Also,  atmospheric 
input  data  (precipitation  chemistry  and  rain  depth)  were  collected  for  the  Na- 
tional Atmospheric  Deposition  Program  at  the  Hoh  Ranger  Station, 
approximately  10  km  inland  from  the  West  Twin  Creek  watershed  (Fig.  1).  The 
125-ha  Hoh  Lake  watershed  is  located  50  km  from  the  Pacific  Ocean  (Fig.  1). 

West  Twin  Creek 

Topography 

West  Twin  Creek  is  a  steep,  dissected,  first-order  valley-wall  watershed. 
Elevations  of  the  watershed  range  from  180  m  at  the  weir  to  greater  than  850 
m  on  the  ridgetops  (Fig.  2).  The  watershed  has  steep  valley  walls,  and  the  first- 


Figure  2.  Topography  and  location  of  the  six  permanent  sampling  plots  in  the  West 
Twin  Creek  watershed,  Olympic  National  Park,  Washington. 
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order  stream  runs  southeast.  Slopes  range  from  33%-80%.  It  generally  faces 
southeast,  but  many  aspects  are  represented. 

Climate 

The  climate  of  the  old-growth  forest  of  the  West  Twin  Creek  watershed 
ranged  from  relatively  cool  and  moist  on  the  lower  slopes  to  warm  and  dry  on 
the  ridgetops.  At  the  Hoh  Ranger  Station  (Fig.  1),  located  10  km  east  of  the 
watershed  at  an  elevation  of  183  m,  mean  January  air  temperatures  were  4.5°C 
and  mean  July  air  temperatures  were  16°C  ( Table  1).  Average  daily  minimum 
soil  temperatures  in  the  lower  watershed  did  not  fall  below  0°C,  and  average 
daily  minimum  air  temperatures  only  occasionally  fell  below  0°C  (Fig.  3). 
Between  1984  and  1993,  average  annual  rainfall  at  the  Hoh  Ranger  Station 
was  about  350  cm  (National  Acid  Deposition  Program-National  Trends  Net- 
work 1994).  Most  of  the  rain  fell  between  October  and  May,  and  there  was  a 
pronounced  summer  dry  season.  Snow  rarely  fell  at  the  base  of  the  watershed 
and  generally  lasted  for  only  a  few  days.  A  winter  snow  pack  may  develop 
from  December  to  March  above  700  m  elevation. 
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Figure  3.  Average  daily  maximum  and  minimum  air  and  soil  temperatures  in  the  lower 
West  Twin  Creek  watershed  throughout  the  year,  1984- 1993. 
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Soils 

Soils  are  formed  from  marine  deposits  with  minor  conglomerates  and  brec- 
cias and  contain  some  green  and  black  sandstone  (Tabor  and  Cady  1978).  Soils 
vary  with  the  complex  topography.  Dominant  soils  of  the  watershed  are  loamy 
to  sandy-clay  textured  and  were  classified  as  coarse-loamy,  mixed,  mesic  Typic 
Dystrochrepts.  However,  a  moderate  degree  of  volcanic  ash  or  other  (e.g., 
spodic)  amorphous  material  present  suggested  that  soils  may  alternatively  be 
classified  as  medial  mesic  Andic  Dystrochrepts.  Because  of  windthrow  that  is 
common  in  the  area,  the  mineral  soil  was  mixed  (Appendix  A). 

Vegetation 

Western  hemlock  {Tsuga  heterophylla)  dominated  the  overstory  vegetation 
in  the  lower  watershed.  In  addition,  Douglas-fir  {Pseudotsuga  menziesii),  west- 
em  redcedar  (Thuja  plicata),  and,  at  the  lowest  elevations  in  the  watershed, 
Sitka  spruce  {Picea  sitchensis)  were  present  as  massive  old-growth  trees  (Fig.  4) 
and  sporadic  reproduction.  The  understory  vegetation  was  dominated  by 
swordfem  {Polystichum  mimitum)  and  Oregon  oxalis  {Oxalis  oregana). 

Pacific  silver  fir  [Abies  amabilis)  and  western  hemlock  were  the  dominant 
overstory  species  in  the  upper  watershed,  but  western  redcedar  was  also  present. 
The  main  understory  species  were  oval-leaf  huckleberry  {Vaccinium 
ovalifolium),  red  huckleberry  {Vaccinium parvifolium)  and  Oregon  oxalis.  Salal 
{Gaultheria  shallon)  occurred  on  south-facing  slopes  at  mid  to  upper  eleva- 
tions (Appendix  B). 


Figure  4.  Large,  old-growth  Douglas-fir  tree  in  the  West  Twin  Creek  watershed. 
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Hardwoods  were  not  common  in  the  West  Twin  Creek  watershed.  However, 
riparian  areas  contained  bigleaf  maple  (Acer  macrophyllum),  and  disturbed 
areas,  such  as  slope  failures  and  creep  slopes,  contained  red  alder  (Alnus  rubra). 

The  tree  population  was  composed  of  82%  western  hemlock,  10%  Pacific 
silver  fir,  6%  western  redcedar,  1%  Douglas-fir,  and  1%  Sitka  spruce.  How- 
ever, only  26%  of  the  tree  basal  area  was  composed  of  western  hemlock.  The 
remaining  basal  area  was  composed  of  26%  Pacific  silver  fir,  25%  western 
redcedar,  15%  Douglas-fir,  and  6%  Sitka  spruce. 

The  forest  was  uneven-aged.  The  age  of  Douglas-firs  ranged  from  237  to 
635  years,  western  redcedars  from  123  to  600  years,  western  hemlocks  from 
71  to  262  years,  and  Pacific  silver  firs  from  98  to  248  years.  Maximum  tree 
height  was  greater  than  90  m,  but  the  canopy  was  broken  up  and  uneven.  Di- 
ameter at  breast  height  (DBH)  of  Douglas-firs  ranged  from  10  to  310  cm, 
western  redcedars  from  27  to  302  cm,  western  hemlocks  from  5  to  260  cm, 
and  Pacific  silver  firs  from  9  to  190  cm. 

Hoh  Lake 

Topography 

The  Hoh  Lake  watershed  is  a  high-elevation  subalpine  site  (Fig.  5).  Slopes 
range  from  55%  to  78%,  and  the  watershed  generally  faces  south  to  southeast. 
Elevations  range  from  1250  to  1525  m;  the  lake  is  at  1390  m  (Fig.  6). 


Figure  5.  View  of  Hoh  Lake  and  the  watershed  from  higher  elevation. 
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Climate 

Mean  January  air  temperatures  were  -2.5°C  and  mean  July  air  temperatures 
were  15°C  (Table  1).  Annual  precipitation  ranges  from  380  to  500  cm;  most 
fell  between  October  and  May,  and  there  was  a  pronounced  summer  dry  sea- 
son. Most  of  the  precipitation  fell  as  snow,  and  typical  snowpack  depths  ranged 
from  760  to  1300  cm. 

Soils 

Soils  are  formed  from  marine  deposits  of  grayish-green  and  black  sand- 
stone (Tabor  and  Cady  1978)  and  were  classified  as  Typic  Cryorthods  or  Typic 
Cryumbrepts.  Soils  vary  with  the  topography.  Dominant  soils  of  the  watershed 
have  sandy-loam  texture  (Appendix  A). 

Vegetation 

The  overstory  vegetation  of  the  Hoh  Lake  watershed  contains  three  tree 
species:  Pacific  silver  fir,  mountain  hemlock  (Tsuga  mertensiana),  and  Alaska 
yellow  cedar  {Chamaecyparis  nootkatensis).  The  understory  vegetation  was 
dominated  by  Alaska  huckleberry  {Vaccinium  alaskaense),  big  huckleberry 


Figure  6.  Topography  of  the  Hoh  Lake  watershed  and  the  locations  of  the  permanent 
sampling  plots. 
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(Vaccinium  membranaceum),  avalanche  fawnlily  {Erythronium  montanum),  and 
common  beargrass  (Xemphyllum  tenax\  Appendix  B). 

The  forest  was  uneven-aged.  The  age  of  Pacific  silver  firs  ranged  from  155 
to  452  years,  mountain  hemlocks  from  198  to  482  years,  and  Alaska  yellow 
cedars  from  200  to  438  years.  Maximum  tree  height  was  5 1  m,  and  the  canopy 
was  uneven.  Maximum  DBH  of  Pacific  silver  fir  was  123.9  cm,  mountain  hem- 
lock 122.3  cm,  and  Alaska  yellow  cedar  65.7  cm. 
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Chapter  3:  Vegetation  Communities 

in  the  West  Twin  Creek  and 

Hoh  Lake  Watersheds 

Vegetation  in  old-growth  forests  in  Oregon  and  Washington  and  the  rela- 
tions between  the  environment  and  vegetation  have  been  described  in  some 
detail  (Fonda  and  Bliss  1969;  Fonda  1974;  Zobel  et  al.  1976;  del  Moral  and 
Long  1977;  Grier  and  Logan  1977;  Hemstrom  et  al.  1982;  Agee  and  Kertis 
1987;  Franklin  and  Dyrness  1988;  Henderson  et  al.  1989;  Kirk  and  Franklin 
1992).  Grier  and  Logan  (1977)  presented  a  detailed  vegetation  description  of 
an  old-growth  watershed  in  the  Oregon  Cascades,  but  detailed  descriptions  of 
old-growth  forest  vegetation  on  a  watershed  scale  in  the  Olympic  Mountains 
are  not  available.  The  present  vegetation  is  only  a  remnant  of  the  historic  veg- 
etation, and  most  of  it  is  in  Olympic  National  Park  (Kirk  and  Franklin  1992).  A 
better  understanding  of  the  spatial  variability  in  vegetation  patterns  in  water- 
sheds will  be  helpful  in  developing  ecologically  sound  management  practices 
at  the  landscape  level  (Franklin  1993). 

Our  objectives  were  the  description  and  classification  of  the  vegetation  in 
low  (West  Twin  Creek)  and  high  (Hoh  Lake)  elevation  old-growth  forested 
watersheds  in  the  Hoh  River  Valley,  Washington,  and  the  examination  of  rela- 
tions between  the  environment  and  plant  species  and  communities. 

Methods 

One  hundred  and  sixty-seven  circular  reconnaissance  plots  (500  m^)  were 
established  in  the  West  Twin  Creek  watershed  in  the  summer  of  1984.  From 
approximately  mid-slope  baselines  at  West  Twin  Creek,  transects  were  posi- 
tioned upslope  and  downslope.  Plots  were  spaced  at  100-m  intervals  up  to  the 
ridge  top  and  down  to  creek  level.  This  spacing  ensured  sampling  the  full  range 
of  plant  communities.  Nested  circular  plots  were  used  to  sample  the  vegeta- 
tion in  each  plot.  Seedling  presence  by  species  was  determined  in  four  20-m^ 
subplots.  In  one  50-m^  subplot,  percent  cover  of  herbs,  shrubs,  moss  lichen,  bare 
soil,  litter  and  rocks  (>8  cm)  were  recorded.  Trees  were  measured  and  coded  in 
250  m^  subplots,  and  large  trees  (>50  cm  dbh)  and  all  snags  were  measured  in  the 
whole  reconnaissance  plot.  Using  the  reconnaissance  data,  five  plant  community 
types  and  a  variant  of  one  type  were  identified  in  the  watershed.  A  permanent 
circular  sample  plot  (1000  m^)  was  randomly  located  in  each  of  these  six  types 
(Fig.  7).  Nested  circular  subplots  were  used  to  sample  the  understory  vegetation  in 
each  plot:  four  20-m^  subplots  for  tree  seedlings  and  one  50-m^  subplot  for  percent 
cover  of  herbs,  shrubs,  mosses,  lichens,  bare  soil,  litter  and  rocks  (>8  cm).  Percent 
slope,  aspect,  elevation,  topographic  shape,  and  slope  position  were  determined  in 
each  plot.  At  Hoh  Lake,  five  circular  permanent  sample  plots  (1000  m^)  were 
established  after  delineation  of  the  plant  communities  of  the  closed-canopy  forest 
in  the  basin  in  an  on-the-ground  reconnaissance  (Fig.  6). 
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Tsuga  heterophylla- 
Abies  amabilis/ 
Vaccinium  alaskaense 


Tsuga  heterophylla- 
Abies  amabilis/ 
Depauperate 


Tsuga  heterophylla/ 
Gaultheria  shallon 


Tsuga  heterophylla/ 
Oxalis  oregana 


Tsuga  heterophylla/ 
Polystichum  munitum- 
Oxalis  oregana 


Figure  7.  Major  plant  communities  in  the  West  Twin  Creek  watershed  and  the  loca- 
tions of  the  permanent  sampling  plots. 


All  measured  trees  were  assessed  for  vigor  and  condition.  These  conditions 
included  canopy  class  (dominant,  codominant,  intermediate,  or  suppressed), 
canopy  condition,  bole  condition,  any  signs  of  mechanical  or  animal  distur- 
bance, cavities  that  would  serve  as  courts  for  disease  or  insects,  and  signs  of 
disease  or  insect  damage  or  symptoms.  In  addition,  crown  damage  and  die- 
back  were  used  as  measures  of  the  condition  of  individual  trees. 

Detrended  correspondence  analysis  (DCA),  an  ordination  technique,  was 
used  to  analyze  the  data.  Species  and  sample  ordinations  were  calculated  with 
DECORANA  (Hill  1979a).  The  samples  scores  were  weighted  means  of  the 
species  scores  within  each  sample.  Data  were  classified  with  the  two-way  in- 
dicator analysis  program  TWINSPAN  (Hill  1979b). 


Results  and  Discussion 

Plant  Communities  in  the  West  Twin  Creek  Watershed 

Plots  1,  2,  and  5  were  in  the  lower  watershed  and  plots  3,  4,  and  6  were  in 
the  upper  watershed  (Fig.  7,  Table  2).  Average  plot  DBH  ranged  from  27.4  to 
49.6  cm  and  total  basal  area  ranged  from  56.4  to  135.2  m^/ha.  Understo'ry 
community  composition  and  structure  were  also  highly  variable,  ranging  from 
5  to  18  understory  species/0.1 -ha  plot  (Appendix  C). 
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Figure  8.  Permanent  sampling  plot  1  in  the  lower  West  Twin  Creek  watershed. 


Lower  Watershed 

The  most  widespread  community  type  in  the  entire  watershed,  Tsuga 
heterophylla-Polystichum  munitum-Oxalis  oregana,  was  located  on  the  lower 
valley  walls  adjacent  to  the  creek  (Fig.  7).  The  two  variants  of  this  community 
type  were  represented  in  plot  1  (Fig.  8)  and  plot  2  (Fig.  9).  Plot  1  was  on  a 
south  slope  (205°  aspect)  at  an  elevation  of  342  m  on  a  67%  slope  and  con- 
tained warm  and  dry  indicator  species  salal  and  red  huckleberry  (Brockway  et 
al.  1983;  Lesher  and  Henderson  1989).  Western  hemlock  was  the  most  abun- 
dant tree  species  in  the  upper  canopy  (72%  of  the  total;  Table  2)  and  dominated 
the  regeneration  and  intermediate  size  classes.  Douglas  fir  represented  20%  of 
the  stocking,  including  some  scattered  large  old  trees,  whereas  western  redcedar 
represented  8%  of  the  stocking,  including  some  large  specimens  (Table  2). 
Sword  fern  and  Oregon  oxalis  were  the  most  important  understory  species 
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Figure  9.  Permanent  sampling  plot  2  in  the  lower  West  Twin  Creek  watershed. 


(Appendix  C).  Major  associated  species  in  this  community  type  included  red 
huckleberry  and  Oregon  beaked  moss  {Eurhynchium  oreganum\  Appendix  C). 
Western  hemlock  was  the  dominant  regenerating  species  at  18  000  seedlings/ 
ha  (Table  3).  Some  Sitka  spruce  and  western  redcedar  seedlings  were  also 
present,  but  no  regeneration  of  Douglas  fir  or  Pacific  silver  fir  was  found. 

Plot  2  was  on  a  north  slope  (30  degrees  aspect)  at  an  elevation  458  m  on 
75%  slope  (Table  2)  and  had  a  much  richer  herb  layer  with  75%  cover  by  six 
moss  species:  Eurhynchium  oreganum,  Rhytidiadelphus  loreus,  Hylocomium 
splendens,  Polytrichum  juniperinum,  Dicranum  scoparium  and  Mnium  (Ap- 
pendix C).  Western  hemlock  was  the  most  abundant  tree  species  (Tables  2  and  3). 
Tree  stocking  was  less  than  in  plot  1  (375  versus  421  trees/ha),  and  the  basal 
area  (56.4  m^/ha)  was  less  than  half  of  that  in  plot  1(121.1  m'/ha;  Table  2). 
The  reason  for  this  difference  may  have  been  the  greater  disturbance  in  plot  2 
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Table  3.  Stem  density  (stems/ha)  of  tree  species  <5  cm  diameter  at  breast 
height  (DBH)  in  permanent  sampling  plots  in  the  West  Twin  Creek  watershed, 
1985. 


Species 

Plot  1 

Plot  2 

Plots 

Plot  3 

Plot  4 

Plot  6 

western  hemlock 

18000 

32300 

9200 

12400 

32400 

1500 

Douglas-fir 

0 

200 

0 

700 

400 

0 

Pacific  silver  fir 

0 

200 

0 

4400 

500 

7200 

western  redcedar 

600 

900 

300 

300 

0 

0 

Sitka  spruce 

1100 

500 

0 

0 

0 

0 

as  evidenced  by  the  high  percentage  of  bare  soil  (10%)  and  high  rock  content 
on  the  plot  (Appendix  C).  Douglas  fir  was  also  regenerating  on  this  plot  at  200 
seedlings/ha  (Table  3).  Five  tree  species  were  represented  in  the  <5-cm  diam- 
eter class,  although  western  hemlock  strongly  dominated.  As  in  plot  1,  the 
understory  vegetation  was  dominated  by  swordfem  and  Oregon  oxalis;  major 
associates  were  mosses  and  lichens. 

Plot  5  (Fig.  10)  was  located  at  458  m  on  a  74%  slope  on  a  southerly  aspect 
in  a  Tsuga  heterophylla  Gaultheria  shallon  community  type.  The  plot  con- 
tained 983  trees/ha  and  135.2  m'^/ha  of  basal  area  (Table  2).  Western  hemlock 
was  represented  in  all  size  classes  and  made  up  82%  of  the  trees.  Thus,  it 
dominated  the  large  size  classes  and  made  up  28%  of  the  basal  area;  a  few 
mid-size  Pacific  silver  firs  were  present.  Western  hemlock  also  dominated  the 
<5-cm-diameter  class  at  9,200  trees/ha  (Table  3).  Western  redcedar  at  300  trees/ 
ha  was  the  only  other  species.  Salal  covered  60%  of  the  plot  (Appendix  C). 
Vaccinium  species  and  vine  maple  {Acer  circinatum)  were  shrub  associates; 
only  trace  amounts  of  herbs  were  present  (Appendix  C). 

Upper  watershed 

As  elevation  in  the  watershed  increased,  vegetation  gradually  changed  from 
dominance  by  western  hemlock  to  dominance  by  Pacific  silver  fir.  Plot  3 
(Fig.  1 1)  was  located  in  a  transition  community  type  {Tsuga  heterophylla-Abies 
amabilis  Vaccinium  alaskaense).  This  plot  was  located  on  the  steep  southeast 
headwall  near  the  top  of  the  watershed,  just  below  the  Hoh  Bogachiel  Trail. 
The  elevation  of  the  plot  was  714  m,  and  the  slope  was  80%.  The  overstory 
was  dominated  by  large  Pacific  silver  fir  and  western  redcedar  (Table  2).  West- 
em  hemlock  was  present  in  the  suppressed  and  intermediate  size  classes  and 
had  the  highest  stem  density  (295  trees/ha).  Western  hemlock  and  Pacific  sil- 
ver fir  were  regenerating  in  the  plot;  hemlock  regeneration  dominated  over 
silver  fir  3  to  1  (Table  3).  Some  Douglas  fir  and  western  redcedar  were  also 
regenerating.  The  understory  was  dominated  by  Alaska  huckleberry  {Vaccinium 
alaskaense),  red  huckleberry,  and  salmonberry  {Rubus  spectabilis).  Except 
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Figure  10.  Permanent  sampling  plot  5  in  the  lower  West  Twin  Creek  watershed. 


deerfern  (Blechnum  spicant),  herb  species  were  present  in  only  trace  amounts 
(Appendix  C). 

Plot  4  was  in  the  Tsuga  heterophylla-Oxalis  oregana  community  and  was 
located  on  a  gentle  east-facing  slope  at  586  m  (Fig.  12).  This  community  at 
239  trees/ha  and  a  basal  area  of  72.2  m^/ha  had  a  lower  stem  density  than  the 
other  types  (Table  2).  Western  hemlock  was  the  dominant  tree  in  the  overstory. 
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Figure  11.  Permanent  sampling  plot  3  in  the  upper  West  Twin  Creek  watershed. 


regeneration  (Table  3),  and  intermediate  size  classes.  Pacific  silver  fir  was  the 
only  other  mature  tree  species  but  was  not  abundant  in  any  size  class.  Douglas 
fir  regeneration  was  found  in  gaps  created  by  canopy  openings.  Shrubs  were 
conspicuously  absent;  the  understory  was  covered  by  a  rich  herb  layer.  Oregon 
oxalis  and  deerfoot  vanillaleaf  (Achlys  triphylla)  covered  80%  of  the  plot,  and 
white  trillium  {Trillium  ovatum).  Smith's  fairybells  {Disporum  smithii),  and 
western  springbeauty  (Montia  sibirica)  were  associates  (Appendix  C). 

Plot  6  (Fig.  13)  was  classified  as  Tsuga  hetewphylla-Abies  amabilis  Depau- 
perate type.  This  community  type  was  located  on  a  warm,  dry,  upper-slope 
position.  The  slope  was  gentle  (33%)  and  faced  south.  This  type  had  little  or 
no  understory  development  of  regenerating  conifers  or  understory  vegetation. 
Traces  of  Alaska  huckleberry  and  northwestern  twayblade  {Listera  caurina) 
were  present.  The  lack  of  understory  vegetation  was  probably  due  to  the  dense 
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Figure  12.  Permanent  sampling  plot  4  in  the  upper  West  Twin  Creek  watershed. 


Figure  13.  Permanent  sampling  plot  6  in  the  upper  West  Twin  Creek  watershed. 


stocking  of  mature  trees  (948/ha).  The  south-facing  ridgetop  position  of  this 
community  type  was  particularly  susceptible  to  wind  disturbance  that  prevails 
from  the  south  and  west.  The  accumulation  of  litter  in  this  stand  was  high; 
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2M  Abies  amabilis/Rhododendron  albiflorum 


D 


Abies  amabilis-Tsuga  mertensiana/ 
Vaccinium  membranaceum 

Abies  amabilis-Tsuga  mertensiana/ 
Erythronium  montanum 
Subalpine  meadow  (Vaccinium  deliciosum^ 
Phyllodoce  sp.) 
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n 


Burned  area 


Talus 


Avalanche  tracks 


CB  Flat 


Figure  14.  Major  plant  communities  and  disturbed  areas  in  the  Hoh  Lake  watershed 
and  the  locations  of  the  permanent  sampling  plots. 


large  logs  of  varying  decay  classes  and  recent  fresh  branch  and  twig  litter 
covered  most  of  the  plot. 

Plant  Communities  in  the  Hoh  Lake  Watershed 

The  environment  at  Hoh  Lake  is  predominantly  cold  and  wet.  Gale  force 
winds  buffet  the  ridgetops  above  Hoh  Lake,  and  forests  are  entirely  absent. 
Just  below  the  elevation  of  Hoh  Lake  (1390  m),  forests  were  continuous  ex- 
cept where  avalanche  tracks  occurred.  Three  plant  community  types  occurred 
in  the  closed-cano/?y /or^^r  (Abies  amabilis-Tsuga  mertensiana  Vaccinium 
membranaceum,  Abies  amabilis-Tsuga  mertensiana  Erythronium  montanum, 
and  Abies  amabilis  Rhododendron  albiflorum)  along  with  a  subalpine  com- 
munity (Figs  6  and  14).  Above  the  lake  and  at  CB  Flat,  a  subalpine  wet  meadow 
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Vaccinium  deliciosum-Phyllodoce  spp.  community  occurred.  The  meadow  at 
CB  Flat  covered  approximately  6  ha  of  gentle  ground  (<20%  slope).  Commu- 
nity composition  and  structure  were  slightly  simpler  at  Hoh  Lake  (8-13 
understory  species/0.1 -ha  plot;  Appendix  C)  than  at  West  Twin  Creek  (5-18 
understory  species/0.1 -ha  plot). 

Five  plots  were  established  in  the  closed-canopy  forest.  Plot  1  at  1295-m 
elevation  and  plot  2  at  1257-m  elevation  were  located  in  the  Abies  amabilis- 
Tsuga  mertensiana  Vaccinium  membranaceum  community  on  the  south-facing 
slopes  east  of  CB  Flat  (Fig.  14).  Plot  1  had  a  slope  of  55%;  plot  2  had  a  slope 
of  78%.  Both  plots  had  abundant  exposed  bedrock.  More  than  50%  of  the  plot 
areas  were  covered  by  large  litter  (>10  cm)  and  coarse  woody  debris  (Appen- 
dix C).  The  overstory  in  plot  1  was  dominated  by  Alaska  yellow  cedar,  but 
Pacific  silver  fir  and  Alaska  yellow  cedar  were  codominant  in  plot  2  (Table  4). 
Basal  area  was  121.3  m^/ha  in  plot  1  and  125.5  m^/hain  plot  2  (Table  4).  Tree 
density  ranged  from  799  stems/ha  in  plot  1  to  989  stems/ha  in  plot  2.  Regen- 
eration of  Alaska  yellow  cedar  and  mountain  hemlock  was  also  present,  but 
regeneration  was  clearly  dominated  by  Pacific  silver  fir  at  over  15  000  stems/ 
ha  (Table  5).  Big  huckleberry  was  the  dominant  understory  shrub,  and  there 
was  a  smaller  cover  of  Alaska  huckleberry.  Except  for  small  amounts  of  ava- 
lanche fawnlily  and  beargrass,  only  trace  amounts  of  herbs  were  present. 

Plot  3  (Fig.  15)  at  1290-m  elevation  and  plot  5  at  1300-m  elevation  faced 
east  and  were  in  Abies  amabilis-Tsuga  mertensiana  Erythronium  montanum 
associations.  Slopes  in  both  plots  were  similar  (65%-70%),  and  the  plots  were 


Figure  15.  Permanent  sampling  plot  3  in  the  Hoh  Lake  watershed. 
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Table  5.  Stem  density  (stems/ha)  of  tree  species  <5  cm  diameter  at  DBH  in 
permanent  sampling  plots  in  the  Hoh  Lake  watershed,  1985. 


Species 

Plot  1 

Plot  2 

Plot  3 

Plot  4 

Plot  5 

Pacific  silver  fir 
Alaska  yellow  cedar 
Mountain  hemlock 

15600 
0 
0 

25200 

4800 

0 

25900 

100 

0 

12700 

3700 

200 

22200 

0 

300 

open  with  few  trees,  averaging  only  358  stems/ha  in  plot  3  and  354  stems/ha  in 
plot  5.  The  overstory  was  codominated  by  Pacific  silver  fir  in  plot  3,  but  by 
mountain  hemlock  in  plot  5  (Table  4).  Four  huckleberry  species  dominated  the 
shrub  layer  with  6%  cover,  and  avalanche  fawnlily  averaged  greater  than  5% 
of  the  cover.  Other  forbs  were  only  in  trace  amounts  (Appendix  C). 

Plot  4  at  1310-m  elevation  faced  northeast  and  was  located  in  the  Abies 
amabilis  Rhododendron  albiflorum  association.  Tree  density  was  527  trees/ 
ha;  Pacific  silver  fir,  western  hemlock,  and  Alaska  yellow  cedar  had  similar 
proportions  in  the  overstory  (Table  4).  Basal  area  (58.8  m-^/ha)  was  the  lowest 
of  the  five  plots  at  Hoh  Lake.  However,  the  shrub  cover  was  high,  averaging 
more  than  55%;  35%  was  white  rhododendron  {Rhododendron  albiflorum). 
Avalanche  fawnlily  cover  was  as  high  as  20%,  and  this  was  the  only  plot  where 
pink  mountain  heather  {Phyllodoce  empetriformis)  was  present  (15%;  Appen- 
dix C).  The  arboreal  lichen,  common  witch's  hair  {Alectoria  sarmentosa),  was 
visibly  present  throughout  the  closed-canopy  forest  of  the  Hoh  Lake  Basin, 
particularly  in  plots  3  and  4. 

Plant  Community  Classification  and  Ordination 

Elevation  is  extremely  important  in  explaining  the  vegetation  distribution 
in  the  West  Twin  Creek  and  Hoh  Lake  watersheds  in  the  Hoh  River  valley,  and 
elevation  corresponds  to  DC  A  Axis  I  in  Figure  16.  Temperature  and  other  en- 
vironmental factors  are  integrated  by  elevation.  Fonda  and  Bliss  (1969)  also 
suggested  that  elevation  was  extremely  important  in  explaining  the  vegetation 
distribution  in  the  Olympic  mountains.  Furthermore,  Agee  and  Kertis  (1987) 
found  that  DC  A  Axis  1  corresponded  to  elevation  in  an  analysis  of  forest  types 
in  North  Cascades  National  Park,  and  del  Moral  and  Long  (1977)  found  that 
elevation  was  responsible  for  much  of  the  variation  among  community  types 
in  the  Cedar  River  watershed  in  the  central  Washington  Cascade  mountains. 

Axis  2  in  Figure  16  seems  to  correspond  to  slope.  In  the  Cascade  Moun- 
tains, del  Moral  and  Long  (1977)  and  Agee  and  Kertis  (1987)  found  the  second 
factor  to  be  a  moisture  gradient,  not  slope.  However,  on  the  west  slopes  of  the 
Olympics,  there  is  less  of  a  moisture  gradient,  so  that  the  influence  of  slope 
increases  in  importance. 

The  Hoh  River  valley  contains  a  wide  array  of  species  and  community  types 
(Fig.  17).  Considerable  variability  is  related  to  changes  in  elevation,  reflecting 
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ABAM  -  Abies  amabilis 

CHNO  -  Chamaecyparis  nootkatensis 

PSME  -  Pseudotsuga  menziesii 

THPL  -  Thuja  plicata 

TSHE  -  Tsuga  heterophylla 

TSME  -  Tsuga  mertensiana 


SHRUBS  AND  HERBS 


GASH  -  Gaultheria  shallon 
MEFE  -  Menziesia  ferruginea 
POMU  -  Polystichum  munitum 
OXOR  -  Oxalis  oregana 
RUSP  -  Rubus  spectabilis 
VAAL  -  Vaccinium  alaskaense 
VAME  -  Vaccinium  membranaceum 
VAPA  -  Vaccinium  parvifolium 
XETE  -  Xerophyllum  tenax 


Figure  16.  Species  ordination  of  the  Hoh  River  valley  and  specific  ordination  of  15 
diagnostic  species.  Axis  1  corresponds  to  elevation,  and  Axis  2  corresponds  to  slope. 


a  strong  temperature  gradient.  Henderson  et  al.  (1989)  also  found  that  in  the 
Olympic  Mountains  vegetation  type  strongly  correlated  with  elevation.  In  the 
two  small  watersheds,  variability  in  vegetatiori  was  considerable  and  mostly 
related  to  slope  differences. 
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Figure  17.  Ordinations  of  major  plant  communities  in  the  Hoh  River  valley  based  on 
108  sample  plots.  Axis  1  corresponds  to  elevation,  and  axis  2  corresponds  to  slope. 
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Chapter  4:  Biomass  Distributions 

in  the  West  Twin  Creek  and 

Hoh  Lake  Watersheds 

In  recent  years,  interest  has  been  great  in  organic  matter  accumulations  in 
the  world's  forests  related  to  global  warming.  Organic  matter  accumulations 
in  forests  could  directly  affect  global  climate  change  by  sequestering  C  through 
plant  growth  and  releasing  C  through  decomposition  (Dixon  et  al.  1994).  Com- 
pilation of  organic  matter  budgets  for  ecosystems  is  important  if  we  are  to 
determine  the  effects  of  global  change  and  forest  management  on  C  pools  and 
fluxes  (Harmon  et  al.  1990;  Dixon  et  al.  1994;  Turner  et  al.  1995). 

Interest  in  organic  matter  accumulations  in  Pacific  Northwest  United  States 
forests  extends  beyond  C  pools  and  fluxes.  Dead  and  downed  coarse  woody 
debris  (CWD)  as  a  structural  and  functional  element  of  forests  seems  valuable 
in  promoting  ecosystem  stability,  habitat  diversity,  and  long-term  productivity 
(Franklin  1992;  Swanson  and  Franklin  1992;  Caza  1993).  Forest  managers  are 
being  urged  to  leave  CWD  on  harvested  lands  to  provide  important  habitat  for 
vertebrates,  invertebrates,  and  plants,  to  provide  soil  organic  matter  and  nutri- 
ents in  the  long-term,  and  to  reduce  erosion  (Harmon  et  al.  1986;  Franklin  et 
al.  1989,  Franklin  1992).  However,  how  much  CWD  must  be  left  in  various 
plant  communities  to  fulfill  certain  ecological  functions  remains  unclear.  Natural 
levels  of  CWD  can  vary  greatly  among  different  ecosystem  types,  depending 
on  climatic  influences  on  plant  growth  and  decomposition. 

Old-growth  forests  of  the  Pacific  Northwest  are  among  the  largest  forests  in 
the  world  and  can  accumulate  impressive  amounts  of  living  and  dead  organic 
matter  (Grier  and  Logan  1977;  Franklin  and  Waring  1979).  Forest  communi- 
ties near  the  coast  may  be  expected  to  have  very  high  biomass  accumulations 
because  high  moisture  levels  could  increase  producfivity  and  lengthen  fire  fre- 
quencies. Agee  and  Huff  (1987)  found  accumulations  of  CWD  that  exceeded 
500  Mg/ha  in  old-growth  forests  of  the  western  Olympic  Mountains;  Keenan 
et  al.  (1993)  found  a  similar  maximum  accumuladon  of  400  Mg/ha  in  old- 
growth  sites  on  northern  Vancouver  Island,  Canada.  Although  vegetafion  and 
CWD  biomass  in  coastal  old-growth  forests  have  been  estimated  by  several 
authors  (Graham  and  Cromack  1982;  Agee  and  Huff  1987;  Harcombe  et  al. 
1990),  few  total  organic  matter  budgets  have  been  compiled. 

Similarly,  ecosystem  biomass  budgets  have  been  compiled  of  only  few  high 
elevation  forests  in  the  Pacific  Northwest.  An  organic  matter  budget  was  com- 
piled of  a  high-elevation,  old-growth  Pacific  silver  fir  forest  in  the  Washington 
Cascades  (Grier  et  al.  198 1)  where  total  ecosystem  organic  matter  was  1247  Mg/ 
ha.  A  recent  study  by  Prichard  (1996),  which  provided  the  first  estimate  of  ecosys- 
tem biomass  in  the  Olympic  Mountains,  reported  that  total  ecosystem  biomass 
ranged  from  569  to  957  Mg/ha  in  Pacific  silver  fir-mountain  hemlock  forests. 
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Our  objective  in  this  research  was  the  compilation  of  organic  matter  bud- 
gets of  living  and  detrital  material  in  two  old-growth  watersheds:  the 
low-elevation  West  Twin  Creek  watershed  and  the  high-elevation  Hoh  Lake 
watershed.  In  the  West  Twin  Creek  watershed,  biomass  of  overstory  and  under- 
story,  CWD,  forest  floor,  soil  organic  matter,  and  fme  roots  were  determined;  in 
the  Hoh  Lake  watershed  only  overstory  vegetation  and  CWD  were  determined. 

Methods 

Field  Measurements  and  Calculations 

In  each  permanent  plot  in  the  West  Twin  Creek  and  Hoh  Lake  watersheds, 
overstory  and  understory  vegetation  and  coarse  woody  debris  were  measured 
to  estimate  biomass.  All  live  stems,  logs,  and  snags  were  mapped.  In  the  West 
Twin  Creek  watershed,  the  depth  of  the  forest  floor  was  measured,  and  forest 
floor  and  mineral  soil  samples  were  taken  to  determine  the  C  content.  Fine 
root  biomass  was  estimated  in  only  one  plot  of  the  West  Twin  Creek  water- 
shed. All  values  were  adjusted  for  slope. 

Live  Vegetation  Biomass 

In  each  permanent  plot  in  both  watersheds,  diameter  at  breast  height  (DBH) 
of  trees  greater  than  5  cm  was  measured  and  species  were  recorded  in  1985. 
Biomass  of  the  overstory  component  of  each  species  in  the  plots  was  calcu- 
lated with  equations  developed  by  Gholz  et  al.  (1979).  No  calculations  were 
performed  if  the  DBH  was  less  than  7  cm  because  biomass  was  negative  below 
that  point.  Tree  biomass  was  separated  into  components  of  foliage,  branches, 
roots,  and  bole-and-bark.  Biomass  of  coarse  roots  was  calculated  from  two 
sources:  Pacific  silver  fir  root  biomass  was  calculated  from  equations  in  Keyes 
(1982),  and  all  other  species  were  calculated  with  the  equations  for  Douglas  fir 
in  Gholz  et  al.  (1979).  Allometric  equations  for  root  biomass  have  not  been 
developed  for  the  other  tree  species  in  the  West  Twin  Creek  and  Hoh  Lake 
watersheds. 

In  the  West  Twin  Creek  plots,  understory  vegetation  was  sampled  for  biom- 
ass with  ten  1-m^  nested  circular  plots.  In  these  subplots,  all  understory  plants 
were  clipped  to  ground  level  at  the  height  of  the  growing  season  (late  July 
1985).  Biomass  samples  were  taken  to  the  laboratory  and  stems  and  foliage 
were  separated  by  live  and  dead  categories. 

In  May  1992,  three  soil  cores  measuring  5.08  cm  x  20  cm  were  collected  in 
Plot  1  of  the  West  Twin  Creek  watershed.  Fine  roots  of  less  than  2  mm  in 
diameter  were  sorted  from  the  cores  and  were  dried  at  75°C  for  48  hours.  Fine 
root  biomass  was  estimated  in  Mg/ha. 

Coarse  Woody  Debris 

Decay  class,  physical  dimensions,  and  species  of  snags  and  fallen  logs  were 
recorded  in  each  permanent  plot.  The  DBH  and  height  of  snags  were  recorded. 
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The  length  and  diameter  of  both  ends  of  logs  were  recorded.  When  a  log  ex- 
tended beyond  the  permanent  plot,  only  the  portion  in  the  plot  was  measured. 
These  physical  dimensions  were  converted  to  biomass  in  the  following  manner. 

Calculating  biomass  of  snags  and  logs  involved  correcting  for  wood  density 
of  each  decay  class  and  then  calculating  volume  and  biomass.  Density  values 
were  used  from  Means  et  al.  (1992).  These  values  were  derived  for  Douglas 
fir,  but  were  used  here  for  all  measured  species.  Smalian's  formula  was  used  to 
calculate  the  volume  of  logs  (Hartman  et  al.  1981): 

Volume  =  0.00007854  *  (Dl  +  D2)/2  *  L/10 

where  D 1  and  D2  equal  the  end  diameters  and  L  equals  the  length  of  the 
log.  These  volumes  were  converted  to  biomass  (Mg/ha)  with  the  following 
formula: 

Biomass  =  volume  (m^)  *  density  (g/cm^)  *  (kg/lO^g) 
*  (10^  cm3/m3)  *  (104  m^/ha)  *  (Mg/lO^  kg) 

The  same  density  values  were  used  for  both  logs  and  snags.  Snag  volume 
was  calculated  as  the  frustum  of  a  cone  (Husch  et  al.  1972)  with  an  assumed 
taper  of  1 :96: 

Volume  =  Height/30  *  (Al  +  (square  root(Al  *  A2))  +  A2)/1000 

where  Al  and  A2  equal  the  area  at  the  upper  and  lower  ends  of  the  snag.  In 
order  to  use  this  formula,  the  top  diameter  was  calculated  with  the  following 
formula  (Jim  Agee,  University  of  Washington,  personal  communication): 

Upper  Diameter  (D2)  =  Diameter  1  -  (1/96  *  Height) 

where  diameter  1  equals  DBH.  Snag  volumes  were  converted  to  biomass  (Mg/ha) 
with  the  same  formula  applied  to  logs. 

Forest  Floor  Depth,  Soil  Carbon,  and  Biomass 

In  July  and  August  1993,  sampling  was  conducted  to  estimate  mass  of  C  in 
the  forest  floor  and  mineral  soil  in  the  West  Twin  Creek  watershed.  In  5  of  the 
6  permanent  plots,  sampling  was  conducted  at  three  10-m  transects  along  ran- 
domly placed  radials  from  the  plot  center.  Forest  floor  depth  was  measured 
with  a  1-cm  diameter  steel  probe  at  0.25-m  intervals  along  the  transect.  The 
rod  was  inserted  into  the  forest  floor  until  an  increase  in  resistance  was  felt,  at 
which  point  it  was  assumed  the  transition  between  the  organic  horizon  and 
mineral  soil  was  reached.  Soil  cores  of  2-cm  diameter  were  taken  to  a  depth  'of 
40  cm  at  6  points  along  each  transect.  Each  core  was  separated  into  3  compo- 
nents: O  horizon,  0-4  cm  mineral  soil,  and  4-40  cm  mineral  soil.  The  forest 
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floor  samples  were  oven  dried  at  75°C  and  weighed  for  biomass.  Mineral  soil 
samples  were  air  dried,  ground,  and  analyzed  on  a  Leco  furnace  for  total  C 
concentration.  Percent  organic  matter  was  calculated  as  1.724  times  percent 
total  C. 

Bulk  density  samples  were  taken  with  a  soil  corer  (diameter  =  5.4  cm,  vol- 
ume =  137.4  cm-^)  in  all  permanent  plots  in  the  West  Twin  Creek  watershed. 
Three  small  soil  pits  were  dug  in  each  plot,  and  soil  cores  were  taken  from  the 
0-10  cm  depth  and  10-40  cm  depth.  These  samples  were  oven  dried  at  104°C 
for  48  hours.  Bulk  density  values  were  used  to  convert  percent  organic  matter 
in  mineral  soil  to  Mg/ha. 

Results  and  Discussion 

Biomass  in  the  West  Twin  Creek  and  Hoh  Lake  Watersheds 

The  average  biomass  of  live  trees  was  similar  in  the  West  Twin  Creek  water- 
shed (1044.7  Mg/ha,  Table  6)  and  the  Hoh  Lake  watershed  (1006.8  Mg/ha, 
Table  7)  despite  the  higher  elevation  at  Hoh  Lake.  Standard  deviations  were 
similar  in  both  watersheds.  Hoh  Lake,  however,  had  lower  total  CWD  biomass 
(182.1  Mg/ha)  than  West  Twin  Creek  (263.3  Mg/ha).  The  average  biomass  of 
snags  (standing  dead  trees)  in  the  West  Twin  Creek  watershed  was  171.0  Mg/ 
ha  (Table  6)  compared  to  144.5  Mg/ha  in  the  Hoh  Lake  Watershed  (Table  7). 
Log  biomass  averaged  92.3  Mg/ha  in  the  West  Twin  Creek  watershed  and  37.6 
Mg/ha  in  the  Hoh  Lake  watershed.  Shrub  and  herb  components  were  not  de- 
termined at  Hoh  Lake,  but  they  accounted  for  <  1  Mg/ha  in  the  West  Twin 
Creek  watershed  (Table  6).  Forest  floor  and  soil  organic  matter  were  also  not 
determined  at  Hoh  Lake,  but  they  contributed  considerably  (450.1  Mg/ha)  to 
the  total  ecosystem  biomass  of  1758.6  Mg/ha  in  the  West  Twin  Creek  water- 
shed (Table  6). 

Live  tree,  snag,  forest  floor  and  ecosystem  total  biomass  was  higher  in  the 
lower  West  Twin  Creek  watershed  than  in  the  upper  watershed  (Table  6).  Soil 
biomass  was  greater  in  the  upper  watershed,  while  log  biomass  was  similar  in 
the  upper  and  lower  watershed.  These  differences  reflect  the  cooler  environ- 
ment in  the  upper  watershed. 

Living  biomass  in  the  West  Twin  Creek  watershed  (616.0-1732.8  Mg/ha) 
was  similar  or  slightly  higher  than  that  reported  from  other  studies  in  the  Pa- 
cific Northwest  (Table  8).  For  example,  Agee  and  Huff  (1987)  reported  living 
biomass  values  exceeding  1000  Mg/ha  in  an  old-growth  stand  in  the  Hoh  River, 
and  Grier  (1976,  1978)  found  similar  vegetation  biomass  in  a  Sitka  spruce 
stand  on  the  Oregon  coast.  High  living  biomass  also  occur  in  old-growth  Dou- 
glas fir  stands  in  the  Washington  and  Oregon  Cascades  (Grier  and  Logan  1977; 
Franklin  and  Waring  1979;  Cole  and  Rapp  1980),  and  Franklin  and  Waring 
(1979)  reported  a  value  as  high  as  1318  Mg/ha  (Table  8). 

The  range  of  snag  biomass  (24.0-570.6  Mg/ha)  found  in  the  West  Twin 
Creek  watershed  was  generally  higher  than  previously  reported  in  the  Pacific 


36 


National  Park  Service  Monograph  2 


Table  6.  Average  biomass  distribution  and  standard  deviations  (Mg/ha)^  by 
ecosystem  component  in  the  upper;  lower,  and  entire  West  Twin  Creek  water- 
shed, 1985. 


Lower  watershed 

Upper  watershed 

Total  watershed 

Biomass  Component 

Average 

SD 

Average 

SD 

Average 

;      SD 

Tree 

Foliage 

29.1 

12.4 

32.5 

5.9 

30.8 

8.9 

Branches 

111.2 

26.8 

132.3 

44.6 

121.8 

34.9 

Roots 

306.2 

152.3 

216.1 

41.3 

261.1 

111.4 

Bole  +  Bark 

781.9 

387.3 

478.8 

89.6 

630.4 

301.3 

Total  Live  Tree 

1228.4 

566.3 

859.7 

75.9 

1044.1 

414 

Shrub 

Leaves 

0.2 

0.2 

0.1 

0 

0.1 

0.2 

Stems 

0.4 

0.4 

0.1 

0 

0.3 

0.3 

Herb 

Leaves 

0.2 

0.2 

0.1 

0 

0.1 

0.1 

Stems 

0.2 

0.2 

0.1 

0 

0.1 

0.1 

Total  Live  Biomass 

1229.4 

566.5 

860.1 

76 

1044.7 

414.3 

Detritus 

Standing  Dead  (snags) 

241.9 

289.6 

100.0 

47.6 

171.0 

201.2 

Fallen  Logs 

92.3 

10.5 

92.4 

6.7 

92.3 

7.9 

Total  CWD 

334.2 

279.1 

192.4 

44.8 

263.3 

194.9 

Dead  shrubs 

0.1 

0 

0 

0 

0.1 

0 

Dead  herbs 

0.4 

0.4 

0.1 

0 

0.4 

0.3 

Forest  floor 

97.1 

34.4 

84.9 

21.3 

89.8 

23.8 

Soil 

329.1 

67.4 

391.5 

83.8 

360.3 

76.3 

Ecosystem  Total 

1990.3 

1529.0 

1758.6 

''  Tree  biomass  was  estimated  with  equations  by  Gholz  et  al.  (1979).  Coarse  woody  debris 
dimensions  were  measured  empirically.  The  biomass  of  logs  was  calculated  with  formulas 
from  Hartman  et  al.  (1981),  and  the  biomass  of  snags  with  formulas  from  Husch  et  al.  (1972) 
and  Means  et  al.  (1992).  Density  values  of  coarse  woody  debris  were  taken  from  Means  et  al. 
(1992).  Values  were  corrected  for  slope. 


Northwest  (Table  8).  Log  biomass  (80.2-99.7  Mg/ha)  was  below  the  maxi- 
mum values  reported  from  other  studies  (Table  8).  Previous  estimates  of  CWD 
biomass  in  the  Hoh  River  watershed  range  from  26.4  to  84.5  Mg/ha  in  snags 
and  93.6  to  452.9  Mg/ha  in  logs  (Graham  and  Cromack  1982;  Agee  and  Huff 
1987).  In  old-growth  forests  on  Vancouver  Island,  Keenan  et  al.  (1993)  re- 
ported snag  biomass  ranging  from  1 9.6  to  1 64.3  Mg/ha  and  log  biomass  ranging 
from  175.1  to  397.5  Mg/ha.  Grier  (1976,  1978)  reported  log  biomass  values 
equaling  212  Mg/ha  in  an  old-growth  stand  at  Cascade  Head,  Oregon.  The 
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Table  7.  Average  biomass  distribution  and  standard  deviations  (Mg/ha)^  by 
ecosystem  component  in  the  Hob  Lake  watershed,  1985. 

Biomass  Component  Average  SD 

Tree 

Foliage  45.0  15.2 

Branches  102.5  37.7 

Roots  224.3  78.5 

Bole  +  Bark  635.0  264.7 

Total  Live  Tree  1006.8  383.7 

Detritus  ^ 

Standing  Dead  (snags)  144.5  178.8 

Fallen  Logs  37.6  29.4 

Total  CWD  182.1  188.6 

^Tree  biomass  was  estimated  with  equations  by  Gholz  et  al.  (1979).  Coarse  woody  debris  dimen- 
sions were  measured  empirically.  The  biomass  of  logs  was  calculated,  with  formulas  from  Hartman 
et  al.  (1981),  and  the  biomass  of  snags  with  formulas  from  Husch  et  al.  (1972)  and  Means  et  al. 
(1992).  Density  values  of  coarse  woody  debris  were  taken  from  Means  et  al.  (1992).  Values  were 
corrected  for  slope. 


great  variability  in  snag  and  log  biomass  in  these  coastal  stands  reflects  differ- 
ences in  fire  history,  presence  of  diseases  and  insects,  and  extent  of  windthrow. 

Biomass  estimates  of  forest  floor  (65.7-121.4  Mg/ha)  in  the  West  Twin  Creek 
watershed  (Table  8)  are  in  the  range  of  previously  reported  in  old-growth  for- 
ests of  the  Pacific  Northwest  (Grier  1976,  1978;  Cole  and  Rapp  1980;  Agee 
and  Huff  1987;  Boone  et  al.  1988;  Means  et.  al  1992),  but  soil  organic  matter 
was  sHghtly  higher  (281.4-471.0  Mg/ha).  Keenan  et  al.  (1993)  reported  much 
higher  forest  floor  biomass  on  Vancouver  Island  because  of  the  large  amounts 
of  decaying  wood  included  in  the  estimates.  Young  stands  of  Douglas  fir  and 
western  hemlock  in  coastal  Washington  had  much  lower  forest  floor  biomass 
values  (9.7-13.7  Mg/ha)  than  those  found  in  this  study,  while  soil  organic 
matter  was  similar  (213-532  Mg/ha;  Edmonds  and  Chappell  1994).  Differ- 
ences in  soil  organic  matter  among  sites  could  reflect  either  varied  sampling 
depths  or  different  climatic  conditions  and  decomposition  processes.  Forest 
floor  biomass  in  a  high  elevation  forest  in  Colorado  was  similar  to  the  esti- 
mates found  in  the  West  Twin  Creek  watershed  (Arthur  and  Fahey  1992). 
However,  biomass  estimates  of  organic  horizons  in  high  elevation  forests  in 
the  western  Olympic  mountains  were  lower  (29.8  Mg/ha;  Prichard  1996).  Or- 
ganic matter  in  mineral  soils  in  the  high  elevation  Olympics  was  slightly  lower 
( 1 77  Mg/ha)  than  that  found  in  the  West  Twin  Creek  watershed  (Prichard  1 996). 

The  living  biomass  estimates  in  the  Hoh  Lake  watershed  (479-1392.8  Mg/ 
ha)  are  greater  than  values  reported  from  other  high  elevation  sites  (Turner 
and  Singer  1976;  Grier  et  al.  1981;  Arthur  and  Fahey  1990;  Prichard  1996; 
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Table  8).  In  other  mountain  hemlock-Pacific  silver  fir  stands  in  the  Olympic 
Mountains,  ecosystem  living  biomass  averaged  501  Mg/ha  (Prichard  1996). 
Total  ecosystem  biomass  was  lower  in  high  elevation  forests  in  Colorado  (420 
Mg/ha;  Arthur  and  Fahey  1992).  In  Pacific  silver  fir  stands  in  the  Washington 
Cascades,  vegetation  biomass  was  450  Mg/ha  (Turner  and  Singer  1976;  Grier 
et  al.  1981),  which  is  the  low  end  of  the  vegetation  biomass  reported  in  this 
study.  The  high  living  vegetation  biomass  observed  in  the  Hoh  Lake  water- 
shed and  by  Prichard  (1996)  could  reflect  the  moister  conditions  of  forests  in 
the  western  Olympic  Mountains.  These  forests  receive  abundant  moisture  from 
storms  off  the  Pacific  Ocean  and  may  be  more  productive  than  high  elevation 
sites  in  the  Cascade  and  Rocky  Mountains.  In  contrast,  the  log  biomass  esti- 
mates in  Hoh  Lake  (10.9-84.0  Mg/ha)  are  in  the  same  range  as  other  studies 
(Grier  et  al.  1981;  Arthur  and  Fahey  1990;  Prichard  1996).  Log  and  snag  bio- 
mass combined  ranged  from  2.5  to  51.4  Mg/ha  in  the  Olympic  Mountains 
(Prichard  1996).  In  Pacific  silver  stands  in  the  Cascades,  log  biomass  was  75 
Mg/ha  (Grier  et  al.  1981);  it  was  52  Mg/ha  in  Colorado  (Arthur  and  Fahey 
1990). 

Clearly,  Pacific  Northwest  forest  ecosystems  are  dominated  by  large  trees, 
which  creates  very  high  ecosystem  biomass  values  composed  mostly  of  living 
trees  and  CWD.  Understory  plants  do  not  contribute  greatly  to  ecosystem  bio- 
mass. However,  below-ground  stores,  including  the  forest  floor,  mineral  soil, 
and  roots,  contribute  substantially  to  total  ecosystem  biomass.  In  contrast  to 
temperate  deciduous  forest  ecosystems,  conifer  forests  in  the  Pacific  North- 
west have  very  high  biomass  (Franklin  and  Waring  1979).  For  example,  in  a 
mature  aspen  stand  in  Michigan,  vegetation  biomass  ranged  from  1 18  to  154 
Mg/ha,  snag  biomass  ranged  from  6  to  10  Mg/ha,  and  log  biomass  ranged 
from  3  to  5  Mg/ha  (Alban  and  Perala  1992).  Similarly,  an  old-growth  oak 
forest  in  Kentucky  showed  much  smaller  snag  (5  Mg/ha)  and  log  biomass 
values  (14  Mg/ha;  Muller  and  Liu  1991)  than  those  found  in  this  study.  Mois- 
ture and  temperature  regimes  combined  with  litter  quality  directly  affect 
accumulations  of  detrital  material  by  influencing  decomposition  processes.  Thus, 
CWD  and  forest  floor  material  decomposes  more  rapidly  in  the  warm  deciduous 
forests  of  the  southeast  United  States  than  in  the  Pacific  Northwest  (Muller 
and  Liu  1991).  Long-lived  coniferous  species  combined  with  climatic  condi- 
tions that  slow  decomposition  lead  to  large  total  ecosystem  biomass 
accumulations. 

Variability  within  West  Twin  Creek  and  Hoh  Lake  Watersheds 

There  was  considerable  variability  in  biomass  within  the  West  Creek  water- 
shed. For  example,  the  average  total  live  tree  biomass  in  the  lower  West  Twin 
Creek  watershed  ( 1 228.4  Mg/ha)  was  much  greater  than  that  in  the  upper  West 
Twin  Creek  watershed  (859.7  Mg/ha;  Table  6).  Total  ecosystem  biomass  was 
also  higher  in  the  lower  watershed.  There  was  considerable  variability  in  live 
tree  biomass  among  the  various  plots  in  both  the  upper  and  lower  watershed  in 
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the  West  Twin  Creek  watershed  (Appendix  D).  In  the  West  Twin  Creek  water- 
shed, the  highest  total  biomass  (2393.2  Mg/ha)  was  in  plot  5  in  the  lower 
watershed  in  the  western  hemlock/salal  community  type.  The  highest  total  live 
tree  biomass,  however,  was  in  plot  1  in  the  western  hemlock/swordfern- 
Oregon  oxalis  community  (Appendix  D).  Log  biomass,  however,  was  not  very 
variable  and  only  ranged  from  80.2  to  99.7  Mg/ha  (Appendix  D).  In  contrast, 
snag  biomass  was  highly  variable  and  ranged  from  24.0  to  570.6  Mg/ha  (Ap- 
pendix D).  Forest  floor  biomass  ranged  from  65 .7  to  1 2 1 .4  Mg/ha.  Soil  organic 
matter  ranged  from  281.4  to  471.0  Mg/ha  (Appendix  D). 

There  was  also  a  lot  of  variability  in  biomass  in  the  various  plant  communi- 
ties in  the  Hoh  Lake  watershed  (Appendix  D).  The  highest  ecosystem  biomass 
was  found  in  plot  5  in  the  Pacific  silver  fir-mountain  hemlock/lily  community 
(1785.2  Mg/ha).  This  community  also  had  the  greatest  biomass  in  trees  and 
logs  of  all  the  plots.  The  lowest  biomass  value  in  the  Hoh  Lake  watershed  was 
found  in  plot  4  (501.4  Mg/ha)  in  the  Pacific  silver  fir/rhododendron  commu- 
nity (Appendix  D),  attributable  to  small  tree,  snag,  and  log  biomass  values. 
There  was  great  variability  in  biomass  of  both  logs  and  snags  among  the  Hoh 
Lake  watershed  plots.  Log  biomass  ranged  from  10.9  to  84.0  Mg/ha,  while 
snags  ranged  from  7.2  to  369.3  Mg/ha  (Appendix  D). 

Distribution  of  Biomass  Among  Ecosystem  Components 

Live  biomass  was  59.4%  of  the  total  ecosystem  biomass  in  the  West  Twin 
Creek  watershed.  CWD  made  up  15.0%  of  the  ecosystem  total  and  forest  floor 
and  soil  was  25.6%  (Table  6).  The  live  understory  made  up  only  0.03%  of  the 
total  ecosystem  biomass;  live  and  dead  understory  totaled  only  LI  Mg/ha. 

Snags  had  greater  biomass  than  logs  in  both  the  upper  and  lower  watershed, 
but  the  difference  was  greater  in  the  lower  watershed  (Table  6).  Mineral  soil 
contributed  higher  amounts  of  organic  matter  than  the  forest  floor  because  of 
the  larger  soil  volume  measured. 

No  fine  root  biomass  data  are  presented  in  Table  6  because  fine  root  biom- 
ass was  determined  only  in  plot  1  of  the  West  Twin  Creek  watershed,  where  it 
averaged  3.5  Mg/ha.  This  falls  within  the  range  of  values  from  the  Washington 
Cascades  (2.8-6.0  Mg/ha)  reported  by  Vogt  et  al.  (1980)  and  Vogt  et  al.  (1987). 
Fine  root  biomass  is  typically  greater  in  high  elevadon  forests,  due  to  nutrient 
poor  soils  (Grier  et  al.  1981).  In  an  old-growth  Pacific  silver  stand  in  the  Wash- 
ington Cascades,  fine  root  biomass  was  12.8  Mg/ha  (Grier  et  al.  1981).  This 
estimate  could  be  high  due  to  definifion  of  fine  roots  as  <5  mm,  rather  than  <2 
mm.  However,  Arthur  and  Fahey  (1992)  also  observed  high  fine  root  biomass  (9.0 
Mg/ha)  in  high  elevation  forests  in  Colorado.  Estimates  of  fine  root  biomass  in  a 
mature  hardwood  forest  in  the  Hubbard  Brook  Experimental  Forest  in  New  Hamp- 
shire were  comparable  to  values  in  this  study  and  ranged  from  4.3  to  4.7  Mg/ 
ha  depending  on  the  season  (Fahey  and  Hughes  1994). 

As  previously  menfioned,  understory  vegetation  represented  a  very  small 
component  of  the  total  ecosystem  biomass  in  the  West  twin  Creek  watershed. 
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The  biomass  of  understory  vegetation  was  similar  to  a  mature  forest  in  the 
Cedar  River  Watershed  (Turner  1975).  In  that  study,  vascular  understory  plants 
made  up  0.2  Mg/ha,  while  moss  c6ntributed  an  additional  1.0  Mg/ha  (Turner 
1975).  However,  the  understory  biomass  values  in  this  study  were  lower  than 
previous  estimates  in  old-growth  stands  in  Oregon  and  Washington  (3.7-9.7 
Mg/ha;  Johnson  et  al.  1982). 

We  found  the  average  forest  floor  mass  to  be  89.8  kg/ha  and  soil  values 
averaged  360.3  Mg/ha  (Table  6)  representing  25.6%  of  the  total  ecosystem 
biomass.  The  biomass  estimates  of  forest  floor  and  soil  in  the  West  Twin  Creek 
watershed  were  similar  to  values  found  in  other  coniferous  ecosystems  (Cole 
and  Rapp  1980;  Agee  and  Huff  1987;  Arthur  and  Fahey  1992).  Soil  organic 
matter  can  represent  a  substantial  component  of  the  ecosystem  total,  ranging 
from  9%  to  45%  of  the  total  ecosystem  biomass  and  averaging  33%  in  the 
Pacific  Northwest  (Cole  and  Rapp  1980).  Forest  floor  and  soil  organic  matter 
(to  60  cm  depth  in  mineral  soil)  accounted  for  34%  of  the  total  ecosystem 
biomass  in  an  old-growth  stand  in  the  Washington  Cascades  (Grier  et  al.  1981). 
Similarly,  Arthur  and  Fahey  (1992)  found  that  soil  organic  matter  made  up 
30%  of  the  ecosystem  total.  Prichard  (1996)  found  that  mineral  and  organic 
soils  contributed  18%  of  total  ecosystem  C,  while  roots  contributed  an  addi- 
tional 1 1%  in  forests  of  the  Olympic  Mountains. 

The  large  biomass  values  reported  from  the  West  Twin  Creek  and  Hoh  Lake 
watersheds  have  direct  implications  with  respect  to  global  change.  Productive 
forests  composed  of  large,  long-lived  trees,  such  as  those  found  on  the  west 
side  of  the  Olympic  Peninsula,  provide  a  sink  in  the  global  C  cycle  (Harmon  et 
al.  1990;  Dixon  et  al.  1994).  On  the  other  hand,  large  accumulations  of  CWD 
could  also  act  as  a  C  source  through  decomposition.  Downed  logs  respired 
from  2.7  to  8.3  g  CO2  m^/day  in  old-growth  stands  in  the  Hoh  River  water- 
shed; the  rate  depended  on  log  species,  decay  class,  and  season  (Marra  and 
Edmonds  1994).  Although  decomposition  of  CWD  may  reflect  a  substantial 
release  of  C  to  the  atmosphere,  old-growth  stands  with  large  biomass  accumu- 
lations in  living  vegetation  and  detritus  can  store  greater  amounts  of  global  C 
than  the  conversion  of  these  stands  to  fast-growing,  second-growth  stands 
(Harmon  et  al.  1990;  Dixon  et  al.  1994). 
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Chapter  5:  The  Hydrological  Cycle 
in  the  West  Twin  Creek  Watershed 

Because  water  has  a  vital  role  as  a  transporting  agent,  chemical  solvent,  and 
catalyst  in  ecosystems,  gathering  quantitative  data  on  hydrology  is  extremely 
important  in  understanding  the  biogeochemistry  of  watersheds  (Likens  et  al. 
1977;  Likens  and  Bormann  1995).  A  better  understanding  of  hydrologic  pro- 
cesses in  coastal  forests  on  the  Olympic  Peninsula  will  also  be  helpful  in 
developing  management  schemes  for  Pacific  salmon,  a  topic  of  great  interest 
in  the  Pacific  Northwest. 

There  is  a  considerable  amount  of  hydrologic  data  available  on  the  Olympic 
Peninsula.  The  USGS  operates  many  continuous  stream  gauging  stations,  but 
most  are  on  larger  rivers  and  streams.  There  are  few  on  smaller  streams  and 
only  a  handful  in  Olympic  National  Park  (Amerman  and  Orsborn  1987a  and 
1987b).  Our  study  helps  to  fill  this  data  gap.  Examination  of  annual  hydro- 
logic  budgets  and  hydrograph  patterns  in  small  watersheds  in  natural  systems 
can  provide  much  insight  into  the  functioning  of  forests  (Waring  et  al.  1981). 

One  of  the  problems  in  determining  annual  hydrologic  budgets,  however, 
is  determination  of  the  water  year.  An  ideal  water  year  is  that  successive  pe- 
riod that  most  consistently,  year  after  year,  gives  the  highest  correlation  between 
precipitation  and  stream  discharge.  In  the  Pacific  .Northwest,  the  water  year 
is  considered  to  start  on  1  October  after  the  dry  summer,  and  end  on  30 
September. 

Small  watershed  studies  enable  one  to  calculate  annual  evapotranspiration 
from  annual  precipitation  and  runoff  over  a  weir  at  the  base  of  the  watershed, 
assuming  that  the  geologic  base  of  the  watershed  is  watertight.  Evapotranspi- 
ration refers  to  the  total  water  loss  from  the  soil  as  a  result  of  evaporation  from 
soil  and  plant  surfaces  and  transpiration  from  the  trees.  The  living  components 
of  ecosystems,  particularly  forest  canopies,  strongly  influence  hydrologic  cycles 
in  ecosystems  (Likens  et  al.  1977;  Larson  1979;  SoUins  et  al.  1980),  and  inter- 
ception of  fog  by  an  old-growth  canopy  may  contribute  30%  of  the  annual 
input  to  the  forest  floor  (Harr  1982).  Watersheds  at  Hubbard  Brook  in  New 
Hampshire  have  relatively  impermeable  bedrock  and  little  deep  seepage.  At 
Hubbard  Brook  from  1956  to  1974,  an  average  of  129.5  cm  of  water  fell  on 
each  unit  area  of  watershed  annually  (Likens  et  al.  1977).  Of  this,  62%  became 
stream  discharge  and  the  remaining  38%  was  lost  as  evapotranspiration. 

In  the  Pacific  Northwest,  runoff  is  generally  a  higher  proportion  of  precipi- 
tation than  it  is  in  the  eastern  United  States  because  of  the  higher  precipitation, 
cooler  temperatures,  and  lower  evapotranspiration.  For  example,  Larson  (1979) 
found  annual  runoff  to  be  91%  of  precipitation  in  small  watersheds  in  the  Hoh 
River  valley  in  1975-76.  In  a  drier  year  (1976-77),  runoff  was  only  86%  of 
precipitation.  In  the  drier  Oregon  Cascades,  Sollins  et  al.  (1980)  found  that 
annual  runoff  from  a  small,  old-growth  watershed  was  65%  of  precipitation. 
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The  objectives  of  this  study  were  to  determine  the  hydrologic  budget  for  the 
pristine  old-growth  West  Twin  Creek  watershed  and  to  examine  the  hydrograph 
patterns.  Precipitation,  evapotranspiration,  and  stream  discharge,  and  how  they 
change  seasonally  and  annually  were  characterized.  Typically,  precipitation 
and  runoff  vary  greatly  from  year  to  year,  so  to  understand  an  ecosystem's 
hydrologic  flux,  a  long-term  record  is  needed  (Likens  et  al.  1977).  Therefore, 
we  conducted  our  study  from  1984  to  1994. 

Methods 

Precipitation 

Precipitation  was  collected  at  the  Hoh  Ranger  Station  to  the  east  of  the  study 
site  (Fig.  1).  Using  a  standard  rain  gauge,  precipitation  was  collected  daily 
from  January  1984  to  December  1994. 

Stream  Discharge 

Initially,  stream  discharge  was  measured  with  a  trapezoidal  weir  constructed 
in  September  1984.  In  May  1988,  the  weir  was  changed  to  a  3: 1  V-notch  (Fig. 
18).  This  change  was  made  to  measure  stream  stage  more  accurately  during 
summer  low  flows.  Stream  width  at  the  weir  ranged  from  about  1  to  2  m  de- 
pending on  the  season.  From  October  1984  through  October  1994,  the  stage 
and  discharge  at  the  weir  were  measured  every  two  weeks.  Stage  was  continu- 


Figure  18.  The  3: 1  V-notch  weir  at  West  Twin  Creek,  showing  the  stilling  well  in  the 
foreground  and  the  proportional  flow  sampler  in  the  background 
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ously  recorded  at  the  weir  by  a  Fluid  Data  Systems  float-driven  potentiometer 
and  data  were  collected  with  a  Campbell  Scientific  CR21  data  logger.  Stage- 
discharge  measurements  were  used  to  construct  a  rating  curve  for  the  weir. 
From  1984  to  April  1988,  the  discharge  equation  used  was  y  =  1.12-1-  1.83  x, 
and  from  May  1988  to  1994,  the  equation  used  was  y  =  1.18-1-  2.22  x,  where  y 
=  log  discharge  (in  cubic  feet  per  second)  and  x  =  log  stage  height  (in  feet). 
Annual  discharge  was  determined  in  centimeters  by  converting  discharge  to 
metric  units  and  assuming  the  area  of  the  watershed  was  58  ha.  During  1990- 
91,  stage  data  were  unreliable  and  not  used  in  this  analysis. 

The  following  hydrologic  budget  equation  (Swift  et  al.  1988)  was  used  on  a 
water-year  basis:  Annual  stream  discharge  or  runoff  =  annual  precipitation — 
annual  evapotranspiration  -i-  soil  storage  (assumed  to  be  zero  on  an  annual 
basis).  Fog  input  was  not  included  in  annual  precipitation. 

Throughfall 

Throughfall  is  precipitation  that  reaches  the  forest  floor  after  contacting  the 
tree  canopy.  Throughfall  was  collected  in  seven  throughfall  samplers  placed 
randomly  in  a  circular  0.01 -ha  plot  located  approximately  100  m  below  the 
weir.  The  collectors  consisted  of  20-cm  diameter  polyethylene  funnels  placed 
on  inverted  2-L  polyethylene  containers  with  the  bottom  removed.  The  funnel 
and  containers  were  placed  on  12-L  polyethylene  buckets.  A  2-cm  outside  di- 
ameter tube  was  attached  to  the  bottom  of  the  funnel,  and  throughfall  drained 
into  a  polyethylene  bag  in  the  bucket.  Two  layers  of  nylon  window  screen 
were  placed  over  the  top  of  the  funnel.  Samples  were  collected  every  four 
weeks  from  January  1992  to  December  1994. 

Results  and  Discussion 
Precipitation 

Monthly  precipitation  from  1984  to  1993  at  the  Hoh  Ranger  Station  is  shown 
in  Figure  19.  Precipitation  was  extremely  variable  from  month  to  month  and 
year  to  year.  Rain  fell  in  every  month  of  the  year,  but  there  was  a  distinctly 
seasonal  pattern  of  highest  precipitation  in  the  fall  and  winter  and  lowest  pre- 
cipitation in  the  summer.  Annual  precipitation  ranged  from  a  low  of  267  cm 
during  1992-93  to  475  cm  during  1990-91  (Table  9).  Larson  (1979)  recorded 
an  even  higher  annual  precipitation  of  537  cm  during  1975-76  in  a  nearby 
watershed.  The  high  precipitation  during  1990-91  was  due  to  intense  storms 
in  the  fall  (Fig.  19).  Fog  does  occur  in  the  watershed,  but  no  estimates  of  fog 
drip  as  a  contribution  to  precipitation  reaching  the  forest  floor  were  made. 
Harr  (1982)  estimated  that  fog  drip  contributed  as  much  as  30%  of  the  annual 
precipitation  in  the  Bull  Run  watershed  near  Portland  in  the  Oregon  Cascades. 

Stream  Discharge 

Annual  stream  discharge  ranged  from  175  cm  during  1986-87  to  347  cm 
during  1989-90  (Table  9).  Stream  hydrographs  can  be  separated  into  two  main 
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components:  direct  runoff  (or  quickflow)  and  baseflow,  which  represents  the 
groundwater  contribution  (Gordon  et  al.  1992).  The  stream  discharge  hydrograph 
for  January  1992  to  December  1993  is  shown  in  Figure  20.  Large  peaks  are 
related  to  storms  in  February  1992,  January  1993,  May  1993,  and  December 
1993.  The  stream  at  West  Twin  Creek  flows  throughout  the  year,  but  it  almost 
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Figure  19.  Monthly  rainfall  at  the  Hoh  Ranger  Station  (National  Acid  Deposition 
Program/National  Trends  Network  1994),  1984-1993. 
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Figure  20.  Stream  discharge  hydrograph  at  the  weir  on  the  West  Twin  Creek, 
1993. 


1992- 


stopped  flowing  in  August  and  September  1993.  Summer  flow  in  1992  was 
much  greater  than  in  1993.  Figure  21  shows  the  detailed  hydrograph  for  the 
storm  on  25  January  1993.  Baseflow  was  about  4500  mVday.  During  the  storm, 
stream  discharge  rose  quickly  to  about  10  000  m^/day  and  fell  quickly  to  about 
half  the  peak  discharge,  or  7000  m-^/day.  It  then  took  several  days  to  return  to 
the  previous  discharge  rate. 

Figure  22  shows  the  hydrograph  for  a  storm  on  8  September  1992,  during 
the  dry  season.  Baseflow  was  about  2700  mVday.  During  the  storm,  discharge 
increased  quickly  to  about  4000  mVday,  and  like  the  winter  storm,  it  fell  rap- 
idly to  about  half  the  peak  discharge.  However,  recovery  to  initial  levels  was 
slower,  and  elevated  flows  were  still  recorded  two  weeks  after  the  storm. 

Hydrograph  shape  is  influenced  by  many  factors  including  the  slope  and 
length  of  the  upland  surface,  catchment  shape,  microtopography,  permeability 
and  moisture  content  of  the  soil,  subsurface  geology,  vegetation  cover,  drain- 
age density,  channel  characteristics,  and  storm  patterns  (Gordon  et  al.  1992). 
The  West  Twin  Creek  hydrograph  has  a  rapid  response  and  rises  quickly.  Such 
hydrographs  are  typical  of  steeply  sloping  watersheds  (Swift  et  al.  1988).  Short, 
wide  catchments  cause  streams  to  rise  faster  than  longer  narrow  ones  because 
of  short  travel  times.  Catchment  factors,  however,  may  not  always  correlate 
well  with  hydrograph  shape  because  they  do  not  reflect  orientation  of  the  catch- 
ment to  storm  tracks  that  may  move  up,  down,  or  across  a  stream.  In  West 
Twin  Creek  the  storms  typically  move  from  west  to  east  across  the  stream. 
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Figure  21.  Hydrograph  of  one  storm  (25  January  1993)  during  the  wet  season,  demon- 
strating a  typical,  high- volume  stream  discharge  at  the  weir  on  West  Twin  Creek. 
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Figure  22.  Hydrograph  of  one  storm  (8  September  1992)  during  the  dry  season,  demon- 
strating a  typical,  low-volume  stream  discharge  at  the  weir  on  West  Twin  Creek. 
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Vegetation  cover  and  land  use  also  strongly  influence  hydrograph  shape 
(Gordon  et  al.  1992;  Jones  and  Grant  1996;  Wemple  et  al.  1996).  Old-growth 
forests  moderate  peak  discharge  and  reduce  the  effect  of  storms.  Jones  and 
Grant  (1996)  illustrated  this  in  a  study  in  the  western  Oregon  Cascade  Moun- 
tains, where  forest  removal  and  road  construction  changed  the  shape  of  the 
hydrograph.  Complete  clearcutting  of  the  old-growth  forest  not  only  increased 
peak  discharge,  but  also  increased  total  storm  volumes  compared  to  the  uncut 
forest.  The  suppression  of  evapotranspiration  after  removal  of  deep-rooted 
conifers  led  to  increased  subsurface  moisture  storage  for  many  years  after  cut- 
ting, even  though  shallow  soil  moisture  returned  to  pretreatment  levels  one 
year  after  harvesting.  Watersheds  with  just  6%  of  their  surface  in  roads  had 
much  higher  peak  storm  discharges,  but  the  hydrograph  fell  rapidly  in  contrast 
to  that  in  clearcut  watersheds. 

Hydrologic  Budget 

A  large  amount  of  the  precipitation  falling  on  the  West  Twin  Creek  water- 
shed leaves  the  watershed  via  the  stream.  From  1985  to  1993,  discharge  ranged 
from  54%  to  89%  (Table  9).  This  is  typical  for  Northwest  forests.  Larson  (1979) 
found  discharge  to  range  from  86%  to  91%  in  nearby  old-growth  watersheds 
in  the  Hoh  River  Valley.  In  the  Cascades  near  Seattle,  he  found  discharge  to 
range  from  73%  to  79%  of  precipitation  in  a  watershed  consisting  of  10-20 
and  40-50-year-old  Douglas  fir  stands  with  some  residual  old  growth  (Larson 
1973).  In  the  eastern  United  States,  Swift  et  al.  (1988)  also  found  that  water- 
sheds at  Coweeta  with  steep  slopes  and  shallow  soils  yielded  75%  of  rainfall 
as  runoff. 

Discharge  was  positively  related  to  precipitation  in  the  six  years  for  which 
we  had  reliable  data  (Fig.  23,  r^  -  0.72,  p  =  0.06).  Fog  inputs,  however,  are  not 
taken  into  account  in  Figure  23,  and  as  mentioned  previously,  fog  can  contrib- 
ute as  much  as  30%  of  the  annual  precipitation  (Harr  1982).  A  data  point 
showing  discharge  and  precipitation  from  a  small  watershed  in  the  Oregon 
Cascades  (Sollins  et  al.  1980)  is  also  included  in, Figure  23,  and  it  falls  very 
close  to  the  regression  line.  At  Hubbard  Brook  in  New  Hampshire,  annual 
stream  discharge  was  also  highly  correlated  with  annual  precipitation  (Likens 
etal.  1977). 

Evapotranspiration  also  varied  from  year  to  year  and  ranged  from  36  to  150 
cm  (Table  9)  or  1 1  %  to  46%  of  precipitation.  Figure  24  shows  average  monthly 
precipitation,  potential  evapotranspiration,  and  actual  evapotranspiration  from 
1913  to  1965  for  Forks,  Washington,  located  50  km  to  the  northwest  of  our 
study  site.  There  was  a  water  surplus  except  for  July  and  August,  and  there 
was  very  little  water  deficit.  At  Forks,  actual  evapotranspiration  averaged  20% 
of  precipitation  from  1913  to  1965  (Phillips  and  Donaldson  1972)  and  was 
similar  to  evapotranspiration  at  the  Quinault  Ranger  Station,  which  is  located 
in  an  old-growth  forest  to  the  south  of  our  site  (Phillips  and  Donaldson  1972). 
At  West  Twin  Creek,  evapotranspiration  calculated  from  precipitation  and 
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Figure  23.  Annual  stream  discharge  at  the  West  Twin  Creek  weir  and  precipitation  at 
the  West  Twin  Creek  from  Water  Year  85-86  to  93-94. 
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Figure  24.  The  relations  among  average  monthly  precipitation,  potential  evapotrans- 
piration,  and  actual  evapotranspiration  (in  cm)  measured  in  Forks,  Washington, 
1913-1965. 
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stream  discharge  averaged  112%  of  precipitation  (Table  9),  which  is  similar 
to  actual  evapotranspiration  calculated  for  Forks  from  1913  to  1965. 

Relations  between  Throughfall  and  Precipitation 

Although  precipitation  varied  from  year  to  year,  the  percentage  of  annual 
precipitation  reaching  the  forest  floor  in  throughfall  was  remarkably  constant, 
ranging  from  56.4%  to  63.1%  from  1992  to  1994  and  averaging  60.4%  (Table 
10).  Similar  percentages  were  observed  by  Fahey  et  al.  (1988)  and  Blew  et  al. 
(1993)  in  drier  Rocky  Mountain  forests.  On  a  monthly  basis,  however,  the 
percentage  of  throughfall  to  precipitation  ranged  widely  from  0  to  nearly  89%. 
Monthly  throughfall  is  plotted  against  precipitation  in  Figure  25,  revealing 
that  as  precipitation  increased,  throughfall  also  increased,  but  only  to  a  monthly 
precipitation  of  about  50  cm.  If  monthly  precipitation  totals  exceeded  50  cm, 
throughfall  remained  constant  at  about  30  cm.  It  is  presumed  that  a  large  pro- 
portion of  the  water  reaching  the  forest  floor  is  via  stemflow  when  precipitation 
is  high,  and  during  large  storm  events,  stemflow  was  commonly  observed. 

Fahey  et  al.  (1988)  proposed  that  throughfall  could  be  predicted  using  the 
following  equation:  throughfall  =  0.634  x  rainfall^ ^^^  Similarly,  Blew  et  al. 
(1993)  suggested  that  throughfall  could  be  predicted  using  this  equation: 
throughfall  =  0.558  x  precipitation' ^^^.  The  range  of  their  data,  however,  was  for 
precipitation  events  <10  cm.  These  equations  would  vastly  overpredict  throughfall 
for  monthly  precipitation  >  10  cm  in  the  West  Twin  Creek  watershed. 
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Figure  25.  Relation  between  monthly  throughfall  and  precipitation  (in  cm)  in  the  West 
Twin  Creek  watershed. 
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Table  10.  Estimated  monthly  precipitation  and  throughfall  (in  cm)  and  the 
percentage  of  precipitation  recovered  as  throughfall  in  the  West  Twin  Creek 
watershed,  1992-1994. 


Date 

Precipitation 

Throughfall 

Percent  Throughfall 

01-13-92 

28.4 

17.8 

62.7 

02-10-92 

64.0 

29.5 

46.1 

03-11-92 

27.4 

16.3 

59.5 

04-08-92 

7.7 

5.1 

66.2 

05-06-92 

33.3 

20.0 

60.1 

06-03-92 

4.8 

4.0 

83.3 

07-01-92 

10.7 

5.9 

55.1 

07-29-92 

1.8 

0.6 

33.3 

08-26-92 

5.9 

4.5 

76.3 

09-23-92 

14.0 

11.4 

81.4 

10-21-92 

20.7 

18.4 

88.9 

11-18-92 

35.6 

24.7 

69.4 

12-16-92 

39.3 

27.1 

69.0 

Total  1992 

293.6 

185.3 

63.1 

01-02-93 

ND^ 

ND 

ND 

02-09-93 

ND 

22.7 

ND 

03-09-93 

15.1 

9.9 

65.8 

04-07-93 

23.6 

17.0 

72.0 

05-02-93 

42.6 

23.1 

54.2 

06-02-93 

22.5 

12.6 

56.0 

06-30-93 

14.2 

5.0 

35.2 

07-30-93 

6.7 

1.7 

25.3 

08-26-93 

5.9 

4.7 

79.7 

09-22-93 

0.5 

0.0 

0.0 

10-21-93 

6.2 

2.3 

37.1 

11-17-93 

17.8 

15.1 

84.8 

12-15-93 

60.0 

29.8 

49.7 

Total  1993 

215.0 

121.2 

56.4 

01-12-94 

32.8 

27.0 

82.3 

02-07-94 

12.6 

10.4 

82.5 

03-07-94 

60.0 

30.6 

51.0 

04-03-94 

18.6 

12.8 

68.8 

05-03-94 

24.6 

14.8 

60.2 

05-31-94 

13.5 

7.7 

57.0 

06-28-94 

16.5 

10.2 

61.8 

07-26-94 

4.0 

3.2 

79.4 

08-24-94 

3.8 

1.1 

28.9 

09-21-94 

12.5 

8.7 

69.5 

10-21-94 

6.6 

3.7 

56.1 

11-22-94 

53.1 

29.6 

55.7 

12-16-94 

36.6 

22.5 

61.5 

Total  1994 

295.2 

182.3 

61.7 

ND  =  Not  determined. 
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Chapter  6:  Precipitation  Chemistry 
along  a  Gradient  Inland 
from  the  Pacific  Ocean 

The  oceans  serve  as  a  source  of  aerosols  necessary  for  cloud  formation 
(Charlson  et  al.  1987)  and  precipitation.  These  aerosols  contain  several  ions 
(Na,  CI,  SO4  and  Mg)  that  are  the  main  constituents  of  precipitation  in  areas 
that  have  a  marine  influence  (Bormann  et  al.  1989;  McDowell  et  al.  1990).  In 
such  areas,  the  relative  amounts  in  precipitation  of  these  and  other  solutes  of 
sea  salt  origin  remain  quite  constant  and  similar  to  amounts  in  seawater.  For 
example,  several  studies  have  suggested  that  the  ratio  of  CI  to  Na  in  rainwater 
collected  within  several  hundred  kilometers  of  the  ocean  is  the  same  as  in 
seawater  (Junge  and  Werby  1958;  Vong  1989).  Likewise,  the  marine  contribu- 
tion to  SO4  and  Ca  deposition  can  be  estimated  from  their  ratios  with  Na  or  CI 
in  sea  salt  (Junge  1963;  Schwartz  1988). 

Marine  influences  on  precipitation  chemistry  are  especially  important  in  the 
northwest  United  States  due  to  the  dominant  westerly  flows  that  bring  air  masses 
inland  from  the  Pacific  Ocean.  Several  studies  have  investigated  marine  influ- 
ences on  the  chemistry  of  atmosphere  and  precipitation  along  the  Pacific  coast 
of  Washington  state  (Ludwick  et  al.  1977;  Vong  et  al.  1985;  Andreae  et  al. 
1988;  Quinn  et  al.  1988;  Vong  et  al.  1988a  and  1988b;  Vong  1989;  Edmonds  et 
al.  1991).  Studies  of  precipitation  and  atmospheric  chemistry  on  the  Olympic 
Peninsula  generally  have  had  one  of  three  designs:  investigations  of  on-shore 
precipitation  that  included  analysis  of  back-trajectories  of  marine  air  masses 
(Ludwick  et  al.  1977;  Vong  et  al.  1988a  and  1988b;  Quinn  et  al.  1988),  estima- 
tion of  chemical  inputs  to  ecosystems  using  multiple  sampling  sites  that  were 
considered  equivalent  enough  to  be  used  as  replicates  (Vong  1989;  Edmonds  et  al. 
1991),  and  investigations  of  precipitation  chemistry  at  several  coastal  sites  to  as- 
sess natural  or  background  concentrations  of  pollutant  constituents  (Vong  1989). 

Several  studies  have  demonstrated  the  influence  of  oceans  on  the  chemistry 
of  precipitation  falling  on  land  masses  and  changes  in  precipitation  chemistry 
across  land  masses.  Junge  (1958)  and  Junge  and  Werby  (1958)  reported  on  the 
distribution  and  sources  of  NH^,  CI,  Na,  K,  Ca,  and  SO4  concentrations  in 
precipitation  over  the  continental  United  States.  McColl  (1982)  and  McColl  et 
al.  (1982)  examined  the  relations  between  distance  from  the  Pacific  Ocean  and 
chemical  composition  of  wet  and  dry  deposition.  However,  such  studies  in- 
volving spatial  distribution  of  precipitation  chemistry  had  resolutions  of  only 
tens  to  hundreds  of  kilometers.  Also,  they  were  concerned  with  sources  and 
long  range  transport  of  atmospheric  constituents  across  land  masses  and  oceans 
rather  than  examination  of  changes  in  precipitation  chemistry  associated  with 
ocean/continent  interfaces.  Little  is  known  about  the  changes  in  precipitation 
chemistry  across  relatively  short  distances  inland  from  the  coast. 
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Our  study  assessed  changes  in  precipitation  chemistry  at  the  Pacific  Ocean/ 
Olympic  Peninsula  interface  on  a  finer  spatial  scale  than  has  previously  been 
reported.  Objectives  were  to  determine  the  effect  of  distance  inland  from  the 
Pacific  Ocean  on  concentrations  of  major  cations  and  anions  in  bulk  precipita- 
tion, to  determine  spatial  and  seasonal  variation  in  deposition  rates  of  major 
cations  and  anions  to  coastal  temperate  rain  forests,  and  to  examine  oceanic 
influences  on  precipitation  chemistry  on  the  Pacific  coast. 

Methods 

The  precipitation  transect  was  located  in  the  Hoh  River  valley  immediately 
west  of  the  Olympic  National  Park  boundary  (32  km  from  the  Pacific  coast) 
and  extended  nearly  to  the  coast  (Fig.  1).  Four  sampling  sites  were  located  in 
clearcuts  that  were  large  enough  to  have  no  sampling  interference  from  the 
forest  canopy.  Sampling  sites  were  located  4,  13,  24,  and  3 1  km  from  the  ocean 
at  elevations  of  120,  100,  170,  and  200  m,  respectively.  At  each  sampling  loca- 
tion there  were  two  bulk  precipitation  collectors  and  a  Belfort  standard 
raingauge.  Each  bulk  collector  consisted  of  a  20-cm  diameter  polyethylene 
funnel  with  a  netted  top  to  exclude  debris.  The  precipitation  drained  from  the 
funnel  into  a  polyethylene  bag.  Samples  were  poured  into  triple-rinsed  1-L 
polyethylene  bottles  and  transported  to  the  Analytical  Laboratory  at  the  Col- 
lege of  Forest  Resources,  University  of  Washington,  within  48  h  of  collection 
and  stored  at  4°C.  Samples  used  in  this  analysis  were  collected  between  Feb- 
ruary 1989  and  September  1992.  Samples  were  collected  at  two-week  intervals 
through  June  of  1990  and  then  at  four- week  intervals  thereafter.  Sampling  at 
the  intermediate  sites  (13  and  24  km  from  ocean)  was  discontinued  in  October 
1991. 

Electrical  conductivity,  pH,  and  alkalinity  (reported  as  HCO3)  were  deter- 
mined on  non-filtered  samples  within  a  few  days  after  arrival  at  the  laboratory. 
Electrical  conductivity  was  determined  with  a  YSI  Model  3 1  Conductivity 
Bridge  (Yellow  Springs  Instrument  Co.,  Yellow  Springs,  OH)  and  corrected  to 
25°C.  pH  was  measured  with  a  Radiometer  PHM85  pH  meter  (Radiometer, 
Copenhagen,  Denmark).  Alkalinity  was  determined  by  titration  (Rhoades  1982) 
to  an  end  point  of  pH  5. 

Remaining  samples  were  filtered  through  Whatman  GF/A  filters  (1.6  |nm 
pore  size),  stored  at  4°C,  and  analyzed  within  a  month  after  arrival  at  the  labo- 
ratory. Solutions  were  analyzed  for  Ca,  Mg,  K,  and  Na,  using  a  IL  95 1  Atomic 
Absorption  Analyzer  (Instrumentation  Laboratory,  Wilmington,  MA).  Ammo- 
nium was  measured  by  the  automated  phenate  method  (EPA  Method  350.1) 
with  a  Technicon  Autoanalyzer  II  (Technicon,  Tarrytown,  NY).  Sulfate,  CI, 
NO3,  and  PO4  were  determined  with  a  Dionex  2010  Ion  Chromatograph  (AS4A 
and  AG4A  columns;  Dionex,  Sunnyvale,  CA).  Macro-Kjeldahl  digests  were 
done  on  all  samples  having  sufficient  volume  and  those  digests  were  analyzed 
for  N  and  P. 
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Data  screening  was  a  three  step  process  based  on  the  procedures  of  Vong 
(1989).  First,  samples  were  eliminated  if  field  notes  made  at  the  time  of  sample 
collection  indicated  any  obvious  problems  with  sample  quality.  The  next  screen- 
ing step  excluded  those  samples  with  a  charge  balance  error  greater  than  40%. 
Charge  balance  error  was  calculated  as  the  difference  between  anions  and  cat- 
ions divided  by  the  sum  of  anions  and  cations  multiplied  by  100.  Finally,  samples 
were  screened  on  the  basis  of  the  Cl:Na  ratio.  To  determine  amounts  of  SO4 
and  Ca  in  precipitation  that  were  in  excess  of  that  derived  from  sea  salt  (excess 
SO4  and  excess  Ca)  based  on  Na  ratios  in  sea  salt,  we  compared  the  Cl:Na 
ratio  in  precipitation  to  that  found  in  sea  salt.  Precipitation  collected  within 
several  hundred  kilometers  of  the  ocean  generally  has  the  same  Cl:Na  ratio  as 
seawater  (Junge  1963).  Samples  with  Cl:Na  molar  ratios  that  deviated  from 
the  sea  salt  ratio  of  1.17  (1.8  mass  ratio)  by  of  factor  greater  than  2  were  ex- 
cluded from  the  data  set. 

Concentrations  of  excess  SO^  and  excess  Ca  were  calculated  based  on  the 
ratios  of  Na  with  SO4  and  Ca  in  sea  salt  such  that: 

Excess  SO4  =  SO4  -  (Na  •  0.1 17)  and 
Excess  Ca  =  Ca  -  (Na  •  0.0434) 

where  Excess  SO4,  SO^,  Excess  Ca,  Ca,  and  Na  are  ion  concentrations  in  |Lieq/ 
L.  Correlation  coefficients  (Wilkinson  1990)  were  calculated  to  investigate 
relations  between  concentrations  of  measured  chemical  parameters  and  pre- 
cipitation amount. 

The  effect  of  distance  from  the  ocean  on  ion  concentrations  in  precipitation 
was  investigated  using  volume-weighted  mean  concentrations  from  each  of 
the  two  collectors  at  each  site.  Volumes  were  based  on  precipitation  depth  in 
Belfort  rain  gauges  at  each  site.  Statistical  analysis  for  the  effect  of  distance 
from  the  ocean  used  one-way  analysis  of  variance  (ANOVA;  Wilkinson  1990). 

Mean  deposition  rates  were  calculated  by  totaling  deposition  for  each  of  13 
four- week  periods  each  year  and  then  averaging,  equivalent  periods  from  each 
year.  Annual  deposition  rates  at  each  site  were  determined  by  the  13  four-week 
means.  This  method  of  calculation  was  necessary  due  to  missing  data  for  sev- 
eral four-week  periods  that  did  not  allow  for  calculation  of  annual  totals  for 
each  individual  year.  Unfortunately,  this  calculation  allowed  only  for  an  esti- 
mate of  the  average  annual  deposition  and  did  not  allow  for  any  estimation  of 
the  variability  associated  with  annual  deposition.  Wet  deposition  data  from  the 
Hoh  Ranger  Station  (National  Atmospheric  Deposition  Program/National 
Trends  Network  1992)  was  used  for  comparison  to  the  data  collected  in  this 
study  and  as  an  extension  of  the  transect  to  41  km  inland  (Fig.  1). 

Screening  of  data  based  on  field  notes  and  charge  balance  resulted  in  the 
removal  of  36  of  a  total  of  327  samples.  An  additional  41  samples  were  re- 
moved based  on  Cl:Na  ratio  screening.  The  remaining  250  samples  were 
included  in  the  analysis. 
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Results 


Correlations  Between  Constituents 

There  were  no  significant  positive  correlations  between  precipitation  depth 
and  any  chemical  parameter,  but  there  was  a  negative  correlation  with  excess 
SO4  (Table  1 1).  Strong  correlations  were  found  between  CI,  SO^,  Mg,  and  Na 
(Table  11).  Electrical  conductivity  was  most  closely  correlated  with  concen- 
trations of  those  same  ions  (Table  11). 

Correlations  between  TKN  and  TKP,  PO^,  and  NH^  were  significant  (Table 
11).  These  correlations  appeared  to  be  due  to  influential  outliers.  Outliers  were 
also  reflected  in  the  high  standard  deviations  associated  with  these  variables 
(Table  12),  especially  at  the  sites  13  and  24  km  from  the  ocean.  These  outliers 
may  represent  the  input  of  organic  contaminants  to  the  precipitation  collector 
during  a  few  sampling  periods.  This  suggests  that  the  source  of  PO^  and  NH^ 
on  some  occasions  may  have  been  from  mineralization  of  organic  material  in 
the  collector.  As  a  result,  the  data  presented  here  may  be  overestimates  of  the 
atmospheric  deposition  of  these  constituents. 

Distance  Effects 

Significant  effects  (p  <  0.05)  of  proximity  to  the  Pacific  Ocean  were  found 
for  electrical  conductivity  and  for  concentrations  of  H,  CI,  SO4,  Mg,  Na,  and 
excess  Ca  (non-sea-salt  Ca;  Table  12).  (ANOVA  indicated  a  significant  dis- 
tance effect  for  H  (p  =  0.050),  but  the  conservative  nature  of  Tukey's  HSD  to 
control  experiment-wise  error  rate  rather  than  comparison-wise  error  rate,  failed 
to  demonstrate  significant  differences  between  any  of  the  group  means.)  Con- 
centrations of  CI,  SO4,  Mg,  and  Na  were  highest  nearest  to  the  coast.  These 
higher  ion  concentrations  were  reflected  in  significantly  higher  electrical  con- 
ductivity in  precipitation  falling  closest  to  the  coast.  Concentrations  of  other 
ions  did  not  differ  significantly  at  different  distances  from  the  coast.  However, 
NO3,  H,  and  K  tended  to  be  highest  at  the  coastal  site. 

Annual  deposition  of  Na,  K,  Mg,  Ca,  excess  Ca,  NO3,  and  CI  tended  to  be 
lower  at  the  collection  site  farthest  from  the  ocean  (Table  13).  Sea  salt  Ca 
(total  minus  excess)  deposition  decreased  from  0.64  at  4  km  to  0.29  kg/ha/year 
at  31  km  from  the  ocean.  Even  though  total  SO^  deposition  remained  rela- 
tively unchanged  as  distance  from  the  ocean  increased  (12.6  kg/ha/year  at  4 
km  and  12.2  kg/ha/year  at  3 1  km),  the  amounts  attributable  to  sea  salt  and  non- 
sea-salt  sources  tended  to  change.  Excess  SO^  (non-sea-salt  SO^)  deposition 
increased  from  3.81  at  4  km  to  4.59  kg/ha/year  at  3 1  km  from  the  ocean.  Depo- 
sition of  SO4  of  sea  salt  origin  (total  minus  excess)  decreased  from  8.79  kg/ha/ 
year  at  4  km  to  7.61  kg/ha/year  at  3 1  km  from  the  ocean. 

Deposition  of  total  P  (0.21  kg/ha/year)  and  PO^  (0.05  kg/ha/year)  remained 
unchanged  as  distance  from  the  ocean  increased  (Table  13).  Nitrate  (NO3)  in 
bulk  deposition  tended  to  decrease  from  1 .86  kg/ha/year  at  4  km  to  1 .58  kg/ha/ 
year  at  31  km  from  the  coast.  Total  Kjeldahl  N  deposition  (TKN)  also  ap- 
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Table  12.  Volume-weighted  mean  concentrations  of  ions  and  electrical  con- 
ductivity (EC)  in  bulk  precipitation  along  a  transect  inland  from  the  Pacific 
Ocean,  1989-1992.  Numbers  in  parentheses  are  standard  deviations. 


Distance  (km) 

Analyte 

4 

13 

24 

31 

4ia 

EC  (mS/m) 

1.02a 

0.66ab 

0.47b 

0.51b 

0.60 

(0.03) 

(0.06) 

(0.21) 

(0.02) 

ILieq/L'' 

HCO3 

5.28a 

6.68a 

3.93a 

3.99a 

ND'^ 

(0.77) 

(0.51) 

(1.47) 

(0.18) 

CI 

55.9a 

35.0b 

25.4b 

28.9b 

27.3 

(1.7) 

(0.6) 

(9.9) 

(0.3) 

NO3 

0.95a 

0.56a 

0.42a 

0.49a 

1.21 

(0.08) 

(0.17) 

(0.20) 

(0.11) 

PO4 

0.06a 

1.08a 

0.02a 

0.02a 

ND 

(0.08) 

(1.46) 

(0.02) 

(0.12) 

SO4 

9.11a 

6.23ab 

4.40b 

4.92b 

5.81 

(0.34) 

(0.94) 

(1.73) 

(0.47) 

Excess  SO4 

2.81a 

1.88a 

1.37a 

1.58a 

ND 

(0.28) 

(0.65) 

(0.38) 

(0.43) 

H 

5.07a 

3.01a 

2.75a 

3.15a 

4.02 

(0.32) 

(0.35) 

(1.07) 

(0.12) 

NH^ 

0.08a 

1.85a 

0.43a 

0.03a 

0.98 

(0.02) 

(2.32) 

(0.10) 

(0.02) 

Ca 

5.51ab 

5.89a 

3.93ab 

3.32b 

1.84 

(0.01) 

(0.78) 

(1.05) 

(0.15) 

Excess  Ca 

4.16a 

4.54a 

3.19ab 

1.55b 

ND 

(0.08) 

(0.61) 

(0.77) 

(0.11) 

K 

2.95a 

1.25a 

1.70a 

1.12a 

0.5 

(1.93) 

(0.28) 

(0.36) 

(0.17) 

Mg 

7.38a 

3.51b 

2.22b 

2.55b 

5.25 

(0.10) 

(0.79) 

(0.70) 

(0.17) 

Na 

53.9a 

38.5ab 

26.3b 

28.8b 

22.9 

(0.53) 

(1.00) 

(11.6) 

(0.26) 

'""  Data  from  Hoh  Ranger  Station  (National  Atmospheric  Deposition  Program/National  Trends 
Network  1992). 

^  Numbers  followed  by  the  same  letter  within  a  row  are  not  significanlly  different  (/>  <  0.05). 
"  ND  =  Not  determined. 
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Table  13.  Average  of  annual  deposition  rates  (kg/ha/yr)  at  two  distances  from 
the  Pacific  Ocean,  1989-1992. 

Distance  (km) 


Analyte  4.00  31.00 

CI  57.20  49.30 

NO3  1.86  1.58 

PO4  0.05  0.05 

SO4  12.60  12.20 

Non-sea-salt  SO4  3.81  4.59 

Sea  salt  SO4  8.79  7.61 

H  0.16  0.19 

NH4  0.04  0.04 

Ca  2.86  2.14 


57.20 

1.86 

0.05 

12.60 

3.81 

8.79 

0.16 

0.04 

2.86 

2.22 

0.64 

3.03 

2.41 

36.10 

3.39 

0.21 

3.42 

78.00 

Non-sea-salt  Ca  2.22  1.85 

Sea  salt  Ca  0.64  0.29 

K  3.03  2.23 

Mg  2.41  0.88 

Na  36.10  31.20 

TKN  3.39  3.08 

TKP  0.21  0.21 

Total  N  3.42  3.11 

Precipitation  (cm)  278.00  351.00 


peared  to  decrease  as  distance  from  the  ocean  increased  (3.39  kg/ha/year  at  4 
km  and  3.08  kg/ha/year  at  31  km). 

Annual  precipitation  increased  from  278  cm  at  4  km  to  351  cm  at  13  km 
from  the  coast  and  remained  constant  to  41  km  inland  (Fig.  26). 

Annual  and  Seasonal  Patterns  of  Deposition 

Precipitation  was  highest  in  the  late  autumn  and  winter  months  and  lowest 
during  the  summer  (Fig.  27).  The  rain  season  in  the  Hoh  River  valley  begins  in 
early  November  and  continues  through  April.  Deposition  rates  of  many  chemical 
constituents,  particularly  NO3  and  H,  had  seasonal  patterns  similar  to  that  of 
precipitation  (Fig.  27).  Seasonal  patterns  of  deposition  rates  for  K  and  excess 
Ca  were  not  as  apparent. 

Discussion 

Volume-weighted  mean  pH  of  bulk  precipitation  was  5.3  at  4  km,  5.5  at  13 
km,  5.6  at  24  km,  and  5.5  at  31  km  from  the  ocean.  Precipitation  pH  reported 
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4  13  24  31  41 

Distance  from  Pacific  Ocean  (l(m) 

Figure  26.  Mean  annual  precipitation  at  each  of  the  collection  sites  and  at  the  Hoh 
Ranger  Station  (National  Atmospheric  Deposition  Program/National  Trends  Net- 
work 1992),  1989-1992. 

in  this  study  (5.3  to  5.6)  was  within  the  range  reported  for  other  sites  along  the 
Pacific  Northwest  coast  and  at  remote  sites  elsewhere  in  the  world  (Vong  1989). 
Farther  inland,  at  sites  experiencing  anthropogenic  influences  on  atmospheric 
chemistry,  lower  precipitation  pH  has  been  reported.  Vong  et  al.  (1985)  re- 
ported volume-weighted  pH  of  4.41  in  the  metropolitan  Seattle  area  (approxi- 
mately 140  km  east  of  the  Hoh  River  valley). 

We  found  that  precipitation  chemistry  in  the  Hoh  River  valley  was  domi- 
nated by  sea  salt  ions,  Na  and  CI,  and  was  consistent  with  that  previously 
reported  for  the  Hoh  River  valley  by  Edmonds  et  al.  (1991).  Mean  electrical 
conductivity  of  precipitation  collected  at  each  site  ranged  from  0.47  to  1 .02 
mS/m  (Table  12).  The  correlations  among  the  major  sea  salt  ions  (CI,  SO4, 
Mg,  and  Na),  and  between  these  ions  and  electrical  conductivity,  are  indica- 
tive of  the  important  role  of  oceans  on  the  composition  and  ionic  strength  of 
bulk  precipitation  in  coastal  areas.  With  the  exception  of  Mg  and  Ca,  the  rela- 
tive contribution  of  sea  salt  ions  to  precipitation  chemistry  (Fig.  28)  remained 
constant  along  the  length  of  the  transect  even  though  concentrations  of  most 
constituents  decreased  as  distance  from  the  ocean  increased  (Table  12).  Much 
of  the  change  in  precipitation  chemistry  occurred  between  the  sampling  sites 
at  4  km  and  at  13  km.  These  changes  in  precipitation  chemistry  may  have  been 
dilution  effects  resulting  from  increased  precipitation  rates  farther  inland  along 
the  transect.  The  greatest  increase  in  rainfall  also  occurred  between  the  sites  at 
4  km  and  13  km  from  the  ocean  (Fig.  26). 

McColl  et  al.  (1982)  demonstrated  a  decrease  in  the  concentration  of  CI 
with  increasing  distance  inland  from  the  Pacific  coast  in  northern  California. 
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Figure  27.  Seasonal  deposition  rates  of  chemical  constituents  of  rainfall  (kg/ha)  and 
precipitation  depth  (mm)  at  each  end  of  transect,  1989-1992.  Bars  represent  deposi- 
tion during  each  4-week  collection  period  averaged  across  similar  time  periods  during 
1989-1992.  Error  bars  are  standard  errors.  Excess  Ca  and  Excess  SO^  are  Ca  and 
SO4  in  excess  of  those  expected  from  sea  salt. 


National  Park  Service  Monograph  2 


63 


In  that  study,  CI  concentration  decreased  rapidly  between  0  and  80  km  inland 
and  decreased  more  slowly  farther  inland.  The  results  of  our  study  also  indi- 
cate that  most  ion  concentrations,  including  CI,  decreased  as  distance  from  the 
coast  increased.  However,  the  ion  concentrations  tended  to  stabilize  closer  to 
the  coast  (13  to  24  km  inland)  than  was  reported  for  CI  (80  km  inland)  by 
McColl  et  al.  (1982).  Also,  McColl  (1982)  predicted,  using  CI  ratios,  that  K 
and  Ca  would  stabilize  20  to  30  km  inland  and  Na  and  Mg  at  80  to  100  km 
inland  at  the  northern  California  site. 

After  sea  salt  ions  (CI,  SO4,  Mg,  Na),  the  main  chemical  constituents  of 
precipitation  were  excess  Ca  and  excess  SO^.  Excess  Ca  concentrations  found 
in  bulk  precipitation  in  our  study  were  similar  to  those  reported  at  other  Pa- 
cific Northwest  sites  (Vong  1989).  The  presence  of  excess  Ca  is  indicative  of 
terrestrial  influences  on  precipitation  chemistry.  The  source  of  excess  Ca  in 
precipitation  was  most  likely  wind-blown  dust  formed  from  crustal  material. 
More  of  this  material  would  be  collected  in  bulk  deposition  collectors  used  in 
our  study  than  in  the  wet  deposition  collectors  used  in  other  studies,  because 
of  the  inclusion  of  dryfall  in  bulk  deposition  collectors.  The  high  concentra- 
tion and  deposition  rates  of  excess  Ca  at  some  of  the  sites  in  this  study  was 
most  likely  due  to  their  proximity  to  gravel  roads.  This  effect  may  have  been 
greatest  at  the  4  km  site  because  this  site  was  within  300  m  of  a  gravel  road. 
The  presence  of  Mg  in  local  bedrock  (Larson  1979)  suggests  that  some  of  the 
Mg  deposited,  like  Ca,  may  have  been  from  terrestrial  sources.  This  may  ac- 
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Figure  28.  Relative  contribution  of  ions  to  total  sea  salt  cations  in  precipitation  at 
four  distances  along  a  transect  from  the  Pacific  Ocean  to  Olympic  National  Park, 
February  1989-September  1992. 
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count  for  the  much  greater  concentration  and  annual  deposition  of  Mg  at  4  km 
than  at  3 1  km  and,  as  previously  noted,  the  higher  relative  abundance  of  Mg  com- 
pared to  other  sea  salt  ions  in  precipitation  collected  at  the  4  km  site  (Fig.  28). 

Using  the  Junge  equation  (Junge  1963)  and  known,  fine  SO^  aerosol  con- 
centrations (Schwartz  1988),  Vong  (1989)  determined  that  background  SO^ 
concentrations  in  precipitation  should  be  in  the  range  of  3  to  16  |ieq/L.  Mean 
SO4  concentrations  in  our  study  were  within  this  range  (Fig.  29).  However, 
excess  SO^  concentrations  were  lower  than  those  reported  for  other  sites  in  the 
Pacific  Northwest  (Vong  1989).  In  remote  coastal  areas,  the  source  of  this 
excess  SO4  is  the  marine  emission  of  (CH^)  2S  (dimethylsulfide  or  DMS)  which 
is  excreted  by  certain  marine  phytoplankton  (Andreae  and  Raemdonck  1983). 
DMS  reacts  in  the  atmosphere  to  form  SO4  and  methane  sulfonate. 

It  is  important  to  note  that  excess  SO4  concentration  was  the  only  measured 
parameter  that  correlated  with  precipitation  depth.  There  are  several  possible 
explanations  for  this  negative  correlation.  It  is  possible  that  increased  concen- 
trations of  SO4  in  the  atmosphere  coincided  with  the  dry  season.  This  may 
reflect  the  seasonal  production  and  emission  of  DMS  from  the  ocean.  Bates  et 
al.  (1987)  reported  higher  concentrations  of  DMS  in  temperate  oceans  during 
the  summer.  Also,  Bates  et  al.  (1987)  calculated  greater  flux  rates  of  DMS 
from  ocean  to  atmosphere  during  the  summer  than  in  the  winter  in  temperate 
oceans.  This  seasonal  change  corresponds  to  lower  rainfall  amounts  during  the 
summer  months  in  our  study  (Fig.  27). 


ONDJ    FMAMJ    JASONDJ    FMAMJ    JAS 
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Time 


Figure  29.  Average  monthly  SO^  concentrations  in  precipitation  near  West  Twin  Creek 
from  October  1986-September  1988. 
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The  correlation  between  excess  SO^  and  rainfall  depth  may  be  related  to  the 
role  of  non-sea-salt  SO^  as  cloud  condensation  nuclei  (CCN).  This  non-sea- 
salt  SO4  aerosol  has  been  reported  to  be  the  principal  nuclei  for  cloud 
condensation  in  the  marine  atmosphere  (Charlson  et  al.  1987).  Wameck  (1988) 
suggested  that  in  clouds,  essentially  all  nuclei  having  a  radius  greater  than  0.2 
|am  would  become  condensation  nuclei.  If  all  available  nuclei  have  been  en- 
gaged and  assuming  CCN  density  to  be  constant,  then  the  rate  of  condensation 
of  water  vapor  to  the  nuclei  would  control  the  amount  of  H2O  condensed  and 
subsequently  deposited.  The  correlation  between  excess  SO4  and  rainfall  depth 
may  simply  reflect  greater  rates  of  water  vapor  condensation  on  nuclei,  result- 
ing in  greater  rainfall.  This  would  result  in  decreased  excess  SO^  concentrations 
associated  with  increased  rainfall  depth.  During  the  summer  this  effect  would 
be  amplified  by  the  greater  CCN  density  in  the  atmosphere,  as  suggested  by 
the  calculated  flux  rates  of  DMS  from  ocean  to  atmosphere  reported  by  Bates 
et  al.  (1987).  If  the  suggestion  that,  in  clouds,  all  CCN  have  undergone  nucle- 
ation  (Wameck  1988)  is  correct,  then  there  would  be  a  potential  for  greater 
SO4  concentration  in  precipitation  during  the  summer,  even  with  the  reduced 
potential  for  high  rates  of  water  vapor  condensation. 

Another  explanation  for  the  negative  correlation  is  the  possibihty  that  a  large 
part  of  that  SO^  was  deposited  independently  of  precipitation  in  the  form  of 
dryfall.  Dry  deposition  of  SO^  was  reported  to  be  most  important  at  drier  sites 
(Lindberg  1992),  and  it  may  be  inferred  to  have  greater  impact  during  the  dry 
season  in  this  forest. 

Sulfate  deposition  (12.2  and  12.6  kg/ha/year)  was  similar  to  that  reported 
by  Larson  (1979)  in  the  Hoh  River  valley  for  the  water  year  1976-77  (11  kg/ 
ha/year).  The  deposition  rates  reported  by  Larson  (1979)  and  in  our  study  are 
greater  than  those  measured  in  wet  deposition  between  1988  and  1991  at  the 
Hoh  Ranger  Station  (8.7  kg/ha/year;  National  Atmospheric  Deposition  Pro- 
gram/National Trends  Network  1992).  This  difference  may  represent  the 
addition  of  SO^  as  dry  deposition  in  the  bulk  deposition  collectors  used  in  this 
study  and  by  Larson  (1979),  but  excluded  by  the  wet  deposition  collectors 
used  at  the  Hoh  Ranger  Station. 

The  SO4  deposition  rates  determined  in  this  study  were  also  similar  to  those 
reported  for  Washington  forests  in  the  Integrated  Forest  Study  (IFS),  but  were 
lower  than  SO^  deposition  in  IFS  forests  in  the  eastern  United  States  (Lindberg 
1992).  Vong  et  al.  (1988a)  suggested  that  on  the  Olympic  Peninsula,  precipita- 
tion chemistry  showed  little  influence  from  anthropogenic  sources,  especially 
during  periods  of  on-shore  air  flow. 

Total  N  deposition  in  bulk  precipitation  found  in  our  study  (3.11  to  3.42  kg/ 
ha/year)  was  less  than  that  reported  for  Washington  forests  in  IFS.  Lovett  (1992) 
reported  total  N  deposition  of  4.76  kg/ha/year  to  a  forest  approximately  180 
km  east  of  the  Hoh  River  valley.  That  IFS  study  used  wet  and  dry  deposition 
estimates  to  estimate  total  N  deposition.  The  dry  deposition  estimate  in  that 
study  included  deposition  of  coarse  and  fine  particulates  and  deposition  of 
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HNO3  vapor.  Because  total  N  deposition  in  our  study  was  based  on  bulk  depo- 
sition alone,  some  dry  deposition  components  would  have  been  included  by 
sedimentation.  However,  no  effort  was  made  in  our  study  to  actively  deter- 
mine total  particulate  impact  or  vapor  input. 
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Chapter  7:  Solution  Chemistry 
in  the  West  Twin  Creek  Watershed 

There  has  been  a  great  deal  of  interest  in  nutrient  retention  mechanisms  in 
forest  ecosystems  (Likens  et  al.  1970;  Vitousek  and  Reiners  1975;  Likens  et 
al.  1977;  Vitousek  etal.  1979;  Sollinsetal.  1980;  Johnson  and  Lindberg  1992; 
Van  Miegroet  et  al.  1992).  The  study  of  element  inputs  and  losses  and  internal 
transformations  in  solutions  as  they  flow  through  small  watersheds  is  impor- 
tant in  understanding  nutrient  retention  mechanisms.  Forest  canopies  strongly 
modify  precipitation  chemistry  (Tarrant  et  al.  1968;  Feller  1977;  Larson  1979; 
Sollinsetal.  1980;  Binkley  etal.  1982;  Cronan  and  Reiners  1983;  Parker  1983; 
Lindberg  et  al.  1986;  Edmonds  et  al.  1991 ;  Lindberg  1992;  Lovett  1992).  After 
leaving  the  canopy,  water  flowing  through  the  forest  floor  and  mineral  soil  is 
further  chemically  modified  so  that  by  the  time  it  reaches  a  stream,  its  ionic 
composition  may  be  vastly  different  from  that  of  precipitation  (Feller  1977; 
Likens  et  al.  1977;  Binkley  et  al.  1982;  Martin  and  Harr  1988;  Sollins  et  al. 
1980). 

Vitousek  and  Reiners  (1975)  related  nutrient  losses,  especially  for  N,  to 
ecosystem  succession.  In  undisturbed  ecosystems,  they  hypothesized  that  nu- 
trient inputs  equal  outputs.  Immediately  after  a  disturbance,  losses  of  limiting 
elements,  such  as  N,  exceed  inputs.  During  later  succession,  N  losses  are  low 
because  biomass  is  accumulating.  As  the  system  approaches  steady  state,  bio- 
mass  accumulation  ceases  and  inputs  equal  outputs  again.  Vitousek  and  Reiners 
(1975)  noted  that  mature  undisturbed  watersheds  in  New  Hampshire  had  rela- 
tively high,  constant  stream  NO3  outputs,  but  that  successional  ecosystems 
had  lower  stream  NO3  concentrations  with  pronounced  seasonal  cycles  in- 
volving growing  season  minima.  From  this  evidence,  they  concluded  that  plant 
uptake  regulated  N  losses  from  these  ecosystems.  They  also  suggested  that  K, 
Mg,  and  Ca,  which  are  essential  but  not  limiting  plant  nutrients,  were  more 
effectively  retained  in  successional  than  in  mature  ecosystems.  Mobile  ele- 
ments such  as  Na  and  CI  were  not  retained. 

Unlike  undisturbed,  mature  ecosystems  in  New  Hampshire  (Likens  et  al. 
1970;  Vitousek  and  Reiners  1975),  pristine,  old-growth  Douglas  fir  ecosys- 
tems in  Oregon,  had  low  stream  NO3  concentrations  (Sollins  et  al.  1980;  Martin 
and  Harr  1988).  Total  biomass  was  still  accumulating  despite  the  advanced 
age  of  these  ecosystems  (Sollins  et  al.  1980).  Old-growth  forests  in  water- 
sheds in  the  Pacific  Northwest  may  not  be  in  steady  state  because  disturbances 
caused  by  fire,  wind,  insects,  and  diseases  are  common.  Disturbance  patch 
size  may  also  influence  nutrient  retention  in  these  ecosystems  (Vitousek  and 
Reiners  1975). 

Van  Miegroet  et  al.  (1992)  recently  further  examined  NO3  retention  mecha- 
nisms in  forest  ecosystems  and  determined  that  the  degree  of  retention  depended 
on  relative  N  source  and  sink  strengths.  Sink  strengths  are  determined  both  by 
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natural  ecosystem  characteristics  (e.g.,  age  and  type  of  vegetation,  soil  age 
and  organic  matter  accumulation),  and  anthropogenic  influences  (e.g.,  man- 
agement practices  and  pollution).  Ecosystems  residing  on  more  recently 
developed  soils  (e.g.,  after  volcanic  ash  deposition  or  glaciation)  are  still  in  an 
active  stage  of  C  and  N  accumulation.  Nitrogen  retention  should  be  at  or  near 
maximum  in  these  ecosystems.  On  the  other  hand,  ecosystems  that  have  accu- 
mulated large  amounts  of  C  and  N,  because  of  age,  site  disturbance,  or  N 
fixation,  usually  have  large  NO3  leaching  losses.  Cole  et  al.  (1992)  suggested 
that  sites  with  a  narrow  soil  C:N  ratio  (<20:1)  and  with  a  total  soil  N  content 
>5000  kg/ha  appear  to  be  susceptible  to  nitrification  and  NO3  leaching. 

The  chemistry  of  solutions  flowing  through  the  West  Twin  Creek  watershed 
is  strongly  related  to  ecosystem  function.  Thus,  the  objectives  of  this  study 
were  to  determine  concentrations  of  major  cations  and  anions  and  dissolved 
organic  C  in  bulk  precipitation,  throughfall,  stemflow,  and  soil  solutions,  and 
nutrient  retention  mechanisms  in  this  pristine  watershed. 

Methods 

Sample  Collection 

Precipitation.  Precipitation  in  the  Hoh  River  valley  was  collected  bi-weekly 
from  October  1984  through  September  1993  in  bulk  precipitation  collectors. 
Precipitation  was  collected  in  an  autoclavable  polyethylene  bag.  From  1984  to 
1987,  four  collectors,  two  at  each  of  two  sites,  were  located  at  elevations  from 
229  to  640  m  in  clearcuts  among  a  16  km^  area  on  Washington  State  Depart- 
ment of  Natural  Resources  land  adjacent  to  Olympic  National  Park  and  as 
close  as  2  km  west  of  the  watershed  boundary  (Fig.  1).  Since  there  was  little 
variability  with  elevation  among  samplers,  two  collectors  at  only  one  site  were 
used  from  1988  to  1993.  Each  collector  was  a  20-cm  diameter  polyethylene 
funnel  with  a  netted  top  to  exclude  debris.  If  precipitation  fell  as  snow,  it  was 
collected  in  10-L  snow  buckets  lined  with  acid-rinsed  linear  polyethylene  (LPE) 
bags.  Samples  were  transported  in  triple-rinsed  1-L  polyethylene  bottles  in  a 
cooler  to  the  Analytical  Laboratory  at  the  College  of  Forest  Resources,  Uni- 
versity of  Washington,  within  48  hours  of  collecfion  and  stored  at  4°C. 

Throughfall,  Stemflow,  and  Soil  Solution.  Plots  for  studying  throughfall, 
stemflow,  and  soil  solution  analyses  were  established  in  the  West  Twin  Creek 
watershed  in  areas  representing  the  two  major  vegetation  zones  of  the  basin:  at 
the  base  (180  m  elevation)  and  top  (800  m  elevation)  of  the  watershed.  Sam- 
pling began  October  1986  and  was  conducted  through  December  1993  in  the 
plot  at  the  base  of  the  watershed,  but  only  through  September  1988  in  the  plot 
at  the  top  of  the  watershed.  Samples  were  collected  biweekly. 

Throughfall  and  stemflow  were  collected  around  four  tree  species.  The  trees 
were  of  different  diameters  and  ages;  at  lower  elevation,  two  western  redcedar 
(40  and  1 12  cm  dbh),  two  Douglas  fir  (145  and  247  cm  dbh),  and  two  western 
hemlock  (49  and  77  cm  dbh)  trees,  and  at  higher  elevation,  two  western  hem- 
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lock  (73  and  81  cm  dbh)  and  two  Pacific  silver  fir  (59  and  141  cm  dbh)  trees 
were  monitored.  The  trees  were  in  a  2.2-ha  area  in  the  upper  watershed  and  a 
3.3-ha  area  in  the  lower  watershed. 

Throughfall  collectors  were  constructed  of  10-cm  PVC  rain  gutters  extend- 
ing from  the  bole  of  a  sampling  tree  to  the  edge  of  the  crown.  One  collector  per 
tree  was  used;  they  were  placed  such  that  throughfall  would  be  collected  from 
that  tree  only.  Collectors  were  lined  with  acid-washed  and  triple-rinsed  linear 
polyethylene  and  covered  with  2-mm  nylon  chiffon  mesh  to  keep  out  litterfall. 
Stemflow  collars  were  made  from  2.5-cm  diameter  Tygon  tubing  that  was 
halved  lengthwise  and  fastened  to  the  lower  boles  of  the  trees,  encircling  the 
stem.  Throughfall  and  stemflow  were  collected  in  120-L  polyethylene  con- 
tainers lined  with  polyethylene  bags.  Volumes  were  determined  in  the  field 
and  a  1-L  sample  from  each  collector  was  returned  to  the  laboratory  for  chemi- 
cal analysis. 

Soil  solutions  were  collected  using  tension  lysimetry  (Cole  1968).  Sixty- 
millimeter  diameter,  fritted  glass  Buchner  funnels  were  modified  by  cutting 
off  the  funnel  1  cm  above  the  fritted-glass  surface  and  bending  the  stem  at 
right  angles.  The  modified  funnels  were  placed  beneath  the  forest  floor  at  5  cm 
and  40  cm  below  the  top  of  mineral  soil.  Tension  of  0.01  MPa  was  applied 
using  a  1-m  hanging  water  column.  In  the  lower  watershed,  at  the  same  loca- 
tions as  the  throughfall  and  stemflow  collectors,  three  replicate  samples  for 
each  soil  horizon  were  taken.  No  lysimeter  samples  were  collected  at  the  up- 
per watershed  plot.  Soil  solutions  were  collected  in  10-L  polyethylene  bottles, 
and  volumes  were  recorded  in  the  field.  If  individual  sample  volumes  were 
insufficient  for  analysis,  samples  were  composited  by  horizon.  Soil  solutions 
were  sampled  monthly,  except  during  dry  periods,  and  they  were  handled  in 
the  same  manner  as  other  solutions.  Also,  in  the  same  locations  as  the  throughfall 
and  stemflow  collectors  in  upper  and  lower  watershed  plots,  soil  samples  of 
the  O  horizon,  0  to  5  cm  depth  (BA  horizon)  and  6  to  40  cm  depth  (Bw  hori- 
zon) of  the  mineral  soil  were  obtained.  Samples  were  placed  in  a  cooler  and 
transported  to  the  Analytical  Laboratory  at  the  College  of  Forest  Resources, 
University  of  Washington,  where  they  were  air  dried. 

Solution  Chemical  Analysis 

Specific  conductance,  pH  and  alkalinity  were  determined  on  non-filtered 
samples  within  a  few  days  after  arrival  at  the  laboratory.  Specific  conductance 
was  determined  with  a  YSI  Model  31  Conductivity  Bridge  (Yellow  Springs 
Instrument  Co.,  Yellow  Springs,  OH)  and  pH  with  a  Radiometer  PHM  85  pH 
meter  (Radiometer,  Copenhagen,  Denmark).  Alkalinity  was  determined  by  ti- 
tration (Rhoades  1982)  to  an  end  point  of  pH  5. 

After  these  measurements,  the  water  samples  were  filtered  through  a 
Whatman  GF/A  filter,  stored  at  4°C  and  analyzed  within  a  month  after  arrival 
at  the  laboratory.  Solutions  were  analyzed  for  Ca,  Mg,  K,  and  Na  using  an  IL 
95 1  Atomic  Absorption  Analyzer  (Instrumentation  Laboratory,  Wilmington, 
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MA).  Ammonium  and  NO3  were  determined  with  a  Technicon  Autoanalyzer 
II  (Technicon,  Tarrytown,  NY).  Sulfate,  CI,  and  PO^  were  determined  with  a 
Dionex  2100  Ion  Chromatograph  (Dionex,  Sunnyvale,  CA).  Dissolved  organic 
carbon  (DOC)  was  determined  using  persulfate  oxidation  at  100  °C  and  an  OI 
Model  700  Organic  Carbon  Analyzer  (OI  Corporation,  College  Station,  TX). 
Precipitation,  stemflow,  throughfall,  and  soil  solution  quantities  were  used  to 
determine  volume-weighted  average  concentrations. 

Soil  Chemistry  Analysis 

Available  N  (NH^  and  NO3)  was  determined  on  moist  field  samples  imme- 
diately upon  arrival  at  the  laboratory  using  1  M  KCl  extraction  (Keeney  and 
Nelson  1982)  and  a  Technicon  Autoanalyzer  II  (Technicon,  Tarrytown,  NY). 
Samples  were  then  air-dried  for  other  analyses;  moisture  correction  was  ap- 
plied by  drying  a  subsample  at  105  °C.  Total  C  was  determined  using  a  Leco  C 
analyzer  (Leco  Corp.,  St.  Joseph,  MI).  Total  N  was  determined  using  an  H2SOy 
LiS04/H202/Se  digestion  (Parkinson  and  Allen  1975)  and  a  Technicon  Atomic 
Absorption  Analyzer  II  (Technicon,  Tarrytown,  NY).  Phosphorus  (Bray  No.  1) 
was  determined  using  the  method  of  (Jackson  1958)  and  a  Technicon 
Autoanalyzer  II  (Technicon,  Tarrytown,  NY).  Sodium,  K,  Ca,  Mg,  and  Al  were 
determined  with  1  M  NH^Cl  extraction  and  an  IL  95 1  Atomic  Absorption  Ana- 
lyzer (Instrumentation  Laboratory,  Wilmington,  MA).  Cation  exchange  capacity 
was  determined  by  extraction  with  1  M  NH^Cl  and  a  saturation  wash  with  1  M 
KCl  (Oak  Ridge  National  Laboratory,  Oak  Ridge,  TN,  EPRI  Forest  Response 
Program  Protocol).  A  1:1  soil:H20  paste  was  used  for  determination  of  pH 
with  a  Radiometer  PHM  85  pH  meter. 

Quality  Control  (QC) 

Sampling  and  analysis  of  water  samples  from  precipitation,  throughfall, 
stemflow,  and  the  stream  followed  QA/QC  procedures  outlined  by  Environ- 
mental Protection  Agency  (EPA).  Field  samples  were  collected  in  acid-washed 
and  triple-rinsed  containers,  labeled,  placed  in  coolers  as  soon  as  possible,  and 
transported  to  the  Analytical  Laboratory  at  the  University  of  Washington  within 
48  hours  of  collection  and  stored  at  4  °C.  Samples  were  generally  analyzed  for 
all  constituents  within  4  weeks  of  collection.  Standard  laboratory  protocol  in- 
cluded replications  and  10%  National  Bureau  of  Standards  rain-water  standards 
and  blanks  during  chemical  analysis. 

Results  and  Discussion 

Precipitation  Chemistry 

Precipitation  pH  averaged  5.3  (Table  14),  which  is  typical  for  relatively  pris- 
tine sites  in  the  western  United  States  (Larson  1979;  Martin  and  Harr  1988; 
Bormann  et  al.  1989).  Slightly  lower  pH  (4.5-4.8)  have  been  reported  in  sites 
near  Vancouver,  British  Columbia  (Feller  1977;  Binkley  et  al.  1982)  and  Se- 
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Table  14.  Average  pH,  electrical  conductivity,  and  dissolved  organic  carbon  (DOC) 
of  precipitation,  stemflow  and  throughfall  of  Douglas  fir,  western  hemlock,  western 
redcedar,  and  Pacific  silver  fir,  and  soil  solutions  in  the  West  Twin  Creek  watershed 
October  1986-September  1988.  Numbers  in  parentheses  are  standard  deviations. 


pH 

Electrical  conductivity^ 

DOC^ 

mS/m 

mg/L 

Precipitation 

5.3a 

1.15a 

1.5a 

(0.3) 

(0.63) 

(1.1) 

Throughfall 

Lower  watershed 

Douglas-fir 

5.1a 

2.24ab 

10.5a 

(0.2) 

(0.13) 

(17.9) 

western  hemlock 

5.1a 

1.82ab 

7.5ab 

(0.2) 

(0.06) 

(5.9) 

western  redcedar 

5.1a 

2.1  lab 

7.  lab 

(0.2) 

(1.97) 

(5.1) 

Upper  watershed 

Pacific  silver  fir 

4.9ab 

2.76ab 

15.6b 

(0.2) 

(1.28) 

(6.6) 

western  hemlock 

5.1a 

1.67a 

8.5b 

(0.2) 

(0.79) 

(4.7) 

Stemflow 

Lower  watershed 

Douglas-fir 

4.  led 

7.02c 

25.5c 

(0.4) 

(2.68) 

(9.3) 

western  hemlock 

4.5bc 

4.03bc 

26.5c 

(0.4) 

(1.71) 

(20.4) 

western  redcedar 

4.0d 

5.25bc 

26.0c 

(0.2) 

(1.78) 

(12.2) 

Upper  watershed 

Pacific  silver  fir 

4.5bc 

9.80c 

33.8c 

(0.4) 

(5.20) 

(25.9) 

western  hemlock 

4.4cd 

4.56bc 

26.1c 

(0.3) 

(2.42) 

(9.9) 

Soil  solution  (lower  watershed) 

Forest  floor 

5.7ae 

1.78a 

10.8b 

(0.1) 

(0.82) 

(1.1) 

A 

5.8e 

1.53a 

9.0b 

(0.1) 

(0.67) 

(4.6) 

B 

6.2f 

2.6ab 

2.9a 

(0.1) 

(0.65) 

(1.1) 

^  Numbers  followed  by  a  different  letter  are  significantly  different  (/?  =  0.5). 
^'DOC  determined  only  in  1987-1988. 
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attle,  Washington  (Van  Miegroet  1986),  probably  because  of  anthropogenic 
inputs.  Electrical  conductivity  (EC)  of  precipitation  was  low  (1.15  mS/m;  Table 
14),  which  suggests  low  concentrations  of  ions  in  the  precipitation. 

Concentrations  of  total  cations  and  total  inorganic  anions  in  precipitation 
were  relatively  even,  with  only  a  small  anion  deficit  (Table  15).  Cations  minus 
anions  yielded  a  deficit  of  8.4  |amol^/L  (Table  15),  and  the  cation:anion  ratio 
was  1.10,  both  indicating  that  organic  anions  were  fewer.  This  is  supported  by 
the  low  concentration  of  dissolved  organic  carbon  (DOC)  in  precipitation  (1.5 
mg/L;  Table  14). 

Sodium,  at  an  average  concentration  of  54.3  joeq/L  (Table  15),  was  the  ma- 
jor cation  in  precipitation  followed  by  Ca,  Mg,  and  K.  Reflecting  the  proximity 
of  the  site  to  the  ocean,  the  high  Na  concentration  is  2  to  5  times  higher  than 
those  reported  from  inland  sites  (Feller  1977;  Sollins  et  al.  1980;  Wolfe  1988). 
High  Na  concentrations  in  precipitation  have  been  found  in  other  studies  near 
salt  water  bodies  (Attiwill  1966;  Carlisle  et  al.  1966;  Larson  1979;  Bormann  et 
al.  1989).  At  sites  in  the  Appalachian  and  Cascade  Ranges,  located  more  than 
50  km  from  an  ocean,  cation  concentrations  rank  Ca  >  Na  >  K  >  Mg  (Henderson 
et  al.  1978).  Likewise,  chloride  was  the  dominant  anion  in  precipitation  (57.6 
\AQqfL)  because  of  proximity  to  the  ocean.  Chloride  concentrations  were  1.5  to 
4  times  higher  than  those  reported  from  sites  farther  inland  (Feller  1977;  Sollins 
et  al.  1980;  Wolfe  1988)  but  were  similar  to  those  at  other  coastal  sites  (Larson 
1979;  Bormann  et  al.  1989).  Anion  concentrations  ranked  CI  >  HCO3  >  S04> 
NO3  >  PO4  (Table  15). 

Sulfate  concentration  averaged  11.2  joeq/L,  of  which  47%  of  the  SO4  was 
non-sea-salt  SO4  (NSS;  Edmonds  et  al.  1991).  NSS  SO4  is  derived  from  oxida- 
tion of  dimethyl  sulfide  produced  by  marine  algae  (Charlson  et  al.  1987),  dust, 
and  wood  smoke  from  slash  burning.  The  SO^  concentration  in  precipitation 
was  lower  than  that  in  precipitation  in  Puget  Sound  (Wolfe  1988)  and  in  south- 
western British  Columbia  (Feller  1977)  where  anthropogenic  inputs  could  be 
expected,  and  much  lower  than  those  observed  in  the  eastern  United  States 
(Vitousek  and  Reiners  1975;  Likens  et  al.  1977;  Johnson  and  Lindberg  1992). 
However,  SO^  concentration  in  precipitation  was  similar  to  that  observed  in 
the  Oregon  Cascades,  (Sollins  et  al.  1980)  and  at  coastal  sites  (Larson  1979; 
Vong  and  Larson  1983;  Bormann  et  al.  1989).  From  1986  to  1988,  highest 
concentrations  of  SO4  were  observed  in  spring,  late  summer,  and  fall  (Fig.  29). 

In  the  Hoh  River  valley,  concentrations  of  NH^  in  precipitation  were  low 
(0.4  )Lieq/L;  Table  15).  Higher  NH^  concentrations  in  precipitation  have  been 
found  at  study  sites  in  and  near  regions  of  human  development  and  anthropo- 
genic sources  in  Puget  Sound  (Wolfe  1988)  and  southwestern  British  Columbia 
(Feller  1977).  Nitrate  concentration  in  precipitation  (1.4  |ieq/L;  Table  15)  was 
higher  than  that  of  NH^.  Similar  to  the  pattern  of  NH^  found  in  precipitation, 
these  NO3  concentrations  are  comparable  to  those  observed  in  precipitation  in 
the  Oregon  Cascades  (Sollins  et  al.  1980),  coastal  sites  (Bormann  et  al.  1989), 
lower  than  NO3  concentrations  in  precipitation  in  Puget  Sound  (Wolfe  1988), 
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and  much  lower  than  those  observed  in  the  eastern  United  States  (Vitousek 
and  Reiners  1975;  Likens  et  al.  1977;  Johnson  and  Lindberg  1992).  There 
were  no  obvious  seasonal  patterns  in  NH^  and  NO3  concentrations. 

Throughfall  Chemistry 

Most  coniferous  canopies  show  a  tendency  toward  net  acidification  of  bulk 
precipitation  inputs  (Cronan  and  Reiners  1983).  Acidification  processes  may 
involve  uptake  of  NH^,  nitrification,  washing  off  of  dry  deposition,  and  leach- 
ing of  plant-derived  acids.  In  conifers,  the  acidification  processes  tend  to 
dominate  over  neutralization  processes.  Neutralization  processes  may  involve 
plant  surface-ion  exchange,  and  release  of  organic  or  bicarbonate  salts  (Cronan 
and  Reiners  1983). 

Throughfall  pH  in  West  Twin  Creek  averaged  near  5.0  (Table  14),  only 
slightly  lower  than  the  incoming  precipitation  pH  of  5.3.  While  a  similar  trend 
was  found  in  southwestern  British  Columbia  (Feller  1977)  and  coastal  Wash- 
ington (Larson  1979),  the  opposite  trend  was  found  in  Cedar  River  watershed 
in  the  Puget  Sound  region,  Washington  (Cole  and  Johnson  1977;  Van  Miegroet 
1986),  and  Cascade  Head,  Oregon  (Tarrant  et  al.  1968).  The  increase  in  pH  at 
Cedar  River  and  Cascade  Head  from  precipitation  to  throughfall  may  be  due 
to  appreciable  quantities  of  HCO3  in  throughfall  (Parker  1983).  Concentra- 
tions of  HCO3  in  throughfall  in  the  Cedar  River  watershed  were  considerably 
higher  (>40  |ieq/L;  Cole  and  Johnson  1977)  than  that  in  the  throughfall  (5.4 
|Lieq/L)  and  precipitation  (8.0  |ieq/L)  in  West  Twin  Creek  watershed  (Table  15). 

Although  there  was  no  significant  difference  between  precipitation  and 
throughfall  pH,  there  was  a  trend  towards  acidification  in  throughfall.  The 
dominant  acidification  mechanism  in  the  watershed  appeared  to  be  leaching  of 
plant-derived  acids.  Overall,  dissolved  organic  carbon  (DOC)  in  throughfall 
increased  (Table  14),  indicating  organic-acid  leaching.  Throughfall  DOC  con- 
centrations were  significantly  different  from  precipitation  DOC  only  under 
Pacific  silver  fir  and  western  hemlock  in  the  upper  watershed  (Table  14).  Pa- 
cific silver  fir  throughfall  pH  was  consistently  lower  than  that  of  western 
hemlock  (Fig.  30).  The  high  DOC  in  Pacific  silver  fir  throughfall  (Table  14) 
may  contribute  to  low  pH  values,  in  particular,  because  Pacific  silver  fir  retain 
needles  longer  than  other  species  (Kimmins  1987),  and  the  needles  may  de- 
compose on  the  branch,  contributing  organic  acids  to  throughfall.  There  were 
no  significant  differences  in  throughfall  pH  or  DOC  among  species  (Table  14). 
McDowell  and  Likens  (1988),  Ugolini  et  al.  (1988),  and  Koprivnjak  and  Moore 
(1992)  also  reported  increases  in  DOC  concentration  in  throughfall  relative  to 
that  in  the  precipitation. 

Throughfall  pH  displayed  slight  seasonal  patterns.  Generally,  higher  pH 
values  were  observed  in  spring  and  summer  and  the  lower  pH  values  were 
observed  in  fall  and  winter  (Fig.  30).  The  seasonal  pattern  may  be  related  to 
decomposition  activities  in  the  canopy  and  the  lower  rainfall  volume  in  spring 
and  summer. 
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Figure  30.  Average  monthly  pH  of  throughfall  beneath  (A)  western  hemlock,  western 
red  cedar,  and  Douglas-fir  in  the  lower  West  Twin  Creek  watershed,  and  (B)  western 
hemlock  and  Pacific  silver  fir  in  the  upper  West  Twin  Creek  watershed,October 
1986-July  1988. 


There  was  a  trend  for  both  electrical  conductivity  (EC,  Table  14)  and  cation 
concentrations  (Table  15)  to  be  higher  in  throughfall  than  in  precipitation,  but 
not  significantly.  As  for  precipitation,  the  major  cation  in  throughfall  was  Na 
(76.5  |aeq/L;  Table  15).  Throughfall  cation  concentrations  ranked  Na  >  K  > 
Mg  >  Ca.  Throughfall  was  enriched  relative  to  precipitation  in  the  order  K  > 
Na  >  Mg  >  Ca,  indicating  that  K  was  most  readily  leached  from  the  canopy. 
Cronan  and  Reiners  (1983)  noted  that  most  forest  canopies  show  a  marked 
enrichment  in  ions  as  rainfall  leaches  through  the  canopy. 
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Leaching,  washing-off  of  deposited  particles,  and  evaporation  from  tree 
crowns  all  contribute  to  cation  enrichment  in  throughfall.  Leaching  is  prob- 
ably the  dominant  mechanism  at  West  Twin  Creek  watershed  for  both 
throughfall  and  stemflow,  which  is  more  enriched  in  cations  than  throughfall 
(Table  15).  Washing  off  of  dry  deposition  or  marine  aerosols  containing  sea 
salt  appears  important,  based  on  the  increased  concentrations  of  Na  and  CI  in 
throughfall  and  stemflow  relative  to  precipitation  (Table  15).  Evaporation  is 
probably  less  important  since  SO^  concentration  changed  little  with  passage 
through  the  canopy  (Table  15).  Sulfate  commonly  shows  little  interaction  with 
foliage  (Lindberg  et  al.  1986). 

Often  through  foliar  uptake,  NH^  concentrations  in  precipitation  decrease 
after  interaction  with  the  canopy  (Feller  1977;  Parker  1983).  However,  NH^ 
concentrations  in  throughfall  were  not  lower  than  that  in  precipitation  at  West 
Twin  Creek  (Table  15).  In  fact,  there  was  a  slight,  but  not  significant,  increase. 

Dominant  anions  in  throughfall  were  CI  (74.7  |aeq/L)  and  SO^  (1 1.6  )Lieq/L; 
Table  15).  The  mean  concentrations  of  all  anions  were  not  significantly  differ- 
ent from  those  in  precipitation.  Sulfate  concentration  in  throughfall  was  lower 
than  that  in  throughfall  in  southwestern  British  Columbia  (Feller  1977)  and 
Cedar  River  watershed,  Washington  (Johnson  1975),  but  it  was  similar  to  that  in 
throughfall  in  the  Oregon  Cascades  (Sollins  et  al.  1980).  Again,  higher  SO^  con- 
centrations in  the  Puget  Sound  area  near  Seattle  are  due  to  anthropogenic  inputs. 

There  were  no  significant  differences  in  the  concentrations  of  any  of  the 
cations  or  anions  in  throughfall  among  species  (Table  16).  The  highest  Ca 
concentration  in  throughfall  occurred  under  western  redcedar,  a  known  Ca  ac- 
cumulator (Bledsoe  and  Zasoski  1981).  Pacific  silver  fir  tended  to  influence 
throughfall  chemistry  to  a  greater  extent  than  the  other  species,  producing 
throughfall  with  the  highest  mean  concentrations  of  K,  Mg,  Na,  CI,  SO^,  and 
PO4  (Table  16). 

In  summary,  the  chemistry  of  throughfall  at  West  Twin  Creek  watershed 
was  similar  to  that  of  precipitation  and  was  dominated  by  oceanic  influences. 
There  were  trends  toward  increases  in  acidity  and  enrichment  of  DOC  and 
cations,  probably  due  mostly  to  leaching,  but  in  no  case  were  concentrations 
of  cations  in  throughfall  significantly  greater  than  those  in  precipitation. 

Stemflow  Chemistry 

Stemflow  solutions  typically  had  higher  ionic  concentrations  than  precipi- 
tation and  throughfall.  This  is  usual  for  most  forest  ecosystems  (Parker  1983). 
There  was  also  a  much  larger  anion  deficit  in  stemflow  (225.5  jLimol^/L)  than 
in  throughfall  (56. 1  |amol^/L;  Table  15).  The  stemflow  cation/anion  ratio  ranged 
from  1.91  to  2.60,  and  the  throughfall  cation/anion  ratio  ranged  from  1.30  to 
1.42.  This  was  largely  due  to  the  presence  of  organic  anions  in  stemflow,  as 
indicated  by  the  high  DOC  concentrations  (Table  14).  Koprivnjak  and  Moore 
(1992)  noted  a  similar  trend  in  throughfall  and  stemflow  DOC  in  a  lichen  wood- 
land in  Quebec. 
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The  concentrations  of  cations  and  anions,  except  HCO3  and  NO3,  tended  to 
be  higher  in  stemflow  solution  than  in  throughfall  and  precipitation  solution, 
but  only  in  the  case  of  K  was  stemflow  concentration  significantly  higher  (Table 
15).  Potassium  ( 1 19.6  |ieq/L)  and  Na  ( 1 1 8.0  |neq/L)  were  the  dominant  cations 
in  stemflow,  and  the  dominant  anion  was  CI  (160.2  |Lieq/L;  Table  15).  The 
ranking  of  cation  concentrations  was  K>Na>Mg>Ca.  Similar  to  throughfall, 
stemflow  was  enriched  relative  to  precipitation  in  the  order  of  K>Ca>Mg>Na, 
indicating  that  K  was  easily  leached  from  the  canopy.  The  dominant  mechanism 
for  cation  enrichment  for  both  throughfall  and  stemflow  is  probably  leaching. 

Commonly,  NH^  and  NO^  in  precipitation  are  taken  up  by  the  canopy,  re- 
sulting in  negative  net  throughfall  and  stemflow  (Parker  1983).  This  was 
observed  of  NH^  in  throughfall  in  southwestern  British  Columbia  (Feller  1977) 
and  of  NO3  in  Oregon  (Sollins  et  al.  1980)  and  Washington  (Johnson  1975).  In 
the  West  Twin  Creek  watershed,  NH^  concentrations  did  not  differ  among 
stemflow  (2.7  |Lieq/L),  precipitation  (0.4  |aeq/L),  and  throughfall  (0.6  |LXeq/L; 
Table  15).  However,  the  trend  toward  a  higher  NH^  concentration  in  both 
throughfall  and  stemflow  suggests  that  N  was  fixated  in  the  old-growth  canopy. 
Denison  (1973)  noted  that  the  epiphytic  lichen  Lobaria  oregana  fixes  N  in  the 
canopy  of  old-growth  trees  in  the  Oregon  Cascades,  and  this  species  occurs  in 
the  West  Twin  Creek  watershed.  Nadkami  (1984)  also  found  that  in  the  Olym- 
pic rainforest  the  amounts  of  N  were  higher  in  throughfall  beneath  branches 
with  epiphytes  than  beneath  stripped  branches. 

The  pH  was  lower  and  DOC  and  EC  were  higher  in  stemflow  than  in 
throughfall  (Table  14).  Stemflow  pH  ranged  from  4.0  to  4.5  (average  of  4.3) 
and  was  significantly  lower  than  precipitation  and  throughfall  pH.  Plant-de- 
rived acids  strongly  controlled  stemflow  acidity.  Stemflow  pH  also  depended 
on  species  (Table  14).  Stemflow  pH  was  lowest  in  western  redcedar  (4.0)  and 
Douglas  fir  (4. 1 )  and  highest  in  western  hemlock  (4.5).  Stemflow  pH  was  similar 
in  western  hemlock  in  both  the  upper  and  lower  watershed  (Table  14).  The 
bark  is  rougher  of  western  redcedar  and  Douglas  fir  than  of  western  hemlock, 
providing  a  greater  surface  area  for  contact  and  chemical  reaction.  Western 
redcedar  bark  also  contains  high  concentrations  of  phenolics  (Barton  1963). 
The  pH  of  stemflow  solution  from  Pacific  silver  fir  was  relatively  higher  (4.5) 
despite  the  comparably  high  DOC  concentration  (Table  14).  Seemingly,  the 
high  cation  concentrations,  especially  Ca,  in  Pacific  silver  fir  stemflow  (Table 
16)  had  a  neutralizing  influence  on  the  solution. 

Tree  species  also  influenced  cation  and  anion  concentrations  in  stemflow. 
The  concentrations  of  cations  and  anions  were  generally  highest  in  stemflow 
solution  of  Pacific  silver  fir  and  relatively  low  in  that  of  western  redcedar 
(Table  16).  Phosphate  concentrations  were  generally  high  in  stemflow  of  Pa- 
cific silver  fir  and  Douglas  fir  (Table  16),  suggesting  that  P  is  being  leached 
from  the  canopies  of  these  species. 

Cronan  and  Reiners  (1983)  suggested  that  cation  enrichment  in  solutions 
flowing  through  forest  canopies  is  mostly  from  plant  recycling,  whereas  most 
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of  the  inorganic  anions  are  from  atmospheric  sources.  In  the  West  Twin  Creek 
watershed,  the  same  mechanisms  are  operating,  although  wash-off  of  sea-salt- 
derived  aerosols  is  probably  responsible  for  much  of  the  Na  enrichment.  The 
anions  CI,  SO4,  and  NO3  were  dominantly  from  atmospheric  sources,  whereas 
PO4  enrichment  was  more  likely  from  plant  sources. 

In  summary,  cations  tended  to  be  more  concentrated  in  stemflow  than  in 
throughfall  and  precipitation.  The  pH  of  stemflow  also  was  much  lower,  which 
was  caused  by  organic  acids  from  plant  sources.  Organic  anions  from  plant 
sources  dominated  inorganic  anions  in  stemflow  and  were  associated  with  cat- 
ion leaching  in  the  canopy.  Although  only  5%-20%  of  the  volume  of 
precipitation  reaches  the  forest  floor  in  stemflow  (Parker  1983),  it  is  a  major 
source  of  plant  nutrients  to  fine  roots  near  the  base  of  trees  and  could  be  an 
important  source  of  nutrients  to  plant  roots  in  the  West  Twin  Creek  watershed. 

Soil  Solution  Chemistry 

Whereas  the  pH  of  solutions  moving  through  the  West  Twin  Creek  forest 
canopy  decreased  from  5.3  in  precipitation,  5.0  in  throughfall,  and  4.3  in 
stemflow,  it  increased  in  the  forest  floor  solution  to  5.7  (Table  14).  Soil  solu- 
tion pH  increased  further  at  15 -cm  depth  to  5.8,  and  at  40-cm  depth  to  6.2 
(Table  14).  Soil  pH  also  increased  with  depth  from  4.2  in  the  forest  floor  to  4.9 
in  the  lower  mineral  soil  (Table  17).  The  pH  is  commonly  lower  in  soil  than  in 
soil  solution.  Adsorption  of  organic  anions  onto  soil  particles  or  weathering  of 
silicate  minerals  and  increase  in  HCO3  may  be  responsible  for  the  increasing 
pH  in  solution  with  soil  depth.  Concentrations  of  DOC  and  anion  deficit  de- 
creased with  soil  depth  (Tables  14  and  15).  The  cation: anion  ratio  decreased 
from  1.50  in  the  forest  floor  solution  to  1.05  at  the  40-cm  soil  depth,  indicating 
that  adsorption  of  organic  anions  is  an  important  mechanism  for  neutralizing 
the  natural  acidity  in  stemflow  and  throughfall  in  this  ecosystem. 

The  chemical  composition  of  a  soil  solution  depends  on  a  complex  series  of 
equilibrium,  adsorption,  immobilization,  plant  uptake,  displacement,  weather- 
ing, and  decomposition  reactions  (Feller  1977).  The  leaching  of  cations  depends 
on  the  amount  of  anions  in  the  solution.  Organic  anions  and  HCO3  are  the 
dominant  anions  associated  with  leaching  in  undisturbed  Pacific  Northwest 
forests,  and  organic  anions  dominate  at  low  pH  (McColl  1972;  Johnson  et  al. 
1977;  Ugolini  et  al.  1977).  Organic  anions  dominated  the  canopy  leaching,  but 
as  the  solution  moved  through  the  forest  floor,  organic  anions  became  less 
important  and  HCO3  dominated  at  the  40-cm  depth  (Table  15). 

Concentrations  of  cations  were  generally  lower  in  forest  floor  leachates  than 
in  stemflow  yet  slightly  higher  than  in  throughfall  (Table  15).  Binkley  et  al. 
(1982)  also  found  that  cations  were  higher  in  forest  floor  solutions  than  in 
throughfall  in  coastal  forests  in  British  Columbia.  Cation  concentrations  tended 
to  be  lower  in  the  mineral  soil  than  in  the  forest  floor  (Table  15),  perhaps 
reflecting  plant  uptake.  The  dominant  cation  in  the  forest  floor  and  mineral 
soil  solutions  was  again  Na  (Table  15).  Cation  concentrations  generally  de- 
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creased  in  the  order  of  Na  >  Ca  >  Mg  >  K  in  the  forest  floor  and  in  mineral  soil. 
Anion  concentrations  decreased  in  the  order  of  CI  >  HCO3  >  SO^  >  PO^  > 
NO3.  Because  chloride  moves  with  Na  in  solution,  it  had  little  influence  on  the 
leaching  of  nutrient  cations. 

The  dominant  anion  in  the  soil  solution  away  from  coastal  influence  is  HCO3, 
particularly  in  low-elevation  forests  (Feller  1977;  Johnson  et  al.  1977;  Sollins 
et  al.  1980).  HCO3  concentration  was  relatively  higher  in  soil  solution  than  in 
throughfall  and  stemflow  solutions,  probably  because  of  decomposition  and 
root  respiration.  The  difference  of  SO^  concentrations  was  little  among  pre- 
cipitation, throughfall,  stemflow,  and  soil  solutions  (Table  15),  but 
concentrations  tended  to  increase  slightly  with  depth  in  the  mineral  soil,  per- 
haps as  a  result  of  weathering.  Phosphate  tended  to  be  lower  in  the  soil  solution 
than  in  stemflow  and  throughfall  (Table  15)  because  of  plant  uptake  and  soil 
retention  mechanisms. 

Nitrate  concentrations  in  the  soil  solution  and  bulk  soil  NO3  were  low  and 
changed  little  with  depth  (Tables  15  and  17).  In  the  forest  floor  and  soil  solu- 
tions, NO3  concentrations  were  generally  lower  (0.1-0.3  |ieq/L)  than  NH^ 
(0.6-3.1  |ieq/L;  Table  15).  A  trend  in  the  bulk  soil  was  similar  (Table  17). 
Sollins  et  al.  (1980)  and  Feller  (1977)  found  similar  ranges  of  NO3  in  the  soil 
solution.  Feller  (1977)  also  found  similar  soil  NH^  concentrations  (2.77  ^teq/ 
L)  in  forest  floor  and  soil  solutions  in  southwestern  British  Columbia.  Thus,  N 
was  retained  in  the  ecosystem  despite  the  very  high  rainfall.  However,  increased 
atmospheric  NO3  inputs  in  the  1993-96  period  were  not  retained  as  discussed 
in  Chapter  8,  and  it  seems  that  the  ecosystem  is  on  the  verge  of  N  saturation. 
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Chapter  8:  Long-Term  Trends  in 

Precipitation  and  Stream  Chemistry 

in  the  West  Twin  Creek  Watershed 

The  influence  of  anthropogenic  sources  of  H,  SO4,  NH^,  and  NO3  on  for- 
ested watersheds  has  been  the  subject  of  much  research  (Swank  and  Waide 
1988;  Driscoll  et  al.  1988;  Federer  et  al.  1989;  Johnson  and  Van  Hook  1989; 
Aber  et  al.  1991;  Johnson  and  Lindberg  1992;  Hedin  et  al.  1994;  Likens  and 
Bormann  1995).  Human  induced  long-term  changes  in  precipitation  chemis- 
try have  occurred  at  many  sites  in  eastern  North  America  and  Europe  (Driscoll 
et  al.  1989;  Likens  et  al.  1990;  Hedin  et  al.  1994),  and  similar  changes  in 
stream  chemistry  have  also  been  noted  at  Hubbard  Brook  in  New  Hampshire 
(Driscoll  et  al.  1989),  in  the  Catskill  Mountains  of  New  York  (Stoddard  1991), 
and  in  eastern  Canada  (Clair  and  Ehrman  1995). 

Strong  seasonal  patterns  in  precipitation  and  stream  chemistry  have  also 
been  noted  in  eastern  North  America  (Swank  and  Waide  1988;  Likens  and 
Bormann  1995),  but  few  long-term  monitoring  studies  of  precipitation  and 
stream  chemistry  have  been  conducted  in  coastal  or  near-coastal  sites  in  west- 
em  North  America,  which  are  considered  to  be  relatively  pollution  free.  Most 
studies  have  been  short  term  (Larson  1979;  Sollins  et  al.  1980;  Martin  and 
Harr  1988;  Bormann  et  al.  1989;  Edmonds  et  al.  1995),  although  precipitation 
and  stream  chemistry  were  monitored  from  1971  to  1986  in  a  small  forested 
basin  in  southwestern  British  Columbia  (Feller  1987).  This  site,  however,  was 
influenced  by  mildly  acid  precipitation.  Galloway  et  al.  (1995)  predicted  that 
atmospheric  N  on  a  global  scale  will  continue  to  rise  because  of  increased 
fossil-fuel  combustion  and  fertilizer  use;  about  two-thirds  of  the  increase  are 
in  Asia.  Nitrogen  deposition  to  forests  in  western  North  America  could  in- 
crease if  there  is  long-range  transport  across  the  Pacific  Ocean.  Thus, 
examination  of  long-term  trends  in  atmospheric  inputs  and  outputs  of  N  and 
other  elements  in  relatively  pristine  areas  to  detect  long-range  transport  and 
global  change  is  particularly  important. 

Our  objectives  were  the  determination  of  trends  in  precipitation  and  stream 
chemistry  over  a  9-year  period  in  a  pristine,  old-growth  forest  watershed  and 
the  determination  of  seasonal  patterns  in  precipitation  and  stream  chemistry. 
This  chapter  is  focused  on  data  collected  from  1984  to  1993,  but  trends  in  NO3 
in  precipitation  and  the  stream  between  1993  and  1997  are  also  presented. 

Methods 

Sample  Collection  and  Analytical  Procedures 

Samples  were  collected  biweekly  from  November  1984  through  June  1990 
and  at  4-week  intervals  thereafter  through  December  1997.  Bulk  precipitation 
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was  collected  in  a  clear-cut  immediately  adjacent  to  the  Olympic  National 
Park  boundary  and  1-2  km  from  the  West  Twin  Creek  watershed  (Fig.  1).  Two 
bulk  precipitation  collectors  were  located  in  the  sampling  site.  Each  bulk  pre- 
cipitation collector  consisted  of  a  20-cm  diameter  polyethylene  funnel  with  a 
netted  top  to  exclude  debris.  The  precipitation  drained  from  the  funnel  into  a 
polyethylene  bag.  Samples  were  collected  in  triple-rinsed  1-L  polyethylene 
bottles  and  within  48  hours  of  collection  transported  in  a  cooler  to  the  Analyti- 
cal Laboratory  at  the  College  of  Forest  Resources,  University  of  Washington, 
where  they  were  stored  at  4°C.  Collection  volumes  were  determined.  A  Na- 
tional Acid  Deposition  Program  (NADP)  station  is  also  located  at  the  Hoh 
Ranger  Station  (Fig.  1;  National  Acid  Deposition  Program-National  Trends 
Network  1994). 

Stream  water  samples  were  collected  with  a  Qualimetric  stage  proportional 
sampler  just  upstream  from  a  V-notch  weir  at  the  base  of  the  West  Twin  Creek 
watershed.  Samples  were  returned  to  the  laboratory  for  chemical  analysis  in 
the  same  manner  as  precipitation  samples. 

Electrical  conductivity  (EC),  pH,  and  alkalinity  (reported  as  HCO3)  were 
determined  on  nonfiltered  samples.  EC  was  determined  with  a  YSI  Model  31 
Conductivity  Bridge  (Yellow  Springs  Instrument  Co.,  Yellow  Springs,  Ohio) 
and  corrected  to  25°C.  The  pH  was  determined  with  a  Radiometer  PHM85  pH 
meter  (Radiometer,  Copenhagen,  Denmark).  Alkalinity  was  determined  by  ti- 
tration to  an  end  point  of  pH  5  (Homann  et  al.  1992).  Remaining  sample  material 
was  filtered  through  Whatman  GF/A  filters  and  stored  at  4°C  and  analyzed 
within  1  month  after  arrival  at  the  College  of  Forest  Resources  laboratory  at 
the  University  of  Washington. 

Solutions  were  analyzed  for  Ca,  Mg,  K,  and  Na  with  a  IL  95 1  Atomic  Ab- 
sorption Analyzer  (Instrumentation  Laboratory,  Wilmington,  Massachusetts) 
prior  to  July  1993  and  by  Inductively  Coupled  Plasma  atomic  emission  spec- 
trometry thereafter.  Samples  were  analyzed  by  both  methods  for  6  months  to 
ensure  that  a  change  in  technique  did  not  influence  the  data.  Ammonium  was 
determined  with  a  Technicon  Autoanalyzer  II  (Technicon,  Tarrytown,  New 
York).  Sulfate,  CI,  NO3  and  PO4  were  determined  with  a  Dionex  2100  Ion 
Chromatograph  (AS4A  and  AG4A  columns;  Dionex,  Sunnyvale,  California). 
Macro-  Kjeldahl  digests  were  done  on  all  samples  of  sufficient  volume,  and 
those  digests  were  analyzed  for  N  and  P  (reported  here  as  Total  Kjeldahl  N  and 
P  (TKN  and  TKP)). 

QA-QC  Procedures 

Standard  chemical  analysis  QA-QC  laboratory  procedures  were  followed 
including  use  of  replications.  National  Bureau  of  Standards  rain  water  stan- 
dards and  blanks.  Data  were  transferred  from  laboratory  notebooks  to  computer 
spreadsheets  and  edited.  Outliers  were  checked  and  data  were  eliminated  if 
field  notes  at  the  time  of  sampling  indicated  obvious  problems  with  the  sample's 
quality.  Precipitation  samples  were  also  screened  on  the  basis  of  Cl:Na  ratio. 
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Samples  with  Cl:Na  ratios  that  deviated  from  the  sea  sah  ratio  of  1.17  by  a 
factor  of  2  were  excluded.  Only  a  few  samples  were  excluded,  however. 

Statistical  Analyses 

Time  trends  and  seasonal  differences  were  examined  with  Autoregressive 
Integrated  Moving  Average  Analysis  (ARIMA;  SPSS  1994).  Data  were  log 
transformed.  Because  of  the  irregular  sampling  schedule,  missing  values  were 
filled  with  the  mean  of  two  nearby  points. 

Results  and  Discussion 

Precipitation  Chemistry 

There  was  no  significant  change  in  EC,  pH,  or  concentrations  of  TKN  and 
TKP  in  precipitation  from  1984  to  1993  (Fig.  31).  EC  was  low,  reflecting  the 
low  concentrations  of  ions  in  precipitation  in  this  site,  and  averaged  0.77  mS/m 
(Table  18).  The  pH  averaged  5.3  and  ranged  from  4.3  to  7.1  (Fig.  31),  which  is 
typical  for  relatively  pristine  sites  in  the  western  United  States  (Larson  1979; 
Martin  and  Harr  1988;  Bormann  et  al.  1989),  and  it  was  similar  to  that  deter- 
mined at  the  nearby  NADP  site  at  the  Hoh  Ranger  Station  (5.4;  National  Acid 
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Figure  31.  Trends  in  EC,  pH,  and  concentrations  of  TKN  and  TKP  in  precipitation 
falling  near  West  Twin  Creek,  1984-1993. 
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Deposition  Program-National  Trends  Network  1994).  In  general,  the  chemis- 
try of  precipitation  in  the  NADP  site  and  in  our  collection  site  was  similar 
(Edmonds  et  al.  1995).  Slighdy  lower  pH  (4.5-4.8)  has  been  reported  in  sites 
near  Vancouver,  British  Columbia  (Feller  1977;  Binkley  et  al.  1982;  Feller 
1987)  and  Seattle,  Washington  (Van  Miegroet  1986),  probably  because  of  an- 
thropogenic inputs.  Volume-weighted  annual  average  concentrations  of  TKN 
and  TKP  were  relatively  low,  averaging  0. 16  mg/L  of  TKN  and  0.015  mg/L  of 
TKP  (Table  18).  Total  N  ranged  from  0.33  mg/L  to  0.56  near  Vancouver,  Brit- 
ish Columbia  (Feller  1987). 

Time  trends  in  concentrations  of  anions  in  precipitation  were  absent  (Fig. 
32),  except  in  NO3  concentrations,  which  showed  a  slight  declining  trend  (p 
<0.001).  Between  1984  and  1993,  there  was  little  evidence  yet  of  transport  of 
NO3  or  NH4  across  the  Pacific  Ocean  from  Asia,  and  concentrations  varied 
little  from  year  to  year  (Table  18;  Edmonds  and  Blew  1997).  Continued  data 
collection,  however,  revealed  that  NO3  concentrations  have  increased  since 
1993,  suggesting  that  the  trans-Pacific  transport  of  NO3  predicted  by  Gallo- 
way et  al.  (1995)  may  now  be  occurring.  This  increased  NO3  was  also  reflected 
in  NADP  sites  in  coastal  Oregon  and  California.  Chloride  was  the  dominant 
anion  in  precipitation  because  of  the  proximity  to  the  ocean,  (volume- weighted 
annual  average  45.1  |aeq/L;  Table  18).  Sulfate  concentrations  averaged  8.61 
|ieq/L,  HCO3  6.42  )aeq/L,  NO3  0.70  |ieq/L,  and  PO^  0.07  |aeq/L  (Table  18). 
Non-sea-salt  SO4  can  account  for  as  much  as  47%  of  the  SO^  (Edmonds  et  al. 
1995)  and  is  derived  from  the  oxidation  of  dimethyl  sulfide  produced  by  ma- 
rine algae,  dust,  and  wood  smoke  from  slash  burning.  Martin  and  Harr  (1988) 
reported  HCO3  as  the  dominant  anion  in  precipitation  in  the  Cascade  Moun- 
tains of  Oregon,  160  km  from  the  Pacific  Ocean.  In  sites  with  urban  influences 
on  precipitation  chemistry  in  the  Northwest,  SO^  is  the  dominant  anion  (Binkley 
et  al.  1982;  Vong  et  al.  1985;  Feller  1987).  Increases  in  SO^  and  NO3  were 
observed  from  1971  to  1986  near  Vancouver,  British  Columbia,  reflecting  the 
urban  influence  (Feller  1987). 

Average  annual  SO^  concentrations  in  precipitation  are  higher  in  the  eastern 
United  States  than  we  observed  (i.e.,  typically  >40  )Lieq/L).  However,  at  Hubbard 
Brook  in  New  Hampshire,  annual  volume-weighted  concentrations  of  SO^  in 
precipitation  decreased  from  more  than  60  |ieq/L  in  1964  to  near  40  |Lieq/L  in 
1987  with  a  slight  increase  in  pH  from  about  4.1  to  4.3  (Driscoll  et  al.  1989). 
This  was  associated  with  a  decline  in  SO2  emissions  in  the  eastern  United 
States.  Average  annual  NO3  concentrations  in  precipitation  are  also  much  higher 
in  the  eastern  United  States  (15-35  |ieq/L)  than  we  reported.  Emissions  of 
NOx  in  the  eastern  United  States  influence  NO3  concentrations  in  precipita- 
tion and,  unlike  SO2,  NOx  emissions  have  not  declined  (Likens  et  al.  1990). 

Concentrations  of  anions  in  precipitation  were  without  strong  seasonal  pat- 
terns (Fig.  32),  but  concentrations  of  SO^  were  typically  highest  in  late  spring 
and  summer.  This  may  have  been  due  to  the  influence  of  non-sea-salt  SO^  as 
described  above.  Bates  et  al.  (1987)  calculated  greater  flux  rates  of  DMS  from 
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Figure  32.  Trends  in  concentrations  of  anions  in  precipitation  falling  near  the  West 
Twin  Creek,  1984-1993. 


the  ocean  to  the  atmosphere  during  summer  than  winter,  which  may  result  in 
greater  concentrations  of  SO^  in  the  atmosphere  at  that  time.  Blew  and  Edmonds 
(1995)  reported  SO^  concentration  to  be  negatively  correlated  with  precipita- 
tion amount.  In  addition  to  the  seasonality  of  DMS  emission  from  the  ocean, 
Blew  and  Edmonds  (1995)  suggested  that  this  may  be  the  result  of  increased 
dry  deposition  of  SO^  during  the  summer  dry  season.  This  is  an  important 
consideration  because  these  samples  were  collected  as  bulk  precipitation  rather 
than  only  as  wet  precipitation. 
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The  dominant  cation  in  precipitation  was  Na  followed  by  Ca  >  Mg  >  K  > 
NH4  (Fig.  33;  Table  18).  The  high  Na  concentration  (2-5  times  higher  than 
more  inland  sites  (Feller  1977;  SoUins  et  al.  1980))  reflects  the  proximity  to 
the  ocean.  Average  annual  volume-weighted  concentrations  of  Na,  Ca,  Mg,  K, 
and  NH4  were  42.4,  6.58,  5.44,  1.81  and  0.12  |aeq/L  (Table  18).  Hydrogen 
averaged  4.546  |ieq/L.  Other  studies  near  salt  water  also  report  high  Na  con- 
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Figure  33.  Trends  in  concentrations  of  cations  in  precipitation  falling  near  the  West 
Twin  Creek,  1984-1993. 
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centrations  (Bormann  et  al.  1989).  In  sites  located  more  than  50  km  from  the 
ocean,  the  Ca  concentration  is  greater  than  Na  concentration  (Sollins  et  al. 
1980).  Concentrations  of  NH^  in  precipitation  were  less  in  the  Hoh  River  Val- 
ley than  in  the  Puget  Sound  (Wolfe  1988)  and  British  Columbia  (Feller  1987), 
reflecting  anthropogenic  influence.  In  precipitation  in  our  site,  nitrate  concen- 
trations were  higher  than  NH^  concentrations,  which  is  typical  in  sites  near  the 
ocean  (Likens  and  Bormann  1995). 

We  found  no  seasonal  trends  in  concentrations  of  cations  in  precipitation 
(Fig.  33),  but  the  sum  of  cations  (Ca,  K,  Mg,  Na;  /7<0.001),  especially  of  Ca, 
significant  declined.  In  fact,  the  volume- weighted  annual  averages  of  Ca  steadily 
decUned  from  13.91  |aeq/L  in  1988  to  1.45  |aeq/L  in  1992  (Table  18).  Declines 
of  Mg  and  Na  were  similar.  This  was  probably  associated  with  a  decrease  in 
timber  harvest  in  the  area  after  1990.  The  bulk  precipitation  collectors  were 
located  in  a  clear-cut  in  a  managed  forest  immediately  adjacent  to  Olympic 
National  Park  ajid  were  within  a  few  hundred  meters  of  gravel  logging  roads. 
Dust  from  these  roads  had  probably  influenced  the  chemistry  of  bulk  precipi- 
tation during  periods  of  heavy  use  of  the  nearby  roads. 

Concentrations  of  base  cations  in  precipitation  (sum  of  non-sea-salt  Ca,  Mg, 
K,  and  Na),  especially  Ca,  also  declined  from  1963  to  1991  over  large  areas  in 
western  Europe  and  eastern  North  America  (Hedin  et  al.  1994).  This  decline  is 
attributed  to  reductions  in  point  source  emissions  from  fuel  combustion,  in- 
dustrial processes,  road  paving,  and  changes  in  agricultural  tilling  practices. 
Natural  emissions  from  forest  fires  and  wind  erosion  of  arid  soils  are  also 
important  sources  of  base  cations.  This  reduction  in  base  cations  could  offset 
the  reductions  of  SO^  deposition  in  many  areas  and  could  contribute  to  in- 
creased sensitivity  of  poorly  buffered  ecosystems. 

Stream  Chemistry 

There  were  no  changes  in  pH,  TKN,  and  TKP  in  stream  water  of  West  Twin 
Creek  through  time  or  season,  but  seasonal  differences  in  EC  were  significant 
(/7<0.001;  Fig.  34).  EC  was  highest  in  late  summer  before  the  rainy  season 
began.  Electrical  conductivity  was  1  order  of  magnitude  higher  in  the  stream 
than  in  precipitation,  reflecting  the  high  concentration  of  ions  from  weather- 
ing in  the  stream.  The  average  stream  pH  was  7.3,  and  the  range  was  5.5-9.0. 
Larson  (1979)  recorded  similar  pH  values  in  nearby  streams,  and  Martin  and 
Harr  (1988)  noted  stream  pH  values  above  7.0  in  Alaska,  British  Columbia, 
and  Oregon.  Stottlemyer  (1992)'observed  stream  pHs  as  high  as  8.5  in  Denali 
National  Park,  Alaska.  In  southwestern  British  Columbia,  stream  water  pH 
was  lower,  averaging  6.6  from  1971  to  1985.  However,  pH  values  of  less  than 
6. 1  were  not  recorded,  even  though  mean  annual  precipitation  pH  values  ranged 
from  4.5  to  4.9.  This  indicates  that  the  moderate  acidic  precipitation  in  this  site 
had  not  yet  significantly  affected  stream  water  chemistry  (Feller  1987).  The 
EC  was  1  order  of  magnitude  lower  in  this  British  Columbia  stream  than  in  the 
West  Twin  Creek  watershed.  Annual  average  volume-weighted  concentrations 
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of  TKN  and  TKP  in  West  Twin  Creek  were  relatively  low:  0.15  mg/L  of  TKN 
and  0.017  mg/L  of  TKP  (Table  19)  and  similar  to  those  in  precipitation. 
Edmonds  et  al.  (1995)  had  previously  reported  that  N  and  P  losses  from  West 
Twin  Creek  were  similar  to  inputs. 
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Figure  34.  Trends  in  EC,  pH,  and  concentrations  of  TKN  and  TKP  in  stream  water  of 
the  West  Twin  Creek,  1984-1993. 
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Concentrations  of  anions  in  the  stream  (Fig.  35)  except  SO^,  which  increased 
significantly  (/?<  0.001),  were  without  time  trends.  Bicarbonate  was  the  domi- 
nant anion  in  stream  water  (volume-weighted  annual  average  589.3  |ieq/L), 
followed  by  SO^  (323.6  |ieq/L),  CI  (  79.9  |aeq/L),  NO3  (6.06  |ieq/L),  and  PO4 
(0.33  |ieq/L;  Table  19).  Hydrogen  averaged  only  0.052  |ieq/L.  Larson  (1979) 
in  nearby  streams  and  Stottlemyer  (1992)  in  coastal  Alaska  found  similar  pat- 
terns of  stream  water  chemistry.  Concentrations  of  SO^  were  higher  than 
observed  in  many  streams  in  the  eastern  United  States  (Stoddard  1991)  and 
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Figure  35.  Trends  in  concentrations  of  anions  in  stream  water  of  the  West  Twin  Creek, 
1984-1993. 
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Canada  (Clair  and  Ehrman  1995).  This  strongly  reflects  the  parent  material  in 
our  site,  which  consists  of  uplifted  marine  sediments.  Stream  water  SO^  con- 
centrations are  now  decreasing  in  many  streams  in  the  eastern  United  States  as 
a  result  of  reduced  atmospheric  SO^  inputs  (Likens  and  Bormann  1995; 
Stoddard  1991).  At  Hubbard  Brook,  New  Hampshire,  however,  stream  pH  did 
not  increase  despite  the  decrease  in  SO^,  perhaps  because  of  a  decrease  in 
basic  cations.  Because  all  major  elements  act  together  to  regulate  the  acidity 
of  soils  and  drainage  waters,  deposition  of  all  elements,  including  cations  and 
NO3,  needs  to  be  considered.  A  combination  of  increasing  NO3  in  precipita- 
tion and  decreasing  Ca  and  Mg  seems  to  be  outweighing  the  influence  of 
declining  trends  in  SO4  inputs  and  outputs  in  the  eastern  United  States. 

Concentrations  of  NO3  in  West  Twin  Creek  did  not  increase  with  time  (Fig. 
35)  and  were  1  order  of  magnitude  lower  in  eastern  streams  where  NO3  seems 
to  be  increasing  (Johnson  and  Lindberg  1992;  Stoddard  1991  and  1994)  be- 
cause of  increased  atmospheric  N  deposition  and  possible  N  saturation  (Aber 
et  al.  1991).  However,  continued  monitoring  since  1993  revealed  an  increase 
in  stream  NO3  concentrations,  mirroring  the  increase  in  NO3  concentrations  in 
precipitation.  Stream  pH  has  also  been  depressed  since  1993. 

Seasonal  patterns  were  strong  in  NO3 ,  EC  (Fig.  34),  HCO3 ,  and  SO4  (Fig. 
35;  /?  <  0.001).  The  pattern  associated  with  EC,  HCO3,  and  SO^  in  stream 
water  seems  to  be  related  to  seasonal  patterns  of  stream  discharge.  As  stream 
discharge  decreased  during  summer  and  fall,  concentrations  of  the  dominant 
ions  increased.  Concentrations  of  the  dominant  ions  reached  their  maximum 
prior  to  the  onset  of  the  wet  season  and  fell  drastically  after  the  beginning  of 
the  rainy  season.  Presumably,  as  surface  flow  in  the  watershed  decreased,  a 
greater  proportion  of  the  water  in  the  stream  was  from  subsurface  or  ground- 
water discharge  into  the  stream.  This  water  has  had  greater  interaction  with  the 
parent  material  and  consequently  carried  greater  concentrations  of  the  prod- 
ucts of  weathering  reactions.  We  did  not  determine  weathering  rates  in  the 
West  Twin  Creek  watershed.  However,  stream  Si02  concentration  averaged 
39.9  |ieq/L  from  1991  to  1993  (R.  L.  Edmonds,  unpublished  data),  indicating 
active  weathering.  Concentrations  of  weathering  products  in  the  stream  water 
are  related  to  stream  discharge  rates.  These  results  were  consistent  with  those 
reported  by  Feller  (1977)  in  a  forested  watershed  in  southwestern  British  Co- 
lumbia. In  addition  to  the  weathering  of  carbonate  minerals,  an  alternate  source 
for  HCO3  is  the  dissolution  of  CO2  released  from  respiration  in  the  soil.  Sollins 
et  al.  (1980)  presumed  that  increases  in  HCO3  as  solutions  passed  through  the 
mineral  soil  were  solely  because  of  respiration.  However,  in  the  streams  stud- 
ied by  Sollins  et  al.  ( 1 980),  the  amount  of  HCO3  was  1  order  of  magnitude  less 
than  we  report  here  of  West  Twin  Creek. 

Seasonal  patterns  of  NO3  in  stream  water  were  a  little  different  from  those 
of  the  dominant  ions.  Nitrate  concentrations  increased  slightly  with  decreas- 
ing stream  flow.  But  the  more  prominent  pattern  was  a  large  increase  in  NO3 
concentrations  after  the  onset  of  the  rainy  season  and  a  subsequent  decline  to 
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background  levels  soon  after  (Fig.  35).  Plant  uptake  of  N  is  far  less  during  this 
period  in  late  fall  and  early  winter  than  in  spring  and  summer.  The  loss  of  this 
N  sink  may  have  resulted  in  an  accumulation  of  NO3  in  the  soil  that  would 
then  be  leached  from  the  surface  soil  by  the  first  significant  rain  events  of  the 
wet  season  and  be  subsequently  carried  to  the  stream.  Other  investigators  have 
reported  high  NO3  concentrations  in  stream  water  peaking  in  early  spring 
(Vitousek  and  Reiners  1975;  Stoddard  1994;  Likens  and  Bormann  1995)  after 
the  onset  of  snowmelt  or  spring  rainstorms. 

The  seasonality  of  NO3  concentrations  we  observed  in  West  Twin  Creek 
also  appears  in  the  Hoh,  Queets,  Quinault,  and  Quilayute  rivers  of  Washington 
(U.S.  Geological  Survey  1972,  1973,  1986),  suggesting  that  this  pulse  of  NO3 
in  streams  is  common  over  a  wide  area  on  the  western  Olympic  Peninsula. 
Regardless  of  its  source,  this  pulse  of  NO3,  represented  a  considerable  portion 
of  the  annual  total  N  loss  from  the  West  Twin  Creek  watershed.  Edmonds  et  al. 
(1995)  estimated  annual  NO3  discharge  from  this  stream  to  be  2.25  kg/ha/year 
and  total  N  loss  to  7. 1 3  kg/ha/year.  We  estimated  the  loss  of  N  from  the  water- 
shed as  NO3  during  this  pulse  to  be  approximately  1.8  kg/ha/year,  which 
represented  approximately  25%  of  the  total  annual  discharge  of  N  from  this 
watershed.  This  contrasts  to  reports  of  watershed  NO3  loss  in  other  sites  in  the 
Pacific  Northwest.  For  example,  Martin  and  Harr  (1988)  reported  NO3-N  losses 
from  watersheds  in  the  Cascade  Mountains  of  Oregon  of  only  0. 1  kg/ha/year. 

The  dominant  cation  in  stream  water  was  Ca  (Fig.  36).  Volume-weighted 
annual  concentrations  averaged  733.8  |ueq/L  of  Ca  followed  by  Na  (180.2  |ieq/L) 
>Mg  (154.0  |Lieq/L)  >K  (10.82  |ieq/L)  >NH4  (1.02  jneq/L;  Table  19).  The  high 
Ca  concentration  and  the  high  HCO3  and  SO^  concentrations  in  the  stream 
reflect  active  weathering.  High  Na  concentrations  reflect  the  proximity  to  the 
ocean  as  well  as  weathering.  Annual  deposition  of  Na  in  bulk  precipitation  to 
this  watershed  during  1986-88  (41  kg/ha/year)  was  equal  to  approximately 
one-third  of  that  discharged  from  the  watershed  by  the  stream  ( 1 24  kg/ha/year; 
Edmonds  etal.  1995). 

Seasonal  patterns  of  Ca,  Mg,  and  Na  (Fig.  36)  resembled  those  of  HCO3  and 
SO4  and  were  statistically  significant  (p  <  0.001).  Potassium,  NH^,  and  H  were 
without  seasonal  patterns.  In  the  stream  NH^  concentrations  were  lower  than 
NO3  concentrations.  Plant  uptake  and  microbial  and  soil  processes  are  impor- 
tant processes  for  NH^  retention  in  the  ecosystem.  Similar  general  patterns  of 
stream  water  chemistry  have  been  reported  of  other  watersheds  in  the  North- 
west. Binkley  et  al.  (1982)  in  coastal  British  Columbia  and  Martin  and  Harr 
(1988)  in  the  Oregon  Cascade  Mountains,  reported  Ca  and  HCO3  as  the  domi- 
nant ions  in  stream  water. 

Long-term  Trends 

The  concentration  of  most  analytes  in  precipitation  and  stream  water  showed 
few  statistically  significant  trends  from  1984  to  1993.  However,  NO3  in  pre- 
cipitation declined  slightly,  whereas  SO^  in  stream  water  increased.  There  also 
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Figure  36.  Trends  in  concentrations  of  cations  in  stream  water  of  the  West  Twin  Creek, 
1984-1993. 


appeared  to  be  a  trend  toward  decreasing  concentrations  of  cations  in  precipita- 
tion. Since  1993,  however,  NO3  concentrations  in  precipitation  and  in  the  stream 
have  increased  and  pH  has  decreased  (Fig.  37):  This  seems  to  be  across-Pacific 
transport  of  NO3  from  Asia.  The  sharp  increase  in  NO3  in  the  stream  during  the 
1993-1996  period  suggests  that  the  ecosystem  is  on  the  verge  of  N  saturation. 
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Figure  37.  Concentrations  of  NO3-N  in  (A)  precipitation  falling  near  the  West  Twin 
Creek  watershed  and  (B)  in  the  stream  in  the  watershed,  1984-1997.  Stream  pH  is 
also  shown  (C). 
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Appendix  A:  Soil  Profile  Descriptions. 

Soil  Profile  Description  of  the  Lower  Portion 
of  the  West  Twin  Creek  Watershed 

General  description  of  soil  pit  site 

Classification:  Typic  Dystrochrept 

Deep,  well-drained  soil  derived  from  colluvial  material  on  a  bench  between 
streams.  Slope  is  20%,  aspect  is  160°,  and  elevation  is  244  m  (800  ft).  Potential 
natural  vegetation  is  Picea  sitchensis,  Tsuga  heterophylla,  and  Vaccinium  parvi- 
folium  with  a  ground  cover  dominated  by  mosses.  Area  is  disturbed  by  windthrow. 

Horizons 

Oi     10.2  to  6.4  cm 

Oe     6.4  to  1.3  cm 

Oa     1.3  to  0  cm;  abrupt  wavy  boundary;  pH  3.86 

Bwl  0  to  33  cm,  with  a  thickness  of  12.7  to  40.6  cm;  moist  color  (lOYR  4/ 
4),  clay  loam  texture;  moderate,  medium  subangular  blocky  structure;  firm 
moist  consistence,  slightly  sticky,  slightly  plastic;  30%  gravel,  5%  cobble,  and 
5%  stone;  many  fine,  medium,  and  coarse  roots;  very  fine  and  fine  interstitial 
and  vesicular  pores;  clear  broken  boundary;  pH  4.55 

Bw2  33  to  61  cm,  with  a  thickness  of  20  to  30  cm;  moist  color  (lOYR  5/4); 
clay  loam  texture;  moderate,  medium  subangular  blocky  structure;  firm  moist 
consistence,  slightly  sticky,  slightly  plastic;  some  cutans  present;  25%  gravel, 
5%  cobble,  and  10%  stone;  few  fine,  medium,  and  coarse  roots;  very  fine  and 
fine  interstitial  and  vesicular  pores;  clear  broken  boundary;  pH  5.06;  irregular 
patches  of  Bh  present  (color  lOYR  2/1) 

Bw3  61  to  102  cm,  with  a  thickness  of  38  to  46  cm;  moist  color  (7.5 YR  4/6); 
loam  texture;  moderate,  coarse  subangular  blocky  structure;  firm  moist  consis- 
tence, slightly  sticky,  slightly  plastic;  25%  gravel,  5%  cobble,  and  10%  stone;  very 
few  fine  and  medium  roots;  very  fine  and  fine  intersdfial  and  vesicular  pores;  clear 
smooth  boundary;  pH  5.49;  irregular  patches  of  Bh  present  (color  lOYR  2/1) 

Bw4  102  cm  and  down;  moist  color  (lOYR  5/6);  silty  clay  loam  texture;  mod- 
erate, medium  subangular  blocky  structure;  friable  moist  consistence,  slightly 
sdcky,  slightly  plastic;  20%  gravel,  5%  cobble,  and  50%  stone;  very  few  me- 
dium roots;  pH  5.75 
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Soil  Profile  Description  of  the  Upper  Portion 
of  the  West  Twin  Creek  Watershed 

General  description  of  soil  pit  site 

Classification:  Typic  Dystrochrept 

Moderately  deep,  fairly  well-drained  soil  derived  from  marine  sediments  of 
sandstone  and  siltstone.  Slope  is  20%,  aspect  is  180°,  and  elevation  is  671  m 
(2200  ft).  Profile  is  located  on  a  slight  convexity  in  a  midslope  position.  Natu- 
ral vegetation  is  Tsuga  hetewphylla  and  Abies  amabilis,  with  a  depauperate 
understory. 

Horizons 

Oi     4  to  2  cm 

Oa     2  to  0  cm,  pH  3.80 

EB  0  to  5  cm,  with  a  thickness  of  2.5  to  10.2  cm;  moist  color  (7.5YR  3/3); 
silty  clay  loam  texture;  moderate,  medium  subangular  blocky  structure  parting 
to  crumb;  firm  moist  consistence,  slightly  sticky,  slightly  plastic;  5%  gravel 
with  no  cobble  or  stones;  many  micro,  common  very  fine,  fine,  and  medium, 
and  few  coarse  roots;  many  tubular  and  common  vesicular  pores;  clear  wavy 
boundary;  pH  3.49 

Bsl  5  to  28  cm,  with  a  thickness  of  15  to  28  cm;  moist  color  (7.5YR  4/5), 
clay  loam  texture;  moderate  to  strong,  coarse  subangular  blocky  structure;  very 
firm  moist  consistence,  slightly  sticky,  slightly  plastic;  5%  gravel;  few  micro, 
very  few  fine,  medium,  and  coarse  roots;  common  tubular  pores;  clear  wavy 
boundary;  pH  3.73 

Bs2  28  to  41  cm,  with  a  thickness  of  7.6  to  17.8  cm;  moist  color  (7.5 YR  4/6); 
clay  loam  texture;  moderate,  coarse  subangular  blocky;  firm  to  friable  moist 
consistence,  slightly  sticky,  slightly  plastic;  5%  gravel;  common  micro  and 
fine,  few  very  fine  and  medium  roots;  common  tubular  pores;  clear  wavy  bound- 
ary; pH  4.40 

Bs3  41  to  69  cm,  with  a  thickness  of  10  to  30  cm;  moist  color  (7.5YR  4/6); 
sandy  clay  loam  texture;  moderate,  medium  subangular  blocky  structure  part- 
ing to  crumb;  firm  moist  consistence,  slightly  sticky,  sHghdy  plastic;  20%  gravel; 
few  fine,  medium  and  coarse  roots;  few  tubular  pores;  gradual  wavy  bound- 
ary; pH  4.38 
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BC  69  cm  and  down;  too  rocky  to  sample  for  color,  structure,  and  texture; 
40%  gravel,  20%  cobble;  few  fine  and  coarse  roots,  few  tubular  and  vesicular 
pores;  pH  4.58 

Soil  Profile  Description  of  Soil  near  the  Hoh  Lake  Watershed 
General  description  of  soil  pit  site 

Classification:  Typic  Cryumbrept 

Vegetation  type:  Abies  amabili  si  Rhododendron  albiflorum  association 

Physiographic  position:  mountain  upland,  1249  m  (4100  fi)  elevation 

Topography:  moderately  steep  (40%),  northeast  exposure,  concave  upper  slope 

Drainage:  excessive 

Disturbance:  fire,  wind  throw 

Horizons 

13-11  cm     Organic  horizon;  no  roots 

1 1-0  cm     Organic  horizon;  medium  and  fine  roots  common 

0-9  cm  Gray  brown  (lOYR  5/2)  silt  loam,  5%  gravel;  moderate,  medium 
subangular  blocky  structure;  slightly  hard,  friable,  slightly  sticky,  slightly  plas- 
tic; extremely  acid  (pH  4.0);  few  medium  and  coarse  roots;  diffuse  boundary; 
charcoal  throughout 

9-19  cm  Brown  to  dark  brown  (7.5YR  4/2)  sandy  silt  loam,  5%  gravel;  mod- 
erate, medium  subangular  blocky  structure;  slightly  hard,  friable,  slightly  sticky, 
slightly  plastic;  extremely  acid  (pH  4.0);  few  fine  and  medium  roots;  clear 
smooth  boundary;  charcoal  throughout 

19-35  cm  Brown  to  dark  brown  (7.5 YR  4/4)  sandy  loam,  12%  clay,  20% 
gravel,  5%  cobbles,  reddish  yellow  (7.5YR  6/8)  mottles;  weak,  fine  subangular 
blocky  structure;  slighdy  hard,  friable,  nonsticky,  nonplastic;  very  strongly 
acid  (pH  4.4);  few  fine  and  medium  roots;  diffuse  boundary 

35-40  cm  Dark  reddish  brown  (5YR  3/2)  sandy  loam,  10%  clay;  weak,  fine 
granular  structure;  slightly  hard,  friable,  nonsticky,  nonplastic;  very  strongly 
acid  (pH  4.5);  few  fine  roots;  abrupt  boundary 
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40-58  cm  Brown  to  dark  brown  (7.5YR  4/4)  sandy  loam,  20%  gravel,  10% 
cobble,  30%  stone;  weak,  fine  granular  structure;  slightly  hard,  friable, 
nonsticky,  nonplastic;  pH  unknown;  no  roots;  clear,  smooth  boundary 

58+  cm  Unconsolidated  bedrock 


National  Park  Service  Monograph  2 


115 


Appendix  B:  Vegetation  in 
Permanent  Sampling  Plots  in  the  West 
Twin  Creek  and  Hoh  Lake  Watersheds. 


West  Twin  Creek  Watershed 


Scientific  name  and  authority 


Common  name 


Alpha 


Trees 

Abies  amabilis  (Dougl.)  Forbes 

Acer  macwphyllmn  Pursh 

Alnus  rubra  Bong. 

Picea  sitchensis  (Bong.)  Carr 

Pseudotsiiga  menziesii  (Mirb.)  Franco 

Thuja  plicata  Donn 

Tsuga  heterophylla  (Raf.)  Sarg. 

Shrubs 

Acer  circinatum  Pursh 
Berberis  nen'osa  Pursh 
Gaultheria  shallon  Pursh 
Rubus  spectabilis  Pursh 
Sambucus  racemosa  (T.  &  G.)  Gray 
Vaccinium  alaskaense  How. 
V.  ovalifolium  Smith 
Vparvifolium  Smith 

Herbs  and  Forbs 

Achlys  triphylla  (Smith)  DC. 
Adenocaulon  bicolor  Hook. 
Athyrium  filix-femina  (L.)  Roth 
Blechnum  spicant  (L.)  With. 
Boykinia  elata  (Nutt.)  Greene 
Disporum  smithii  (Hook.)  Piper 
Epilobium  angustifolium  L. 
Galium  triflorum  Michx. 
Goodyera  oblongifolia  Raf. 
Hieracium  albiflorum  Hook. 
Lactuca  muralis  (L.)  Fresen. 
Li  St  era  spp. 

Luzula  parxnflora  (Ehrh.)  Desv. 
Montia  sibirica  (L.)  How. 


Pacific  silver  fir 

ABAM 

bigleaf  maple 

ACMA 

red  alder 

ALRU 

Sitka  spruce 

PISI 

Douglas-fir 

PSME 

western  redcedar 

THPL 

western  hemlock 

TSHE 

vine  maple 

ACCI 

Oregon  grape 

BENE 

salal 

GASH 

salmonberry 

RUSP 

red  elderberry 

SARA 

Alaska  huckleberry 

VAAL 

oval-leaf  huckleberry 

VAOV 

red  huckleberry 

VAPA 

deerfoot  vanillaleaf 

ACTR 

pathfinder 

ADBI 

ladyfem 

ATFI 

deerfem 

BLSP 

coast  boykinia 

BOEL 

Smith's  fairybells 

DISM 

fireweed 

EPAN 

sweetscented  bedstraw 

GATR 

rattlesnake  plantain 

GOOB 

white  hawkweed 

HIAL 

wall  lettuce 

LAMU 

twayblade 

LISSP 

small-flowered  wood-rush 

LUPA 

western  springbeauty 

MOSI 
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Scientific  name  and  authority 


Common  name 


Alpha 


Oxalis  oregana  Nutt.  ex  T.  &  G. 
Polystichum  munitum  (Kaulf.)  Presl. 
Smilacina  racemosa  Wats. 
Tiarella  trifoliata  L. 
Trillium  ovatum  Pursh 

Epiphytes 

Alectoria  sarmentosa 
Dicranum  scoparium  Hedw. 
Eurhynchium  oreganum  (Sull.) 

J.  &S. 
Hylocomium  splendens  (Hedw.) 

B.S.G. 
Lobaria  oregana 
Mnium  spp. 

Polytrichum  juniperinum  Hedw. 
Rhytidiadelphus  loreus  (Hedw.) 

Wamst. 


Oregon  oxalis 

OXOR 

swordfem 

POMU 

false  solomonplume 

SMRA 

three-leaved  coolwort 

TITR 

white  trillium 

TROV 

common  witch's  hair 

ALSA 

broom  moss 

DISC 

Oregon  beaked  moss 


EUOR 


mountain  fern  moss 

HYSP 

lungwort 

LOOR 

MNSPP 

juniper  haircap  moss 

POJU 

little  shaggy  moss 

RHLO 

Hoh  Lake  Watershed 


Scientific  name  and  authority 


Common  name 


Alpha 


Trees 


Abies  amabilis  (Dougl.)  Forbes 

Pacific  silver  fir 

ABAM 

Chamaecyparis  nootkatensis 

(D.  Don)  Spach 

Alaska  yellow  cedar 

CHNO 

Tsuga  mertensiana  (Bong.)  Can* 

mountain  hemlock 

TSME 

Shrubs 

Gaultheria  ovatifolia  Gray 

slender  gaultheria 

GAOV 

Menziesia  ferruginea  Smith 

rustyleaf 

MEFE 

Phyllodoce  empetriformis 

(S.W.)  D.  Don 

pink  mountain  heather 

PHEM 

Rhododendron  albiflorum  Hook. 

white  rhododendron 

RHAL 

Rubus  pedatus  J.E.  Smith 

strawberry-leaf  blackberry 

RUPE 

Sorbus  sitchensis  (Wenzig) 

C.L.  Hitchcock 

Sitka  mountain  ash 

SOSI 

Vaccinium  alaskaense  How. 

Alaska  huckleberry 

VAAL 

V.  deliciosum  Piper 

blueleaf  huckleberry 

VADE 

V^  membranaceum  Dougl.  ex  Hook. 

big  huckleberry 

VAME 

V  ovalifolium  Smith 

oval-leaf  huckleberry 

VAOV 
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Scientific  name  and  authority 

Common  name 

Alpha 

Herbs  and  Forbs 

Erythronium  montanum  Wats. 

avalanche  fawnlily 

ERMO 

Goodyera  oblongifolia  Raf. 

rattlesnake  plantain 

GOOB 

List  era  spp. 

twayblade 

LISSP 

Luzula  parviflora  (Ehrh.)  Desv. 

small-flowered  wood-rush 

LUPA 

Montia  sibirica  (L.)  How. 

western  springbeauty 

MOSI 

Nothochelone  nemorosa  (Dougl.) 

woodland  bear-tongue 

NONE 

Polystichum  munitum  (Kaulf.)  Presl. 

swordfem 

POMU 

Pyrola  secunda  L. 

one-sided  wintergreen 

PYSE 

Pyrola  uniflora  L. 

wax-flower 

PYUN 

Tiarella  trifoliata  L. 

three-leaved  coolwort 

TITR 

Trillium  ovatum  Pursh 

white  trillium 

TROV 

Valeriana  sitchensis  Bong. 

Sitka  valerian 

VASI 

Xerophyllum  tenax  (Pursh)  Nutt. 

bear-grass 

XETE 

Epiphytes 

Alectoria  sarmentosa 

common  witch's  hair 

ALSA 

Die  ran  um  seoparium  Hedw. 

broom  moss 

DISC 

Peltigera  aphthosa  (L.)  Willd. 

PEAP 

Rhytidiopsis  robusta  (Hook.)  Broth. 

robust  moss 

RHRO 

u; 
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Appendix  C.  Estimated  percent  cover  of 

major  vegetation,  rock,  bare  soil,  and 

litter  in  the  permanent  sampling  plots  in 

the  West  Twin  Creek  and 

Hoh  Lake  watersheds, 

Olympic  National  Park,  Washington. 

West  Twin  Creek  Watershed 


Lower  watershed 

Plotl 

Plot  2 

Plot  5 

Tsuga  heterophylla/ 
Polystichum  munitum- 
Oxalis  oregana 
(south  slope) 

Tsuga  heterophylla/ 
Polystichum  munitum- 
Oxalis  oregana 
(north  slope) 

Upper  watershed 

Tsuga  heterophylla/ 
Gaultheria  shallon 

Plots 

Plot  4 

Plot  6 

Tsuga  heterophylla- 
Abies  amabilis/ 
Vaccinium  alaskaense 

Tsuga  heterophylla/ 
Oxalis  oregana 

Tsuga  heterophylla- 
Abies  amabilis/ 
Depauperate 

Lower  watershed 

Upper  watershed 

Lifeform 

Plot  1     Plot  2      Plot  5 

Plot  3     Plot  4    Plot  6 

Rock  (>  8  cm)  3  35  0            10 

Bare  soil  3  10  0                         8 

Litter  30  45  30  60  80 
Tall  shrubs 

Alaska  huckleberry  3  2  7           15 

(Vaccinium  alaskaense) 

red  huckleberry  6  7  8  Tr 

{Vaccinium  parvifolium) 

salmonberry  1  Tr  1  Tr 

{Rubus  spectabilis) 

vine  maple  Tr 

{Acer  circinatum) 


0 

0 

90 

Tr 
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Lower  watershed  Upper  watershed 


Lifeform  Plot  1     Plot  2      Plot  5     Plot  3     Plot  4    Plot  6 

red  elderberry  Tr 

{Sambucus  racemosa) 

Low  shrubs 

salal  5  60  Tr 

{Gaultheria  shall  on) 

Oregon  grape  Tr  Tr 

{Berberis  nervosa) 

Herbs 

swordfem  15  30  1 

{Polystichum  munitum) 

deerfern  1  Tr  Tr  4 

{Blechnum  spicant) 

ladyfem  2  1 

{Athyrium  filix-femina) 

western  springbeauty  Tr 

{Montia  sibirica) 

Oregon  oxalis  40  65  70 

{Oxalis  oregana) 

three-leaved  coolwort  Tr 

(Tiarella  trifoliata) 

twayblade  Tr  Tr 

{Listera  spp.) 

Smith's  fairybells  Tr  Tr  Tr  Tr 

{Disporum  smithii) 

fireweed  Tr 

{Epilobium  angustifolium) 

white  hawkweed  Tr  Tr 

{Hieracium  albiflorum) 

white  trillium  Tr  1  Tr 

(Trillum  ovatum) 

starry  solomonplume  Tr 

{Smilacina  racemosa) 

trail  plant  Tr 

(Adenocaulon  bicolor) 
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Lower  watershed 

Upper  watershed 

Lifeform 

Plot  1 

Plot  2 

Plot  5 

Plot  3    Plot  4    Plot  6 

wall  lettuce 

1 

{Latuca  muralis) 

sweetscented  bedstraw 

Tr 

{Galium  triflorum) 

rattlesnake  plantain 

Tr 

Tr 

{Goody era  oblongifolia) 

small-flowered  wood-rush 

Tr 

Tr 

Tr 

{Luzula  parviflora) 

coast  boykinia 

Tr 

{Boykinia  elata) 

deerfoot  vanillaleaf 

Tr 

11 

{Achlys  triphylla) 

Moss  or  lichen 

Oregon  beaked  moss 

8 

5 

3 

—          Tr            5 

{Eurhynchium  oreganum) 

little  shaggy  moss 

25 

5 

—          Tr 

{Rhytidiadelphus  loreus) 

Mnium  spp. 

2 

1 

mountain  fern  moss 

8 

{Hylocomium  splendens) 

juniper  haircap  moss 

Tr 

{Polytrichum  juniperinum) 

broom  moss 

4 

3 

{Dicranum  scoparium) 

lungwort 
{Lobaria  oregana) 

common  witch's  hair 
{Alectoria  sarmentosa) 


Tr 


—  Tr 


Tr 


Plots  1  and  2  represent  the  same  community  type  but  have  subtle  differences  because  of  aspect. 
Cover  values:  Tr  <1%;  1%-10%  by  1  percent  increments;  15%-100%  by  5  percent  increments. 
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Hoh  Lake  Watershed 
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blueleaf  huckleberry 

{Vaccinium  deliciosum) 

oval-leaf  huckleberry 

{Vaccinium  ovalifolium) 

white  rhododendron 

{Rhododendron  albiflorum) 

rustyleaf 

{Menziesia  ferruginea) 

Sitka  mountain  ash 

{Sorbus  sitchensis) 

Low  shrubs 

slender  gaultheria 
{Gaultheria  ovatifolia) 
pink  mountain  heather 
{Phyllodoce  empetriformis) 
strawberry-leaf  blackberry 
{Rubus  pedatus) 
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Lifeform  Plot  1      Plot  2       Plot  3       Plot  4       Plot  5 

Beargrass  Tr  1 

{Xerophyllum  tenax) 

avalanche  fawnlily  2  2  7  20  3 

(Erythronium  montanum) 

western  springbeauty  Tr 

(Montia  sibirica) 

twayblade  Tr  1 

{Listera  spp.) 

rattlesnake  plantain  Tr 

(Goodyera  oblongifolia) 

three-leaved  coolwort  Tr 

(Tiarella  trifoliata) 

one-sided  wintergreen  Tr 

(Pyrola  secunda) 

wax-flower  Tr 

{Pyrola  uniflora) 

Sitka  valerian  Tr 

{Valeriana  sitchensis) 

wood-rush  Tr 

{Luzula  spp.) 

white  trillium  Tr 

{Trillum  ovatum) 

woodland  bear-tongue  Tr 

{Nothochelone  nemorosa) 

Moss  or  lichen 

common  witch's  hair  3  25  8  3 

{Alectoria  sarmentosa) 

broom  moss  1  4  3  2 

{Dicranum  scoparium) 

freckle  pelt  Tr 

{Peltigera  aphtosa) 

robust  moss  Tr 

{Rhytidiopsus  robusta) 

Plots  1  and  2  represent  the  same  community  type  but  have  subtle  differences  because  of  aspect. 
Plots  3  and  5  represent  the  same  community  type  but  have  subtle  differences  because  of  aspect. 
Cover  values:  Tr  <1%;  1%-10%  by  1  percent  increments;  15%-100%  by  5  percent  increments. 
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Appendix  D.  Biomass  distribution 

(Mg/ha) «  by  plot  in  the  West  Twin  Creek 

and  Hoh  Lake  watersheds, 

Olympic  National  Park,  Washington. 
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As  the  nation's  principal  conservation  agency,  the  Department  of  the  Interior 
has  responsibility  for  most  nationally  owned  public  lands  and  natural  resources. 
This  includes  fostering  sound  use  of  the  land  and  water  resources;  protecting 
fishes,  wildlife,  and  biological  diversity;  preserving  the  environmental  and 
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all  people  by  encouraging  stewardship  and  citizen  participation  in  their  care.  The 
department  also  has  a  major  responsibility  for  American  Native  reservation 
communities  and  for  people  who  live  in  island  territories  under  U.S.  administration. 
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